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Units and Conventions 



In this thesis we study the Universe starting from its early phase and following its 
expansion till cosmological scales. Thus, discussions will span a wide variety of length 
scales, ranging from the Planck length up to the size of the Universe as a whole. 

In any branch of physics the chosen units and conventions are almost always dic- 
tated by the problem at hand. It is not hard to imagine then that also here disparate 
units come into the game. Astronomers will not always feel comfortable with those 
units employed by particle physicists, and viceversa, and so it is worthwhile to spend 
some words in order to fix notation. 

We will employ the so-called natural units, namely h = c = ks = 1, unless otherwise 
indicated, h is Planck's constant divided by 2tt, c is the speed of light, and is 
Boltzmann's constant. Thus, all dimensions can be expressed in terms of one energy- 
unit, which is usually chosen as GeV = lO^eV, and so 

[Energy] = [Mass] = [Temperature] = [Length]"^ = [Tirne]"^ . 

Some conversion factors that will be useful in what follows are 
1 GeV = 1.60x10-3 erg (Energy) 

1 GeV = 1.16x10^3 K (Temperature) 

1 GeV = 1.78x10-2^ g (Mass) 

1 GeV"^ = 1.97x10"^^ cm (Length) 

1 GeV-^ = 6.58x10-25 s (Time) 

The Planck mass mp and associated quantities are given (both in natural and cgs 
units) by 

mp = 1.22 X lO^^GeV = 2.18 X lO^^g 

Ip = 8.2 X 10-2°GeV-^ = 1.62 x lO'^^cm 

tp = 5.39 X 10-^^s 

In these units Newton's constant is given hj G = 6.67 x 10~^cm3g~^sec~2 = rrip^. 

These fundamental units will be replaced by other, more suitable, ones when study- 
ing issues on large-scale structure formation. In that case it is better to employ as- 
tronomical units, which for historical reasons are based on solar system quantities. 



Thus, we use the parsec (the distance at which the Earth-Sun mean distance subtends 
one second of arc), or better the megaparsec, which is more useful when deahng with 
structures on cosmological scales. Regarding masses, the standard unit is given by the 
solar mass Mq. 

Ipc ^ 3.26 light years = 3.1 x lO^^cm 
IMq = 1.99 X lO^^g 
IMpc = 3.1 X lO^^'cm 

It is worthwhile also recalling some cosmological parameters that will be used in 
this thesis, like the present Hubble time (Hq^) and distance (cHq^), and the critical 
density pc- 

Hq^ = 3.09 X W^^h-^ s 
cHq ^ = 2997.9/1-^ Mpc 

Pc = 3H^/{8ttG) = l.SS/i^ X 10"^^ g cm-3 = l.OS/i^ x 10^ eV cm'^ 

We adopt the following conventions: Greek letters denote spacetime indices and 
run from to 3; spatial indices run from 1 to 3 and are given from Latin letters. 

For any given quantity, the subscript zero (o) obviously indicates the present value, 
while the prime indicates that it is referred to the mirror sector instead of the ordinary 
one (example: D,b is the density parameter of ordinary baryons, is the one of mirror 
baryons) . 

Below we list some symbols and abbreviations used in the text. 





total density parameter 




density parameter of vacuum 


— 


density parameter of matter 


fib = 


density parameter of baryons 


UJi, = 


nb 




= 


ordinary 


M = 


mirror 


CMB 




cosmic microwave background 


LSS 




large scale structure 


CDM 




cold dark matter 


WDM 




warm dark matter 


HDM 




hot dark matter 


SIDM 




self-interacting dark matter 


MBDM 




mirror baryon dark matter 



FRW 

L IVVV 


nT*! f^H TYl ClTlTl — r? r\f»T*'l"G;r^n — ^\A/C» 1 k^^T* 

± i ICLlllidilil lL*JUd totJil VVCliA.d 


RA 


— hJcXiL yyjii cxo y liiiiidji y 


BBN 


— Rin R/^np" nnplpocivnf np^ic; 


GUT 


— Errand unififtci thftorv 


EW 


= electroweak 


SM 


= standard model 


MSSM 


= minimal supersymmetric standard model 


RD 


= radiation dominated 


MD 


= matter dominated 


MRE 


= matter-radiation equality 


MRD 


= matter-radiation decoupling 



Chapter 1 

Standard cosmological paradigm 



1.1 The expanding Universe 

On large scales the Universe is homogeneous and isotropic (to a good approximation), 
then it does not possess any privileged positions or directions: this idea is stated as the 
Cosmological Principle (also known as the Copernican Principle: no observer occupies 
a preferred position in the Universe).^ In effect, we know only the observable Universe 
(our present Hubble volume) , but for purposes of description we may assume the entire 
Universe is homogeneous and isotropic. 

The most general space-time metric describing a Universe in which the cosmological 
principle is obeyed is the maximally-symmetric Friedmann- Robertson- Walker metric 
(FRW), written as 



dr 



2 



1 — kr'^ 



where (t, r, 9, 0) are comoving coordinates^, a{t) is the cosmic scale factor, and, 
with an appropriate rescaling of the coordinates, k can be chosen to be +1, —1, or 
for spaces of constant positive, negative, or zero spatial curvature, respectively. The 
coordinate r is dimensionless, while a{t) has dimensions of length. 



"'^ There is a stronger formulation called the Perfect Cosmological Principle, in which the Universe 
is the same not only in space, but also in time; this idea led to develop the Steady State cosmology, 
theory that implies the continuous creation of matter to keep the density of the expanding Universe 
constant. This theory was abandoned for the difficulties encountered to explain the properties of the 
cosmic microwave background, radio sources and the primordial nucleosyntesis. 

^An important concept is that of the "comoving observer" . Loosely speaking, a comoving observer 
follows the expansion of the Universe including the effects of any inhomogeneities that may be present. 
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The time coordinate in eq. (1.1) is just the proper (or clock) time measured by an 
observer at rest in the comoving frame, i.e. (r, 9, (f)) — const.; this coordinate system is 
called the synchronous gauge (see appendix A. 3). 

It is often convenient to express the metric in terms of conformal time r, defined 
by dr — dt/a{t) 



ds^ = a^(r) 



Redshift and Hubble law 



1 — kr^ 



;i.2) 



Consider the equation of geodesic motion (A. 4) of a particle that is not necessarily 
massless. The four-velocity it" of a particle with respect to the comoving frame is 
referred to as its pecuhar velocity. If we choose the a — component of the geodesic 
equation, in the particular case of the FRW metric, we find that u oc a~^, and the 
magnitude of the momentum of a freely-propagating particle "red shifts" as: |p| oc a^^. 

In this way we can see why the comoving frame is the natural frame. Consider an 
observer initially (at time t) moving non-relativistically with respect to the comoving 
frame with physical three velocity of magnitude v. At a. later time to, the magnitude 
of the observer's physical three velocity vq, will be 

Then, in an expanding Universe the free-falling observer is destined to come to rest 
in the comoving frame even if he has some initial velocity with respect to it (in this 
respect the term comoving is well chosen) . 

The light emitted by a distant object can be viewed quantum mechanically as 
freely-propagating photons. Since the wavelength of a photon is inversely proportional 
to its momentum (A = h/p), that changes in proportion to a"^, the wavelength at time 
to, denoted as Aq, will differ from that at time t, denoted as A, by 

A = ^ . (1.4) 
Ao a{to) 

As the Universe expands, the wavelength of a freely-propagating photon increases, just 
as all physical distances increase with the expansion. 

Hence, we introduce a new variable related to the scale factor a which is more 
directly observable: the redshift of an object, z, defined in terms of the ratio of the 
detected wavelength Aq to the emitted wavelength A 

z^^^ l + z=^. (1.5) 

A a[t) 
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Any increase (decrease) in a{t) leads to a red shift (blue shift) of the hght from distant 
sources. Since today astronomers observe distant galaxies to have red shifted spectra, 
we can conclude that the Universe is expanding. 

Hubble's law, the linear relationship between the distance to an object and its 
observed red shift, can be expressed (assuming negligible peculiar velocity) as 

z^Hd^ IQ-^^cm-^ X d , (1.6) 

where d is the proper distance ^ of a source, and H is the Hubble constant or, more 
accurately, the Hubble parameter (because it is not constant in time, and in general 
varies as t~^), defined by 

H = ^. (1.7) 

ait) ^ > 

Now, we can expand a{t) in a Taylor series about the present epoch (i.e., for times 
t close to to): 



a{t) = a(to) 



1 + Ho{t - to) - ^ qoHlit - toY + ... 



;i.8) 



where 



a{to) a^[to) a[to)H^ 



are the present values of the Hubble constant and the so-called deceleration parameter. 

At present time the Hubble parameter Hq is not known with great accuracy, so it 
is indicated by 

Ho = 100 h km s-^Mpc"^ 0.5 < h < 0.8 , (1.10) 

where the dimensionless parameter h contains all our ignorance on Ho- 

The present age and the local spatial scale for the Universe are set by the Hubble 
time^ and radius 



Hq^ ~ 9.778 X lO^h-^ yr ~ 3000 h'^ Mpc . (1.11) 



^The proper distance of a point P from another point Pq is the distance measured by a chain of 
rulers held by observers which connect P to Pq- 

*Note that earlier than some time, say tx, or better for a less than some ax^ our knowledge of the 
Universe is uncertain, so that the time elapsed from a = to a = ax cannot be reliably calculated. 
However, this contribution to the age of the Universe is very small, and most of the time elapsed since 
a = accumulated during the most recent few Hubble times. 
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The Friedmann equations and the equation of state 



GafS = Rap — T^RgaP — —^Taf) + h-QaP , (1-12) 



-) +-,=-Gp. (1.15) 



The expansion of the Universe is determined by the Einstein equations 

1 SttG, 
-Rg^p - — 

where R^is is the Ricci tensor, R is the Ricci scalar, g^p is the metric tensor, T^p is the 
energy-momentum tensor, and A is the cosmological constant^ For the FRW metric 
(1.1), they are reduced to the form (see appendix A.l) 

- = -^G(p + 3p) (1.13) 
a 6 

for the time-time component, and 

^ + 2f^)' + 24 = 47rG(p-p) (1.14) 
a \aj 

for the space-space components, where, if the cosmological constant is present, p and 
p are modified according to eq. (A. 9). From eqs. (1-13) and (1-14) we obtain also 

'a\^ A; _ Stt 
,a/ 3 

The equations (1.13) and (1.15) arc the Friedmann equations; they are not inde- 
pendent, because the second can be recovered from the first if one takes the adiabatic 
expansion of the Universe into account. 

From eq. (1.13), models of the Universe made from fluids with — 1/3 < w < 1 have 
d always negative; then, because today a > 0, they possess a point in time where a 
vanishes and the density diverges; this instant is called the Big Bang singularity. Note 
that the expansion of the Universe described in the Big Bang model is not due in any 
way to the effect of pressure, which always acts to decelerate the expansion, but is a 
result of the initial conditions describing a homogeneous and isotropic Universe. 

The Friedmann equation (1.15) can be recast as 

where fl is the ratio of the density to the critical density pc necessary for closing the 
Universe ^ 



^Thc cosmological constant A was introduced by Einstein for the need to have a static Universe. 
However, now we know that the Universe is expanding for sure. In fact, one of the biggest theoretical 
problems of the modern physics is to explain why the cosmological term is small and not order M|>. 

^The present value of the critical density is poc ~ 1-88 h"^ x 10^^" g cm^'', and taking into account 
the range of permitted values for h, this is ~ (3— 12) x 10~^^ kg m^'', which in either case corresponds 
to a few H atoms/m^. Just to compare, a 'really good' vacuum (in the laboratory) of 10~® N/m^ at 
300 K contains ~ 2 x 10^^ molecules/m^. The Universe seems to be empty indeed! 
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Since iJ^a^ > 0, there is a correspondence between the sign of k, and the sign of (Jl — 1) 

CLOSED A; = +1 =^ Q > 1 

FLAT A; = =^ Q = 1 

OPEN A; = -1 =^ < 1 



Prom equation (1.16) we find for the scale factor today 

3000 h-^ Mpc 
I Qo - 1|V2 




(1.18) 



which can be interpreted as the current radius of curvature of the Universe. If = 1, 
then oo has no physical meaning and can be chosen arbitrarily (it will always cancel 
out when some physical quantity is computed). 

In order to derive the dynamical evolution of the scale factor a{t), it is necessary 
to specify the equation of state for the fluid, p — p{p). It is standard to assume the 
form 

p — wp (1-19) 

and consider different types of components by choosing different values for w (assumed 
to be constant). 

If the Universe is filled with pressureless non-relativistic matter (dust), we are in 
the case where p <^ p and thus w = 0. Instead, for radiation, the ideal rclativistic 
gas equation of state p = l/3p is the most adapt, and therefore w = 1/3. Another 
interesting equation of state is p = —p, corresponding to w = —1. This is the case 
of vacuum energy ^, which will be the relevant form of energy during the so-called 
inflationary epoch; we will give a brief review of inflation in §1.3. 

The a = component of the conservation equation for the energy-momentum 
tensor, T°'^p = 0, gives the 1st law of thermodynamics 

p + 3H{p + p) ^0 , (1.20) 

where the second term corresponds to the dilution of p due to the expansion {H — a/a) 
and the third stands for the work done by the pressure of the fluid. 



''There is also a widely popular idea that the "dark energy" can be variable in time, related to 
the energy density of very slowly rolling scalar field (so-called "quintessence" ) . This situation would 
imply w in general different from -1, and it can be tested by high precision data on CMB and matter 
power spectra (see e.g refs. [8, 10]). However, in this thesis we do not consider this possibility and in 
our calculations assume the dark energy to be constant in time. 



6 



Chapter 1. Standard cosmological paradigm 



Prom equations (1.19) and (1.20) we can obtain the relation 

^q3(1+«;) ^ ^^^^^ _ (^^ 21) 

In particular we have, for a dust or matter Universe 

w = — > pma^ = const — > pm oc pom(l + zf , (1-22) 
for a radiative Universe ^ 

W — 1/3 )• PrO^ — const >■ Pr OC Por(l + ^)^ : (1-23) 

and for a vacuum Universe 

w — —1 — > px = Pox = const . (1-24) 

Recalling the definition of the deceleration parameter q from eq. (1.9) and using 
eqs. (1.13) and (1.15), we may express it as follows 

? = ^Q+^)^^- (1-25) 

This shows that w — —1/3 is a critical value, separating qualitatively different models. 
A period of evolution such that w < —1/3 {q < 0) is called inRationary. For the 
present value go; we find for a matter dominated (MD) model q'q = ^o/2, for a radiation 
dominated (RD) model go — ^o, ^-nd for a vacuum dominated model go — (the 
expansion is accelerating, a > 0). 

If we want to calculate the density parameter Q at an arbitrary redshift 2; as a 
function of the present density parameter Jlo, we can use the expressions (1.17) together 
with the equations (1.21) and (1.15) to obtain the relation 

Note that if > 1, then Q{z) > 1 for all z; likewise, if < 1, then il{z) < 1 for 
all z; on the other hand, if Qo = 1, then il{z) — 1 for all time. The reason for this is 
clear: the expansion cannot change the sign of the curvature parameter k. It is also 
worth noting that, as 2; ^ 00, i.e., as we move closer and closer to the Big Bang, fl{z) 
always tends towards unity: any Universe with 7^ 1 behaves like a flat model in the 
vicinity of the Big Bang. We shall come back to this later when we discuss the flatness 
problem in §1.3. 



^Consider for example photons: not only their density diminishes due to the growth of the volume 
(oc a~^), but the expansion also stretches their wavelength out, which corresponds to lowering their 
frequency, i.e., they redshift (hence the additional factor a a~^). 
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1.2 Flat models 



From equations (1-15) and (1-21) we can obtain for /c = 

l+3w 



-j) + (1 - no) 



(1.27) 



a(t)=ao(-j , (1.28) 



Now we shall find the solution to this equation appropriate to a flat Universe. For 
Q,Q — 1, integrating equation (1.27) one can obtain 

^ ^ 2/3{l+w) 
to, 

which shows that the expansion of a flat Universe lasts an indeflnite time into the 
future; equation (1.28) is equivalent to the relation between cosmic time t and redshift 

t = to(l + -2)~'^'^"^^' ■ (1-29) 
From equations (1.28), (1-29) and (1-21), we can derive 

H ^ ' =Ho'-^ = Ho{l + zr^-^/' , (1.30) 

a 6[1 + w)t t 

_ ad l + 3w 
q = -— = — - — = const. = qo , (1.31) 



2 



3(1 +7^)iJo 



(1.32) 



In appendix A. 2 we report the above relations for the special cases of a Universe 
dominated by matter or radiation. A general property of flat Universe models is that 
the scale factor a grows indefinitely with time, with constant deceleration parameter 
go- The role of pressure can be illustrated again by the fact that increasing w and, 
therefore, increasing the pressure causes the deceleration parameter also to increase. 

A cosmological model in which the Universe is empty of matter and has a positive 
cosmological constant is called the de Sitter Universe. From equations (A. 9) and (1-15) 
we obtain 

'd\^ StiG 

- =^PA, 1.34 
a/ 3 

which for positive A implies the exponentially fast expansion 
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corresponding to a Hubble parameter constant in time. In the de Sitter vacuum Uni- 
verse test particles move away from each other because of the repulsive gravitational 
effect of the positive cosmological constant. 

Finally, the age of flat Universe that contains both matter and positive vacuum 
energy (fi = + 0,a — 1) is 



In 



1 + Q 



1/2 



1/2 



(1.36) 



It is interesting to note that, unlike previous models, a model Universe with Q\ > 
0.74 is older than Hq^; this occurs because the expansion rate is accelerating. In the 
limit Q\ — s> 1, to oo. For this reason the problem of reconciling a young expansion 
age with other independent age determinations (like for example the globular clusters) 
has at several times led cosmologists to invoke a cosmological constant. 



1.3 Inflation 

The inflation is a microphysical mechanism which operates at very early times {t ~ 
10"^"^ s, well above the Planck time^), according to which the Universe, at that stage 
dominated by vacuum energy and filled with the potential energy of a scalar field 
(called inHaton). underwent a brief epoch of exponential expansion of the type (1.35). 
This field is initially displaced from the minimum of its potential, so it decays from 
the false vacuum to the true vacuum (lowest energy) state, and the system has to 
tunnel across the bump and then it slowly rolls down the potential. After reaching 
the minimum of the potential, it executes damped oscillations, during which energy 
is thermalised and entropy is increased (reheating). From that point onwards the 
Universe is well described by the hot Big Bang model. 

Due to the exponential expansion, primordial ripples in all forms of matter-energy 
perturbations at that time were enormously amplified and stretched to cosmological 
sizes and, after the action of gravity, became the large-scale structure that we see 
today, as firstly argued by Guth in 1981 [86]. 

In a Universe dominated by a homogeneous scalar field with potential V{4>), 

^Thc Planck time tp ^ 10~'^^s is the time after the Big Bang for which quantum fluctuations are 
no longer negligible and the theory of General Relativity should be modified in order to take account 
of quantum effects. One can then define Planck mass Mp ~ pp Z|, ~ 1.22 x 10^® GeV, where Ip c^^ ctp 
is the Planck length and pp ~ {Gtp)~^ is the Planck density. 
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minimally coupled to gravity, the equation of motion is given by 

dV 

(j> + 3Hcj) + — = . (1.37) 

0(p 

The criterion that the potential must be sufficiently flat is quantified by the slow- 
rolling conditions: 

1 fV'V 

- 16^ y «^ ' ^'-''^ 



1 V 



<1 , (1.39) 



SttG V 

where primes denote differentiation with respect to (f). 

Due to the period of accelerated expansion of the Universe, all those perturbations 
present before inflation were rendered irrelevant for galaxy formation, since inflation 
washes out all initial inhomogeneities; otherwise, the same mechanism stretched the 
quantum mechanical fluctuations enough to produce scalar (density) and tensor (grav- 
itational waves) perturbations on cosmologically interesting scales. 

An important feature of the inflation is the ability to generate power spectrum of 
scalar and tensor perturbations as power laws with spectral index respectively [114, 
105, 170] 

n = ns = l + 2ri-6e Ut = -2e , (1.40) 

where n — 1 represents scale invariance. Thus, we see that, according to the slow- 
rolling conditions, the inflation naturally produces flat spectra, with n — 1 and rit — 0. 

The great success of the inflation is to solve some shortcomings present in Big Bang 
cosmology, namely the flatness, the horizon, and the unwanted relics problems. 

The horizon problem. The CMB radiation is known to be isotropic with a very 
high precision. In fact, two microwave antennas pointed in opposite directions in the 
sky do collect thermal radiation with AT/T < 10~^, T being the black body temper- 
ature. The problem is that these two regions from which the radiation of strikingly 
uniform temperature is emitted cannot have been in causal contact at the time of last 
scattering [148], since the causal horizon at that time subtends an apparent angle of 
order only 2 degrees today. How then causally disconnected spots of the sky got in 
agreement to produce the same temperature anisotropics? 

The inflation can solve this problem: in fact, two regions initially causally connected 
are moved away one from another by the exponential inflationary expansion on the 
Universe itself and exit from their horizons, so that now we see regions of the sky 
disconnected at decoupling, but that were connected before inflation. 



10 



Chapter 1. Standard cosmological paradigm 



The flatness problem. Prom equation (1.26) we learned that the density parame- 
ter Q rapidly evolves away from 1 while the Universe expands. Given that observations 
indicate that Q is very close to 1 today, we conclude that it should have been much 
closer to 1 in the past. Going to the Planck time we get |Q(10~^^s) — 1| ~ C(10~^^), 
and even at the time of nucleosynthesis we get |Q(ls) — 1| ~ O{10~^^). This 
means, for example, that if Q at the Planck time was slightly greater than 1, say 
Q(10~^^s) — 1 + 10~^^, the Universe would have collapsed millions of years ago. How 
the Universe was created with such fine-tuned closeness to Q — 17 

Inflation can avoid this fine-tuning of initial conditions, because, whatever was the 
curvature of the local space-time, the huge inflationary expansion stretches all length, 
included the curvature radius, so that it becomes large enough that at the end of the 
inflation the Universe is as flat as we want. 

The unwanted reUcs problem. Arising from phase transitions in the early 
Universe, there is an overproduction of unwanted relics. As an example, defects (e.g., 
magnetic monopoles) are produced at an abundance of about 1 per horizon volume 
(Kibble mechanism), yielding a present f^monop far in excess of 1, clearly cosmologically 
intolerable. How can the cosmology rid the Universe of these overproduced rehcs? 

Even now, the key is the huge expansion produced by an early inflationary era, 
which dilutes the abundance of (the otherwise overproduced) magnetic monopoles or 
any other cosmological "pollutant" relic. ^'^ 

1.4 Thermodynamics of the Universe 

In this section we will study the properties of the Universe considered as a thermody- 
namic system composed by different species (electrons, photons, neutrinos, nucleons, 
etc.) which, in the early phases, were to a good approximation in thermodynamic 
equilibrium, established through rapid interactions. Obviously, coming back to the 
past, decreasing the cosmic scale factor we have an increase of the temperature. 

Let us evaluate the total energy density p and pressure p of the cosmological fluid, 
considering all the particles in thermal equilibrium in the Universe, and express them 
in terms of photon temperature T^ = T using the eqs. (A. 35) and (A. 36) 



This also implies, however, that any primordial baryon number density will be diluted away 
and therefore a sufficient temperature from reheating as well as baryon-number and CP-violating 
interactions will be required after inflation (see Kolb & Turner, 1990 [104]). 



P 




(1.41) 
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^ - ' ^Iry' GTT^i, exp{«,-y,}±l ' ^'-^'^ 

where the sums run over all species i (which may have a thermal distribution, but with a 
different temperature than that of the photons) and Ui — Ei/Ti, Xi — rrii/Ti, i/i — /li/Ti. 
Note that, from eqs. (A. 37) and (A. 38), non relativistic particles contribute neghgibly 
to the energy density in the radiation dominated era, since their energy density is 
exponentially suppressed with respect to the case of relativistic particles; thus we can 
neglect the contribution of non relativistic species in the sums above. Assuming all 
species non degenerate, we then get 

P - , (1-43) 

P ^ Ip - "^g^T' , (1.44) 
where in contribute only the relativistic degrees of freedom 

B 8 p 



^* ^ E + 8 ^ [ij ^^-^^^ 



{B — bosons, F — fermions). Note that g^^ is in general a function of T, since the 
number of degrees of freedom becoming relativistic at a given temperature depends on 
T itself. Moreover, at a given time, not all (relativistic) particles in the bath are in 
equilibrium at a common temperature T}^ A particle will be in kinetic equilibrium 
with the background plasma (that is Tj = T) only as long as its interaction with the 
plasma is fast enough; although the conditions for this to occur will be discussed in 
§1.5, it is obvious that these involve a comparison between the particle interaction and 
the expansion rate H. 

Using definition (A.47) and eqs. (1.43), (1.44) we can evaluate the total entropy 
density in the comoving volume, obtaining 

■9*sT^, (1.46) 



45 



where 



^^Some example: for T <^ 1 MeV (the only relativistic species are the three neutrino species and 
the photon) = 3.36; for T > 1 McV, up to 100 MeV, = 10.75 (the electron and positron are 
additional relativistic degrees of freedom); for T > 100 GeV (all the species in the standard model - 
8 gluons, W^, Z, 7, 3 generations of quarks and leptons, and 1 complex Higgs doublet - should have 
been relativistic) gr* = 106.75. 
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and the sums again run only over the relativistic degrees of freedom in equihbrium 
(in the considered approximation) . For most of the history of the Universe all particle 
species had a common temperature, and can be replaced by g^^g- Note that s in eq. 
(1.46) can be parametrized in terms of the photon number density, using the eq. (A. 41), 
as follows 

s = 45((3) ~ l.SOg^snj, (1.48) 

where ({x) is the Ricmann zcta-function and ({3} ~ 1.2021. From entropy conservation 
we can now obtain the scaling law relating the cosmic scale factor with the temperature; 
in fact, from eqs. (A. 46) and (1.46) we get 

T ~ g-'^'a-'. (1.49) 

We stress that only if g^,s is constant we can obtain the familiar result T oc a~^, valid 
for pure expansion. 

Time - temperature relationship 

The useful relationship between the time t and the background (photon) temperature 
T in the radiation dominated era can be obtained straightforwardly by integrating the 
Friedmann equations by means of entropy conservation. In fact, we have 

r*^*^ da 

t ^ --. 1.50 

JQ Ha 

By using the equations (1-15) (neglecting the curvature term in the radiation dominated 
era) and (1.43), considering that G ~ Mp^, we obtain 



from which, considering the entropy conservation, and during the periods when both 
g^ and ^f^^ are approximately constant (i.e. away from phase transitions and mass 
thresholds where relativistic degrees of freedom change), we find 

t ~ 0.301 g:^/"^ ^ — > t{sec) ~ ^-^(MeV) . (1.52) 



1.5 Thermal evolution and nucleosynthesis 

Sufficiently away from the Big Bang event, we can approximate the evolution of the 
Universe as made of several subsequent phases of (different) thermal equilibrium with 
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temperature T mainly varying as a~^. Thermal equilibrium is realized if the reactions 
between the particles in the heat bath take place very rapidly compared with the 
expansion rate, set by H; so the key to understanding the thermal history of the 
Universe is the comparison of the particle interaction rates and the expansion rate. 

Denoting with F — n{av) the thermal average of the given scattering reaction rate 
{n being the number density of target particles, a and v the cross section and the 
particle velocity, respectively), the (approximate) condition for having equilibrium is 
r ^ H, and the decoupling temperature Tn is defined by 

r(r^) = H{Tn) . (1.53) 

If we ignore the temperature variation of g^,, in the RD era the expansion rate is 
approximately, from eq. (1-51) 

while the scattering rate depends on the particular interactions experienced by the 
particles. 



Neutrino decoupling 

Massless neutrinos are maintained in equilibrium by reactions such as vV e'^e", 
ve <-> ue and so on, whose reaction rates are given by 

r - GlT^ , (1.55) 

where Gp — 1-17 x 10~^ GeV~^ is the Fermi coupling constant. From eqs. (1.53), 
(1.54), (1.55) we then find the neutrino decoupling temperature 

~ IMeV. (1.56) 

Thus, above ~ 1 MeV neutrinos are in equilibrium with the plasma of photons with 
— T and from eq. (A. 41) n„ — (3/4)n^, afterwards they decouple from the plasma 
and their temperature T^, scales approximately as a~^. Subsequently, as the temper- 
ature drops below ~ 0.5 MeV, e+ and e~ annihilate, heating the photons but not the 
decoupled neutrinos. Thus, from eq. (1.49) with — g^g we find that g*{aT)^ — const, 
and so across annihilation (aTi,) remains constant, while the entropy transfer 

^^The correct way of proceeding is to solve the Boltzmann transport equations for the given species, 
but here we are interested in a semi-quantitative discussion, and this criteria fit well. 
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increases (aT^) by a factor 



T 

J- V 



g^irrie < T < Tj 



^] - 1-40. (1.57) 



g^{T < me) 

After annihilation there are no other relativistic species that can become non 
relativistic (altering the effective degrees of freedom), so that eq. (1.57) is just the 
relation between the present values of 7 and v temperature (given T = ~ 2.73 K, 

~ 1.96 K). 

Matter - radiation equality 

Remembering the equations (1.22) and (1.23), which relates matter and radiation 
densities to the redshift, we note that they scale with a{t) in different ways, so that 
there must be a moment in the history of the expanding Universe when pressureless 
matter will dominate (if very early on it was radiation the dominant component, as 
is thought to be the case within the hot Big Bang). In fact, PmatI Prad oc a(t). This 
moment, known as the matter-radiation equality time (MRE), is an important time 
scale in the thermal history of the Universe, and is usually denoted as teq 

Priteq) = pm{teq) ■ (1.58) 

Then it follows that the redshift, temperature and time of equal matter and radiation 
energy densities are given by 

l + z,q = ^^a:^'^^^2.32xlO'Q^h' , (1.59) 

Teq = Toil + Zeq) =5.50 n^h'' eV , (1.60) 



^eq 



1.4 X 10='(fio/i')"'/'(^^m/i')"'/' yr . (1.61) 



Photon decoupling and recombination 

In the early Universe matter and radiation were in good thermal contact, because of 
rapid interactions between the photons and electrons. Since the photons continue to 
meet electrons, they do not propagate in straight paths for very long distances; in 
other words, the Universe is opaque to electromagnetic radiation. However, when the 
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electrons become slow enough, they are captured by the protons, forming stable hydro- 
gen atoms. When this process, called recombination is completed, the photons find 
no more free electrons to scatter against, and they do not scatter against the neutral 
hydrogen atoms; thus the matter becomes transparent to radiation, and the decoupled 
photons continue moving along in straight lines (or, more precisely, geodesies) from the 
last scattering surface^^ until today. It is these photons, redshifted by the expansion of 
the Universe, that are known as the cosmic microwave background radiation (CMB) 
(for more details see chapter 3). 

This occurs when the interaction rate of the photons becomes equal to the 
Hubble rate H, with 

r^ = ne(7T, (1.62) 

where Ue is the number density of free electrons, and ar — 6.652 • 10~^^ cm^ is the 
Thomson cross section. 

For simplicity we assume all the baryons in the form of protons. The charge neu- 
trality of the Universe implies Up — Ug, and baryons conservation implies rib — rip + nn- 
In thermal equilibrium, at temperatures less than m^, where i — e,p, H, the equation 
(A. 37) is valid for every species i, and in chemical equilibrium the process p+e — > 
guarantees that /ip + /le — /J'H- By introducing the fractional ionization Xg = Up/rib, 
we can obtain the Saha equation for the equilibrium ionization fraction 



where B = nip + rUe — tuh — 13.6 eV is the binding energy of hydrogen and r) — rih/n^ 
is the baryon to photon ratio. 

If we define recombination as the point when 90% of the electrons have combined 
with protons, we find that it occurred at a redshift Zrec ~ 1100, with temperature 

Tree = To(l + Zrec) ^ 3000 K ~ 0.26 cV . (1.64) 

If the Universe was matter dominated at decoupling, its age was 

tdec ^ tree = ^^o"' ^0 (1 + ^rec)"'/' ^ 180000 (l^o h^^^^ ■ (1-65) 



^■^Thc term recombination is slightly misleading, because the electrons have never been combined 
into atoms before. 

^^The last scattering did not occur to all photons at the same time, so this surface is really a shell 

of thickness Az ~ 0.07z. 

^^Note that recombination occurs at a temperature of about 0.3 eV, not at T ~ S ~ 13.6 eV, 
because the released binding energy reheats the remaining electrons, and the large amount of entropy 
in the Universe favours free photons and electrons. 
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Primordial nucleosynthesis (BBN) 

According to this theory, going back in time, we reach densities and temperatures 
high enough for the synthesis of the Ughtest elements: when the age of Universe was 
between 0.01 sec and 3 minutes and its temperature was around 10 - 0.1 MeV the 
synthesis of D, "^He, ^He, ^Li took place. The nuclear processes lead primarily to 
"^He, with a primordial mass fraction of about 24%, while lesser amounts of the other 
elements are produced: about 10^^ of D and ^He, and about 10~^° of ^Li by number 
relative to H. The prediction of the cosmological abundance of these elements is one of 
the most useful probes of the standard hot Big Bang model (and certainly the earliest 
probe we can attain) [133, 181]. 

The outcome of primordial nucleosynthesis is very sensitive to the baryon to photon 
ratio 1] = rih/n^ and the number of very light particle species, usually quantified 
as the equivalent number of light neutrino species N,y (=3 in the Standard Model 
physics). Then, the nucleosynthesis predictions can be compared with observational 
determinations^'' of the light elements abundances, obtaining 

2.6 X 10-^° < 77 < 6.2 X 10"^° . (1.66) 

This parameter provides a measure of the baryon content of the Universe. In fact, 
fixed rij by the present CMB temperature (from eq. (A.41) n° = (2C(3)/7r^) ~ 
422 cm"^) , we find 

Qbh'^ ~ 3.66 X 10^77 . (1.67) 

The presence of additional neutrino flavors (or of any other relativistic species) at 
Tw ~ 1 MeV (the "freeze-out" temperature of weak interactions) increases g^, hence 
the expansion rate, the value of TV) the neutron to proton ratio, and therefore Ybbn- 
Otherwise, the limits on Nj, can be translated into limits on other types of particles 
that would affect the Hubble expansion rate during BBN. In chapter 4 we will use 
these concepts in the presence of a mirror sector. 



^^Here we describe only few generalities on this phenomenon, referring the reader to §4.5 for a more 
quantitative discussion. 

^''As the ejected remains of stellar nucleosynthesis alter the light element abundances from their 
primordial values, but also produce heavy elements ("metals"), one seeks astrophysical sites with low 
metal abundances, in order to measure light element abundances very close to primordial. 
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1.6 Baryogenesis 

Today the Universe seems to be populated exclusively by matter rather than anti- 
matter. There is in fact strong evidence against primary forms of antimatter in the 
Universe. Furthermore, as seen in the previous section, the density of baryons com- 
pared to the density of photons is extremely small, rj ~ 10~^°. 

It is well known that a non-zero baryon asymmetry (BA) can be produced in the 
initially baryon symmetric Universe if three following conditions are fulfilled: (i) B 
violation, (ii) C and CP violation, and (iii) departure from the thermal equilibrium 
[153]. The first two of these ingredients are expected to be contained in grand unified 
theories (GUT) as well as in the non-perturbative sector of the standard model. 

Generally speaking, the baryogenesis scenarios can be divided in two categories, in 
which the out of equilibrium conditions are provided (a) by the Universe expansion 
itself, or (b) by fast phase transition and bubble nucleation. In particular, the lat- 
ter concerns the electroweak baryogenesis schemes, while the former is typical for a 
GUT type baryogenesis or Icptogcnesis. In §4.3 we will treat in more detail these two 
scenarios in both the standard and mirror sectors. 

At present it is not possible to say which of the known mechanisms is responsible 
for the observed BA. We only know that the baryon to photon number density ratio 
T) — rib/rij is restricted by the BBN constraint given by eq. (1.66). 

1.7 Present status of cosmology 

In the previous sections we gave a brief review of the so-called "standard cosmological 
paradigm" . As we saw, there is some free parameter which describes the state and the 
evolution of the Universe; in this section we try to describe the present status in the 
determination of these parameters. 

First of all, it is worthwhile to remember that during the last years cosmology 
seems to have reached a so-called concordance [172], given that almost all parameters, 
as measured by sometimes very different methods, are converging to the same values. 
This is very important, as it hardly suggests that we are working in the right direction, 
but, at the same time, it confirms the existence of some trouble in the today's picture of 
the Universe, with a strong evidence of the need of new physics or new forms of matter, 
till unknown. In particular, the estimates of the cosmological parameters renews one 
of the most exciting problem in astrophysics, the dark matter (see next chapter). 
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A powerful instrument to estimate all the "dynamical" cosmological parameters is 
the joint analysis of CMB and LSS power spectra (see chapter 3) , as made for example 
by Percival et al. [140, 141], Wang et al. [180], and Sievers et al. [162], all assuming 
that the initial seed fluctuations were adiabatic, Gaussian, and well described by power 
law spectra (for more details see next chapters) . We will study the CMB and LSS for 
a mirror Universe in chapter 6. 

Total density of the Universe. Prom the first results of BOOMERANG experiment 
[49], it seemed immediately clear that the Universe is approximately spatially flat, as 
showed by the CMB anisotropies. Subsequently, other experiments (MAXIMA [111], 
DASI [87, 145], VSA [159], CBI [119, 136], and again BOOMERANG [131]) confirmed 
and are still confirming this result, which, very interestingly, is in agreement with the 
inflationary prediction (see §1.3). Using the result of ref. [162], we have 

Qo = 1-00 ±^:^^ . (1.68) 

Bciryonic density. According to the standard BBN scenario, its value can be ob- 
tained in direct measurements of the abundances of the lights elements (see §1.5), and 
its strongest constraint comes from the primordial deuterium abundance observed in 
Lyman-a feature in the absorption spectra of high redshift quasars, which in a recent 
analysis [38] gives 

u^C^) = 0.020 ± 0.001 , (1.69) 

while, from the cited indirect observations provided by the CMB and LSS analysis, we 
obtain 

a;b(CMB) = 0.022 ± 0.002 ± 0.001 , (1.70) 
values which are in excellent agreement. 

Matter density. Combining the results for of the three cited joint compilations 
of CMB and LSS power spectra, we find 

= 0.29 ± 0.05 ± 0.04 . (1.71) 

This is in agreement with the constraint on the total density of clustered matter in the 
Universe which comes from the combination of x-ray measurements of galaxy clusters 
with large hydrodynamic simulations.^^ 

^^Using measurements of both the temperature and luminosity of the x-rays coming from hot gas 
which dominates the baryon fraction in clusters, under the assumption of hydrostatic equilibrium, 
one obtains the gravitational potential, and in particular the ratio of baryon to total mass of these 
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Vacuum density. Whatever it is, cosmological constant or quintessence, the vacuum 
energy is today the dominant form of energy density of the Universe. In fact, given 
Q,\ — fto — Qm, combining the results (1.68) and (1-71), we obtain 

= 0.71 ±0.07 . (1.72) 

Furthermore, a non-zero cosmological constant is favoured by the type-la supernovae 
measurements (Perlmutter et al, (1998) [142]; Schmidt et al (1998) [157]). 

Spectral index. In the usual assumption of scalar adiabatic perturbations (see next 
chapter) , which provide good fits with experimental data and are exactly what inflation 
predicts, tensorial spectral index is obviously rif — 0, while for the scalar spectral index 
we obtain, from the cited joint compilations. 

Us = 1.02 ± 0.06 ± 0.05 , (1.73) 

again in agreement with the inflationary theory, which predicts a scale invariant spec- 
trum. 

Hubble parameter. The situation for the determination of Hq is slightly different 
from the previous ones, as in this case the concordance of different methods is less 
evident. A great effort in this direction was made by the Hubble Space Telescope (HST) 
Key Project [77], which, using a calibration based on revised Cepheid distances, found 
the value of Hq for five different methods (type la supernovae, Tully-Fisher relation, 
surface brightness fluctuations, type II supernovae, and fundamental plane), with a 
weighed average given by 

/^o = 72 ± 8 km s-^Mpc"^ . (1.74) 

However, significantly smaller values of Hq are obtained by other methods, such 
as measurements of time delays in gravitationally lensed systems, or the Sunyaev- 
Zeldovich (SZE) effect in X-ray emitting galaxy clusters, which bypass the traditional 
"distance ladder" and probe far deeper distances than the objects used by the Key 
Project. Recently, Kochaneck used five well-constrained gravitational lenses, obtaining 

//o = 62 ± 7 km s'^Mpc"^ . (1.75) 

Measurements of the SZE in 14 clusters also indicate a value of -f^o ~ 60 km s^^ Mpc~^ 

with presently a large (~ 30%) systematic uncertainty [39]. Rowan- Robinson has 

systems. Employing the constraint on the cosmic baryon density coming from BBN, and assuming 
that galaxy clusters provide a good estimate of the total clustered mass in the Universe, one can then 
arrive at a range for the total mass density in the Universe [62, 63, 107, 183]. 
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argued that the Key Project data need to be corrected for local peculiar motions using 
a more sophisticated flow model than was actually used, and also for metallicity effects 
on the Cepheid calibration, then lowering the value of Hq inferred from the same 
dataset to 63 ± 6 km s~^ Mpc~^ [151]- Furthermore, it is worth noting that all joint 
analysys of CMB and LSS indicate a value of Hq closer to the lower estimation. 



The figure below shows an histogram which compares the amounts of the different 
components of the Universe. The second component from left is the dark matter 
(this argument will be treated in the next chapter), which together with the third 
component, the baryonic matter, constitutes the total matter content of the Universe. 
We can also note that the vacuum energy density is today comparable to (or better 
higher than) matter density. 

At this point, the raw conclusion is very simple: the Universe is spatially flat, some 
70% of the total energy content is dark, possibly in the form of a cosmological constant, 
and some 25% of gravitating matter is dark and unknown! 



Energy Densities 



Figure 1.1: Histogram of the composition of the Universe, given, from left to right, by: 
energy of the vacuum (70% ± 10%), dark matter (26% ± 10%), ordinary baryonic matter 
(4% lb 1%), and negligible energy density from photons. About 13% of the matter in the 
Universe is baryonic (J7b/^m = 0.04/0.3 ~ 0.13). The baryons can be further divided into 
3% warm invisible gas, 0.5% optically visible stars and 0.5% hot gas visible in the x-rays [80]. 
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Cosmological structure formation 
and dark matter 



2.1 Dark matter 

With the term dark matter cosmologists indicate an hypothetic material component of 
the Universe which does not emit directly electromagnetic radiation (unless it decays 
in particles having this property). 

Historically, the first evidences of the existence of some form of dark matter comes 
from measurements of velocity dispersion of galaxies in clusters (firstly obtained by 
Zwicky in 1933 [189]) and the surprisingly fiat rotation curves of spiral galaxies^. What 
we deduce is that on galactic scales and above the mass density associated with lumi- 
nous matter (stars, hydrogen, clouds, x-ray gas in clusters, etc.) cannot account with 
the observed dynamics on those scales. 

Prom figure 1.1 and from the parameter estimation of the §1.7 we clearly see the 
discrepancy existent between the total matter density and the baryon density. The 
conclusion is that most of the matter in the Universe is dominated by an unknown 
form of matter or by an unfamiliar class of dark astrophysical objects [106]. 

There are essentially two ways in which matter in the Universe can be revealed: by 
means of radiation, by itself emitted, or by means of its gravitational interaction with 
baryonic matter which gives rise to cosmic structures. In the first case, electromagnetic 

"'^ Rotation curves of galaxies arc characterized by a peak reached at distances of some Kpc and a 
behaviour typically flat for the regions at distance larger than that of the peak. A peculiarity is that 
the expected Keplerian fall is not observed. These rotation curves are consistent with extended halos 
containing masses till 10 times the galactic mass observed in the optical [175]. 
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radiation permits to reveal only baryonic matter. In the second case, we can only tell 
that we are in presence of matter that interacts by means of gravitation with the 
luminous mass in the Universe. 

The most widespread hypothesis is that dark matter is in form of coUisionless parti- 
cles, but we still do not exactly know neither the nature nor the masses of them. High 
energy physics theories, however, provide us with a whole "zoo" of candidates which 
are known as cosmic relics or relic WIMPs (Weakly Interacting Massive Particles). 
Typically, we distinguish between "thermal" and "non-thermal" relics. The former are 
kept in thermal equilibrium with the rest of the Universe until they decouple (a char- 
acteristic example is the massless neutrino), while the latter (such as axions, magnetic 
monopoles and cosmic strings) have been out of equilibrium throughout their lifetime. 
Thermal relics are further subdivided into three families. 

• Hot Dark Matter (HDM): "hot" rehcs which are still relativistic when they de- 
couple. A typical hot thermal rehc is a hght neutrino with rrii, ~ 10 eV. 

• Cold Dark Matter (CDM): "cold" relics which go non-relativistic before 
decoupling.^ The best motivated cold relic is the lightest supersymmetric partner 
of the standard model particles^, identified in the "neutrahno" , with m > 1 GeV. 

• Warm Dark Matter (WDM): the intermediate case, given by thermal relics with 
masses around 1 keV. Right-handed neutrinos, axinos and gravitinos have all 
been suggested as potential warm relic candidates. 

The study of dark matter has as its finality the explanation of formation of galaxies 
and in general of cosmic structures. For this reason, in the last decades, the origin of 
cosmic structures has been "framed" in models in which dark matter constitutes the 
skeleton of cosmic structures and supply the most part of the mass of which the same 
is made. In this chapter, after a brief review of linear perturbation theory (for more 
details see appendix B), we will see some features of the most commonly used dark 
matter scenarios. 

^In general non-thermal relics are also indicated as CDM. 

■^Supersymmetric particles were postulated in order to solve the strong CP problem in nuclear 
physics. This problem arises from the fact that some interactions violate the parity P, time inversion 
T, and CP. If these are not eliminated, they give rise to a dipole momentum for the neutron which is 
in excess often order of magnitude with respect to experimental limits (Kolb and Turner, 1990 [104]). 
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2.2 Introduction to structure formation 

We know that the Universe was very smooth at early times (as recorded by the CMB) 
and it is very lumpy now (as observed locally). Cosmologists believe that the reason is 
gravitational instability: small fluctuations in the density of the primeval cosmic fluid 
grew gravitationally into the galaxies, the clusters and the voids we observe today. 
However, there is the lack of the exact scenario, which is to say the lack of the exact 
composition of the Universe ("dark matter" and "dark energy" problems). 

The idea of gravitational instability was flrst introduced in the early 1900s by Jeans, 
who showed that a homogeneous and isotropic fluid is unstable to small perturbations 
in its density [96, 97]. What Jeans demonstrated was that density inhomogeneities 
grow in time when the pressure support is weak compared to the gravitational pull. 

The Newtonian theory, the extension of Jeans theory to an expanding Universe, is 
only applicable to scales well within the Hubble radius, where the relativistic effects 
are negligible, and, even in this context, one can only analyze density perturbations in 
the non-relativistic component. Perturbations on scales over the Hubble radius, or in 
the relativistic matter at all scales, require the full relativistic theory. 

Using these theories, if we want to construct a detailed scenario of structure for- 
mation, we need to know: (i) the composition of the Universe; (ii) the contribution of 
its various components to the total density (namely Jib from ordinary baryons, f^wiMP 
from relic WIMPs, Q-y from relativistic particles, Qa from a potential cosmological con- 
stant, etc.); (iii) the spectrum and the type (adiabatic or isocurvature) of the primeval 
density perturbations. 

Adiabatic and isocurvature perturbations. In general, when deahng with the 
pre-recombination plasma, we distinguish between two types of perturbations, namely 
between "isoentropic" (adiabatic) and "entropic" (isocurvature or isothermal) modes 
(Zeldovich, 1967) [187]. Before recombination, a generic perturbation can be decom- 
posed into a superposition of independently propagating adiabatic and entropic modes, 
but, after matter and radiation have decoupled, perturbations evolve in the same way 
regardless of their original nature. 

"Adiabatic" modes contain fluctuations both in the matter and the radiation com- 
ponents, while keeping the entropy per baryon conserved. So, if we remember eq. 
(A. 48), deflning, as usual, S = (Sp/p), we find the condition for adiabaticity 




(2.1) 
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These perturbations are naturally generated in the simplest inflationary models 
through the vacuum fluctuation of the inflaton field (see Liddle & Lyth (2000) [113]). 
In the "isocurvature" modes we have 

Sp^O^ p^5^ + p^5y^0 ^ -1 = -^. (2.2) 

Oh 

This implies that the geometry of the 3-dimensional spatial hypersurfaces remains 
unaffected, hence the name isocurvature. Unlike adiabatic disturbances, isocurvature 
perturbations are usually absent from the simplest models of inflation. 

In this thesis we study only adiabatic perturbations, which are today the preferred 
perturbation modes, and we leave out the isocurvature modes, which could also have 
a contribution, but certainly cannot be the dominant component [60, 61]. 

2.3 Linear evolution of perturbations 

As long as the inhomogeneities are small, they can be studied by the linear perturbation 
theory. A great advantage of the linear regime is that the different perturbative modes 
evolve independently and therefore can be treated separately. In this respect, it is 
natural to divide the analysis of cosmological perturbations into two regimes. The 
early phase, when the perturbation is still outside the horizon, and the late time 
regime, when the mode is inside the Hubble radius. In the first case microphysical 
processes, such as pressure effects for example, are negligible and the evolution of the 
perturbation is basically kinematic. After the mode has entered the horizon, however, 
one can no longer disregard microphysics and damping effects. 

The Newtonian treatment suffices on sub-horizon scales and as long as we deal with 
fluctuations in the non-relativistic component, while general relativity is necessary on 
scales outside the horizon and also when studying perturbations in the relativistic 
component. Here we remember only the principal results of these theories and some 
important concept, as the Jeans length, referring to appendix B for a more quantitative 
discussion. 

The complete set of the relativistic equations reveals three types of perturbations: 
tensor, vector and scalar modes. Tensor perturbations correspond to the traceless, 
transverse part of 5gaf3- They describe gravitational waves and have no Newtonian 
analogue. Vector and scalar perturbations, on the other hand, have Newtonian coun- 
terparts. Vector modes correspond to rotational perturbations of the velocity field. 
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Table 2.1: Growth of perturbations. 



Epoch 






5b 


a < Oenter < ^eq (A > Ah) 


grows as 


grows as 


grows as 


Oenter < a < tt^q (A < Ah) 


oscillates 


grows as In a 


oscillates 


Oeq < a < Odec (A < Ah) 


oscillates 


grows as a 


oscillates 


Odec < a (A < Ah) 


oscillates 


grows as a 


grows as a 



while scalar modes are associated with longitudinal density fluctuations. We will only 
consider the latter type of perturbations. 

As explained in appendix B, the key scales for the study of the perturbation evo- 
lution are the Jeans length Aj and the related Jeans mass Mj 

Aj = t;s,/^ Mj = ^7rp(:^| , (2.3) 




, Gpo 3 \ 2 

where fg is the adiabatic sound speed defined in eq. (B.ll), po is the unperturbed 
density and p is the density of the perturbed component. The Jeans scale constitutes a 
characteristic feature of the perturbation. It separates the gravitationally stable modes 
from the unstable ones: fluctuations on scales well beyond Aj grow via gravitational 
instability, while modes with A ^ Aj are stabilized by pressure. 




scale factor 



Figure 2.1: Evolution of the perturbation amplitudes on a scale M ~ 10^^ for the dark 
matter component (green line), baryons (red) and radiation (black). 



Starting from the results obtained in appendix B, we summarize the trends of 
density contrasts for radiation, dark matter and baryons in table 2.1, where Oenter 
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indicates the epoch at which a given perturbation enters the Hubble radius Ah- In 
figure 2.1 we depict the perturbation evolution of the same components. 

Before the horizon entry fiuctuations in all components grows as a^. After this 
epoch baryons and photons enter in causal contact and, being tightly coupled, begin 
to oscillate until the moment of decouphng. Meanwhile, given that WIMPs interact 
with matter by means of gravity only, and do not feel the effect of pressure forces 
due to interaction with radiation, dark matter perturbations continue to grow, but 
whit a slower rate, as In a between horizon entry and matter-radiation equivalence and 
as a after equivalence. Finally, after baryon to photon decoupling, photons continue 
to oscillate, while perturbations in the baryonic distribution will be driven by the 
gravitational potential of the coUisionless species, and grow rapidly to soon equalize 
with those in the dark matter {dark matter boost). Subsequently, perturbations in 
both components grow proportionally to the scale factor. Note that the presence of 
dark matter is necessary, because without it baryonic fiuctuations do not have enough 
time to produce the observed structures compatibly with the extreme smoothness of 
the CMB temperature (see next chapter) . In this context the key role played by dark 
matter is to start growing perturbations in its density earlier than those in the baryonic 
component; in chapter 5 we will see that the same result could also be obtained with 
a particular self-coUisional form of dark matter. 

2.4 Non-baryonic structure formation 
2.4.1 Evolution of the Jeans mass 

If the dark component is made up of weakly interacting species, the particles do not feel 
each other's presence via collisions. Each particle moves along a spacetime geodesic, 
while perturbations modify these geodesic orbits. We can treat the coUisionless species 
as an ideal fiuid, so the associated Jeans length is obtained similarly to the baryonic 
one. When dealing with a coUisionless species, however, one needs to replace eqs. (B.l), 
(B.2) with the Liouville equation (see Coles & Lucchin (1995) [47]). Then, 



where now Vdm is the velocity dispersion of the dark matter component. The corre- 
sponding Jeans mass is 




(2.4) 




(2.5) 
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where pdm is the density of the non-baryonic matter. 

Hot thermal relics. They decouple while they are still relativistic, that is, defining 
tor as the time when the particles become non-relativistic, t^ec < tm, where we assume 
that tnv < ieq- Throughout the relativistic regime v^m. ~ 1, p — P7 oc and 
Pdm oc a~^, implying that Aj oc and mS"^ (X a . Once the species have become 
non-relativistic and until matter-radiation equality, Vdm oc (due to the redshifting 
of the momentum p oc a"^ when the particles are non-relativistic), p ~ oc 
and pdm oc a~^. Recall that the particles have already decoupled, which means that 
^dm 7^ T^- Consequently, Aj'^ oc a and Mj^^ — const. After equipartition v^m cc a"^ 
and p ~ Pdm oc a'^, which translates into Aj'^ oc a^/^ and Mj^^ oc a~^/^. Overall, the 
Jeans mass of hot thermal relics evolves as 



m]^) oc <; 



constant Onr < « < Oeq ; (2.6) 

a~^/^ Oeq < a . 



Clearly, m]^^ reaches its maximum at Onr- In fact, the highest possible value corre- 
sponds to particles with = such as neutrinos with rrii, ~ 10 eV. In this case 
(Mj''^)max ^ 3.5 X 10i5(n^/i2-)-2 (ggg Q^^gg ^ Lucchiu (1995) [47]). For a typical 
hot thermal relic (Mf ^)max ~ 10^^ - 10^^ Mq. 

Cold thermal relics. Cold thermal relics decouple when they are already non- 
relativistic (i.e. tnr < tdec < ^eq)- ThuS, for t < tnv WC haVC V^m ~ 1, Pdm OC 

and p ~ p^ oc a~^, implying that Aj^^ oc and Mj^^ (x a?. In the interval between 
inr and tdec the key variables evolve as I'dm oc a'^^'^ (recall that Tdm — until idee)) 
Pdm oc and p ~ p^ oc a"^. As a result, Aj'^^ oc a^/^ and Mf^ oc a^/^. After the 
particles have decoupled Tdm 7^ T^, which means that Vdm oc a~^. At the same time 
Pdm oc and p o: a~^, ensuring that Aj'^^ oc a and Mj^^ = constant. After 

equality I'dm oc a"^ and p oc pdm oc a~^, implying that Aj^^ oc a^/^ and Mj^^ oc a~^/^. In 
short, the Jeans mass of cold thermal relics evolves as 



m]"^ oc i 



Oi ^ Ztox > -Z > ^dec ) ('2 7) 

constant 2;dec > > -^^eq ; 



a > ^ 



Accordingly, the maximum value for M^^ corresponds to species with tdec = ^eq- In 
other words, the sooner the particles cease being relativistic and decouple, the smaller 
the associated maximum Jeans mass. Typically {M'f')^^^ <^ 10^^ M0. 
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2.4.2 Dissipative effects 

The ideal fluid approximation for coUisionless species holds on sufficiently large scales 
only. On small scales, the free geodesic motion of the particles will wipe out any struc- 
ture. This phenomenon, known as "Landau damping" or "free streaming" , consists in 
the smoothing of inhomogeneities in the primordial Universe because of the motion of 
coUisionless particles from overdense to underdense regions. 

Consider the coordinate (comoving) distance traveled by a free streaming particle 

XFS = /' —dt , (2.8) 
Jo a 

where Aps = dxps is the corresponding physical (i.e. proper) distance. Clearly, pertur- 
bations in the dark matter component on scales smaller than Aps will be wiped out by 
free streaming. Integrating the above during the three intervals t < tnv , t^v < t < teq 
and t > teq , we find the total physical free-streaming scale 



Aps = 



2inr 



1 + lnf^ 

Vanr 



1 - 



a . (2.9) 



At late times, when a ^ a^q, the above approaches its maximum value 



Aps (Aps)inax - -r 

Cl-n 



5 + 21n(^ 



(2.10) 



To obtain numerical estimates we need to identify the epoch the species become non- 
relativistic. Assuming that the transition takes place when T ~ mdm/3, we find (see 
Padmanabhan (1993) [132]) 

(Afs)_ - 0.5 (^)"''' {^...hf Mpc , (2.11) 

where is the mass of the coUisionless particles. Accordingly, the minimum scale 
that survives coUisionless dissipation depends crucially on the mass of the dark matter 
species. For neutrinos with ~ 30 eV we find (AFs)max — 28 Mpc and a corresponding 
mass scale (Mps)max ~ 10"*^^ Mq. For a much heavier candidate, say ~ 1 keV, 
we find (Aps)max ~ 0.5 Mpc and (Mps)max ~ 10^ Mq. In general, the hghter the dark 
matter species less power survives on small scales. 

Cold thermal relics and non-thermal relics have very small dispersion velocities. As 
a result, the maximum values of the Jeans mass and of the free streaming mass are very 
low. In this case, perturbations on all scales of cosmological interest grow unimpeded 
by damping processes, although they suffer stagnation due to the Meszaros effect until 
matter-radiation equality. After recombination the potential wells of the coUisionless 
species can boost the growth of perturbations in the baryonic component on scales of 
the order of (Mj)rec ~ 10^ Mq. 
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2.4.3 Evolution of the Hubble mass 

We can build a useful quantity analogous to the Jeans mass by means of the Hubble 
radius Ah, and call it the Hubble mass 

MH = ^7rp(^y, (2.12) 

where p is the energy density of the perturbed species. This is an additional important 
scale for structure formation; Ah and Mh define the scale over which the different parts 
of a perturbation are in causal contact. Note that a mass scale M is said to be entering 
the Hubble radius when M — Mh- 

For t < tnv we have p oc a~^, Ah oc i oc a^. During the interval t^r < t < teq we 
have p (X and Ah oc i oc a^. Finally, when teq < t, p oc and Ah oc i oc a^/^. 
Overall, taking in mind that Onr and Oeq depend on the species, the Hubble mass of 
every component evolves as 



Mh oc < 



a < a 

,3 



nr , 



a"* Onr < a < Oeq , (2-13) 



7.3/2 



Oeq < a . 



Following definitions (2.4), (2.5) and (2.12), one can easily verify that the Jeans mass 
and the Hubble mass are efi:ectively identical as long as the relic species are relativistic, 
namely a < Onr- 

2.4.4 Scenarios, successes and shortcomings 

The study of origin and formation of structures in the Universe has been historically 
fundamentally framed into the two HDM and CDM scenarios, according to what is the 
dominant form of dark matter. Although both theories have the same starting points 
(flat Universe fundamentally constituted by dark matter; small baryonic contribution 
to the mass of the Universe; primordial fluctuations adiabatic, scale invariant and 
Gaussian), structure formation in these scenarios is completely different since hot and 
cold relics are subject to different physical phenomena. 

As shown in figure 2.2, there are essentially three stages in the evolution of a mode 
which enters the Hubble radius between a^r and Oeq: (A) {t < tent < ^eq) the wavelength 
of the perturbation is bigger than the Hubble radius, and the mode grows as a^; (B) 
{tent < t < teq) the wavclcngth is inside the Hubble radius and is bigger than Aj, so 
pressure support cannot stop the collapse, but the perturbation is frozen in due to the 
Meszaros effect; (C) {teq < t) the wavelength is inside the Hubble radius and is bigger 
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than A J, the mode becomes unstable again and grows as a. Note that fluctuations with 
size smaller than Mps are wiped out by neutrino free streaming. 




Figure 2.2: Evolution of the Hubble and Jeans mass for dark matter, with indicated a 
perturbed scale M = M(A), the three stages of the text and the neutrino free streaming 
mass. 



Hot dark matter models. Keeping a light neutrino species with my ~ 10 eV (which 
implies ^ acq), the key feature of the perturbation spectrum (see chapter 3) is 
the cutoff at Afs — 40 Mpc due to the neutrino free streaming. Because of this, 
the first structures to form have sizes of approximately 10^^ M0, which corresponds 
to a supercluster of galaxies (see figure 2.2). Moreover, because the scale is very 
large, collapse must have occurred at relatively recent times (i.e. at z < 3). Thus, in 
a Universe dominated by hot thermal relics, structure formation proceeds in a "top- 
down" fashion. Perturbations on scales as large as 10^^ Mq go no linear in a highly non- 
spherical way. As a result, they collapse to one dimensional objects called "pancakes" 
(Zeldovich (1970) [188]). Once the pancake forms and goes non-linear in one of its 
dimensions, the baryons within start colliding with each other and dissipate their 
energy. Thereby, the baryonic component fragments and condenses into smaller galaxy- 
sized objects. 

These models have problems reproducing the small scale clustering properties of 
galaxies. In particular, the HDM simulations can agree with the observed galaxy- 
galaxy correlation function (see chapter 3) only if the epoch of pancaking takes place 
at z ~ 1 or less. This seems too late to account for the existence of galaxies with 
redshift z > 1 and of quasars with 2; ~ 5. 

Cold dark matter models. For such particles the maximum damping scale is too 
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small (<S 1 Mpc) to be of any cosmological relevance, so that the CDM model has 
a spectrum without a cut-off at short wave-length (at least till scales much smaller 
than the galactic ones) . Another important feature is the weak growth experienced by 
perturbations between horizon crossing and equipartition. It means that the density 
contrast increases as we move to smaller scales, or that the perturbation spectrum has 
more small-scale power. Thus, in standard CDM scenarios the first objects to break 
away from the background expansion have sub-galactic sizes (< 10^ Mq). These struc- 
tures virialize through violent relaxation (Lynden-Bell (1967) [117]; Shu (1978) [161]) 
into gravitationally bound configurations that resemble galactic halos. At the same 
time the baryons can dissipate their energy and condense further into the cores of these 
objects. As larger and larger scales go non-linear, bigger structures form through tidal 
interactions and mergers. Hence, according to the CDM scenario structures form in a 
"bottom-up" or "hierarchical" fashion. 

Dcirk matter model problems. Today the standard model consists of an Universe 
in which all the dark matter is made of CDM, and the remaining energy density 
needed to reach the critical density is provided by the dark energy: this is the so- 
called ACDM scenario. This model obtained noteworthy successes in the description of 
the characteristics of the Universe (clustering statistics of galaxies, peculiar velocities, 
CMB fluctuations) from the galactic scale on [12, 30, 56, 79, 137, 184], but it has 
shown various discrepancies, when compared with data on smaller (subgalactic) scales, 
at least when any bias of the distribution of galaxies relative to the mass (see chapter 3) 
is constant with scale [7, 35, 52, 53]. Some of these problems come from the following. 

1. A^-body CDM simulations give cuspy halos with divergent profiles toward the 
galactic center [130], in disagreement with the galaxy rotation curves [25] and 
with observations from gravitational lensing [66]. In fact, in CDM halos we find 
p(r) ~ with a ~ 1, so that the rotation velocity at the centers of galaxies 
should increase as r^/^, but the data, especially that on dark-matter-dominated 
dwarf galaxies, instead showed a linear increase with radius, corresponding to 
roughly constant density in the centers of galaxies [126]. 

2. One of the most striking features of halos in high resolution CDM simula- 
tions [100] is that they are heavily populated with small subhalos or subclumps, 
but their number is 1-2 orders of magnitude more than what is observed [127], 
leading to an evident excess of small scale structures [156]. Even if there are 
mechanisms which inhibit star formation in small clumps and thereby leave them 
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dark, a large population of clumps may heat the galactic disk and even endanger 
its stabihty [129, 174, 182]. 

3. The distribution of specific angular momentum in dark matter halos has a univer- 
sal profile [36] . But if the baryons have the same angular momentum distribution 
as the dark matter, this implies that there is too much baryonic material with 
low angular momentum to form the observed rotationally supported exponential 
disks [36, 176]. 

4. Bar stability in high surface brightness spiral galaxies also demands low-density 
cores [50]. 

It is conceivable that these discrepancies are due to problems with current simu- 
lations, the quality of present observations or the omission of important astrophysical 
processes in the models. However, if any one of them persists, they may be an indica- 
tion that the dark matter is not collisionless. 

Moreover, in the last years studies on predicted 7-ray emission from CDM annihi- 
lation in clumps are severely constraining the CDM particles parameters in the light 
of recent (EGRET) and future (GLAST) experiments, giving important indication on 
the reliability of present CDM candidates [4, 171]. 

Alternative options. Issues of the ACDM model that have arisen on small scales 
have prompted people to propose alternatives to CDM, such as "warm dark matter" 
(WDM) [31] and "self-interacting dark matter" (SIDM) [169]. 

Warm dark matter particles have relatively high thermal velocities, so free stream- 
ing can suppress the formation of structure on small scales, while on larger scales the 
spectrum is the same as for the CDM case. Simulations [48, 31] show that there are far 
fewer small satellite halos with AWDM than ACDM. Dark matter halos nevertheless 
have density profiles much like those in CDM [93, 128, 58]. Hydrodynamical simula- 
tions also indicate that the disk angular momentum problem may be resolved with this 
suppression of small scale power [168]. 

While CDM assumes that the dark matter particles have only weak interactions 
with each other and with other particles, SIDM assumes that they have strong elastic 
scattering cross sections, but negligible annihilation or dissipation. In this way SIDM 
might suppress the formation of the dense central regions of dark matter halos. 

Finally, another dark matter candidate is the so-called mirror dark matter, which 
in principle could contain the virtues of all the previous CDM, WDM and SIDM. His 
study is just the aim of this thesis! 
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CMB, LSS and evolutionary 
equations 

3.1 The cosmic microwave background (CMB) 

As explained in §1.5, at recombination the photons decoupled from the baryons, and 
began to travel nearly unperturbed from the last scattering surface^ until the present 
day, where we observe them peaked in the microwave region of the spectrum as the 
cosmic microwave background radiation (CMB). Thus, they provide a fossil record of 
our observable Universe when it was roughly 10^ years old and about 10^ times smaller 
than today. 

The existence of this background of cold photons (with a present number density 
n-y ~ 420 cm~^), discovered in 1965 by Penzias and Wilson [139], was predicted several 
years before by Gamow, Alpher and Herman [5, 83] as a consequence of the hot Big 
Bang theory. 

The cosmic background radiation is the most perfect blackbody ever seen, according 
to the FIRAS (Far InfraRed Absolute Spectrometer) instrument [65] of the Cosmic 
Background Explorer (COBE) sateUite, which measured in the early 90's a temperature 
of To = 2.725 ± 0.002 K [120]. 

In addition to the FIRAS results, the other major result of COBE comes from 
the DMR (Differential Microwave Radiometer) instrument, which measured a dipole 
(ATdip = 3.353±0.024 mK) associated with the 300 km/s flow of the Earth in the CMB 
frame, an intrinsic quadrupole {£ — 2) amphtude AT2 = lO.OlU A*K, and provided the 
first unambiguous detection of anisotropies at a level AT/T 10~^ on large angular 
scales (> 7°) [167]. This breakthrough immediately stimulated the realization of many 

^The sphere surrounding us at 2: ~ 1100, which represents the position at which the CMB photons 
seen today last interacted directly with matter, is called the "last scattering surface". 
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new experiments aiming at measuring the CMB angular distribution with increasing 
resolution and sensitivity, revealing temperature anisotropics on smaller angular scales, 
which correspond to the physical scale of observed structures such as galaxies and 
clusters of galaxies. 

3.1.1 Anisotropies and power spectrum 

Theories of the formation of large-scale structure predict the existence of shght in- 
homogeneities in the distribution of matter in the early Universe, which eventually 
underwent gravitational collapse to form galaxies, galaxy clusters, and superclusters. 
These density inhomogeneities lead to temperature anisotropies in the CMB, because 
the radiation leaving a dense area of the last scattering surface is gravitationally red- 
shifted to a lower apparent temperature and viceversa for an under dense region. 

Mapping the temperature on the celestial sphere, the temperature anisotropy at a 
point on the sky {9, (f)) can be expressed in the basis of spherical harmonics as 

-7fr{0A) = Y. E ai^Yi^iOA), (3.1) 

^ 1=0 m=-l 

where i ~ 180°/^, and aem represent the multipole moments that, in the case of a 
Gaussian distribution of density perturbations (as suggested by inflationary models), 
should be characterized [47] by zero mean, {a^rn) — 0, and non-zero variance 

Ce = {\aem\') (3.2) 

(the angle brackets indicate an average over all observers in the Universe; the absence 
of a preferred direction implies that (|a^mP) should be independent of m). The set 
of Ce is known as the angular power spectrum. A cosmological model predicts the 
amplitude of the aim coefficients, which can be expressed in terms of Ce alone in the 
hypothesis of gaussianity. 

Temperature fluctuations in the CMB around a mean temperature in a direction a 
on the sky can be analyzed in terms of the correlation function C{9) 

C(^) = (f^(«)f^(« + ^)) , (3.3) 

which measures the average product of temperatures in two directions separated by an 
angle 9. For small angles 9 the temperature correlation function can be expressed as a 
sum of Legendre polynomials Pi{9) of order £, with coefficients or powers a^, 

-| oo 

C{9) = —Y.a]{2^+l)Pt{cos9) . (3.4) 
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All analyses start with the quadrupole mode £ — 2 because the £ — monopole mode 
is just the mean temperature over the observed part of the sky, and the £ — 1 mode is 
the dipole anisotropy due to the motion of Earth relative to the CMB. The higher the 
angular resolution, the more terms of high £ must be included. 

Anisotropics on scales 0.1° ^ ^ ^ 2° are related to causal processes occurring in 
the photon-baryon fluid until recombination. Photons and baryons are in fact tightly 
coupled and behave like a single fluid. In the presence of gravitational potential forced 
acoustic oscillations in the fluid arise: they can be described by an harmonic oscilla- 
tor where the driving forces are due to gravity, inertia of baryons and pressure from 
photons. Recombination is a nearly instantaneous process and modes of acoustic os- 
cillations with different wavelength are "frozen" at different phases of oscillation. The 
first (so-called Doppler) peak at degree scale in the power spectrum is therefore due to 
a wave that has a density maximum just at the time of last scattering; the secondary 
peaks at higher multipoles are higher harmonics of the principal oscillations and have 
oscillated more than once. A precise measurement of the acoustic peaks can reveal 
information on the cosmological parameters. 

In the short but finite time taken for the Universe to recombine, the photons can 
diffuse a certain distance. Anisotropics on scales smaller than this mean free path 
will be erased by diffusion, leading to the quasi-exponential damping [165, 92] in the 
spectrum at large ts. This is called "Silk damping" (see chapter 5) and becomes quite 
effective ai £> 1000, corresponding to angular scales 9 < 10'. 

In all models the power spectrum contains the bulk of the statistical content in 
CMB maps. A statistically isotropic sky means that all the {2£ + 1) m's are equivalent 
and the power at each £ can be written as {2£ + l)Cf/(47r). For an idealized full-sky 
observation, the variance of each measured Q is [2/{2£+l)]C^. This samphng variance 
(known as cosmic variance), which comes about because each is distributed with 
{2£ + 1) degrees of freedom for our observable volume of the Universe [185], sets the 
ultimate limit on the accuracy of our estimates of the power spectrum.^ 

3.2 The large scale structure (LSS) 

Fourier description and power spectrum 

Let us now turn to the distribution of matter in the Universe, where naturally appear 
random fiuctuations around the mean density p{t) , manifested by compressions in some 

^Of course similar considerations to those discussed above apply when an only patch of the sky is 
available, so temperature maps may have sampling variance too. 
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regions and rarefactions in other regions. 

It is generally assumed that this distribution is given by the superposition of plane 
waves, each with its characteristic wavelength A or comoving wave number k and its 
amplitude 5k, that evolve independently, at least until they are in the linear regime 
(see §2.3). Let we divide the Universe in cells of volume Vu (for example a cube of side 
L) and impose periodic conditions on the surfaces. We define, as usual, the density 
contrast as 

^(x) = . (3.5) 

p 

and we assume this to be expressible as a Fourier expansion 

5(x) = ^ exp(ik • x) = ^ (5k exp(— ik • x) , (3.6) 

k k 

where ^1,^,2 = L and, for the periodicity condition, 5(x, y, L) = S{x, y, 0) (and 

similar conditions for the other components). The inverse relationship of the (3.6) 
gives the Fourier coefficients 6^ 

(5k = — / 5(x) exp(— ik ■ x)(ix , (3.7) 
Vu Jvu 

which are complex quantities 

4 = l^kl exp (i0k) • (3.8) 

The assumption of periodic boundaries results in a discrete k-space representation. 
Conservation of mass in Vu implies (5k=o = and the reality of 5{x) requires (5k = (5_k- 
If we consider n volumes Vu, we have the problem of determining the distribution 
of Fourier coefficients (5k and that of \5\. If we suppose that phases 0k are random, in 
the hmit T^u — ^ 00 it is possible to show that we get \Sf — Z^k l<^k|^- The Central Limit 
theorem leads us to conclude that the distribution for 5 is Gaussian 

m « exp (^-^j . (3.9) 
The quantity a present in eq. (3.9) is the variance of the density field, defined as 

^'-{s')-T.m') = ^T.sl. (3.10) 

k k 

This quantity characterizes the amplitude of the inhomogeneity of the density field. If 
Vu — > 00, we obtain the more usual relation 

(7^ = ^ I p{k)d'k = ^ I p{k)edk . (3.11) 
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The term P(k) — {\5kf) is called power spectrum of perturbations. It is function only 
of k because the ensemble average in an isotropic Universe depends only on x. 

With the lack of more accurate knowledge, one assumes for simplicity that the 
power spectrum of primordial gravitational fluctuations is specified by a power law 

P{k)=Ak''% (3.12) 

where Us is the spectral index of scalar fluctuations and A is the amplitude, which is 
expected to be equal on all scales. Inflationary models also predict that the power 
spectrum of matter fluctuations is almost scale-invariant as the fluctuations cross the 
Hubble radius. This is the Harrison-Zeldovich spectrum, for which Ug = I. 



Correlation function 

An important quantity connected with the spectrum is the two-points correlation func- 
tion ^(r), which describes the real-space statistical properties of spatial density per- 
turbations. It can be defined as the joint probability of finding an overdensity 6 in two 
distinct points of space 

where averages are averages on an ensemble obtained from several realizations of Uni- 
verse. Correlation function can be expressed as the joint probability of finding a galaxy 
in a volume 6Vi and another in a volume 6V2 separated by a distance ru 

5^P = n^[l + ^{ru)]SVi5V2 , (3.14) 

where ny is the average number of galaxies per unit volume. The concept of correlation 
function, given in this terms, can be enlarged to the case of three or more points. 

Correlation functions have a fundamental role in the study of clustering of matter. 
In order to show the relation between perturbation spectrum and two-points correlation 
function, we apply the Fourier machinery to equation (3.13) and arrive at the relation 

e(r)=E(l<^kr)exp(-ik.r), (3.15) 

k 

which, in passing to the limit Vu —>■ 00, becomes 

C(r) = J P{k) exp(-^k ■ r)cik . (3.16) 
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This result shows that the two-point correlation function is the Fourier transform of 
the power spectrum (this is the so-called Wiener-Khintchin theorem) . This is similar 
to the situation in the context of CMB anisotropies, where the waves represented 
temperature fluctuations on the surface of the surrounding sky, and the powers 
were coefficients in the Legendre polynomial expansion. 

In an isotropic Universe, it is |r| = r and then |k| = k, and the spectrum can be 
obtained from an integral on |k| = k, so that correlation function may be written as 



Phenomenological models of density fluctuations can be specifled by the amplitudes 
Sk of the correlation function ^(r). In particular, if the fluctuations are Gaussian, they 
are completely specifled by the power spectrum P{k). The models can then be com- 
pared to the real distribution of galaxies and galaxy clusters, and the phenomenological 
parameters determined. 

Transfer function 

During the evolution of the Universe and after perturbations enter the horizon, the 
spectrum is subject to modulations because of physical processes characteristic of the 
model itself (free streaming for acoUisional components. Silk damping for coUisional 
particles, etc.; see previous and next chapters). Ampliflcation, stagnation or damping 
cause an evolution of the fluctuations different from one spatial scale to another. The 
combined effect of the various processes involved in changing the shape of the original 
power spectrum can be summarized in a single quantity, the transfer function T{k; t), 
which connects the primordial spectrum P{k; tp) at time tp to the flnal spectrum at 
time tf 



where b{t) is the law of growth of perturbations in the linear regime. In the absence of 
other physical effects, this spectrum would simply scale with time in accord with the 
linear growth law for each perturbation mode. 

This processing of the primordial spectrum happens in a way that depends on 
cosmological parameters and the form of any non-baryonic dark matter, so that this 
function can be calculated for every model. 




(3.17) 



Averaging equation (3.15) over r gives 




(3.18) 




(3.19) 
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Biasing factor 

Attempts to confront theories of cosmological structure formation with observations 
of galaxy clustering are complicated by the uncertain and possibly biased relationship 
between galaxies and the distribution of gravitating matter. 

The problem is that one observes the distribution of galaxies, that is of luminous 
objects, which are luminous in a very specific way^, rather than the distribution of 
the total matter, while there is no a priori reason why the galaxy distribution should 
be a good tracer of the mass distribution of the Universe. Indeed, observations show 
that it definitely cannot be; the correlation functions for, to give an example, galaxies 
selected optically and galaxies selected in the infrared are different and hence clearly 
cannot both trace the mass distribution accurately. This effect is known as bias in 
the galaxy distribution, and it seriously impairs our abilities to use it to constrain the 
matter spectrum [51, 147, 166]. 

The literature contains at least three different definitions of the biasing factor or bias 
parameter, b, which are not equivalent to one another, but all represent the possible 
difference between mass statistics and the statistics of galaxy clustering. They are: 



where (t| represents the dimensionless variance in either galaxy counts or mass in 
spheres of radius Sh^^Mpc^ In general b need not be a constant. However, it may well 
be adequately represented by a constant across some range of scales, and indeed there 
is observational evidence supporting this as long as we look to large enough scales, 
which in practice more or less means scales in the linear regime [135]. 

In practice, b parametrizes our ignorance of galaxy formation in the same way as 
the mixing-length parameter does in the theory of stellar convection. To understand 
how this occurs we need to understand not only gravitational clustering, but also star 
formation and gas dynamics. All this complicated physics is supposed to be contained 
in the free parameter b. 



^We might imagine that galaxies should form not randomly sprinkled around according to the local 
density of matter, but at specific locations where collapse, cooling and star formation can occur. 

^This choice is motivated by the observational result that the variance of counts of galaxies in 
spheres of this size is of order unity, so that b ~ 1 /erg (mass). 




.galaxy — 



.mass -I 



(3.20) 
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3.3 Evolutionary equations 

To predict CMB anisotropics one has to solve the equations for the evolution of all 
particle species present (sec e.g. [32, 138, 177]). In this thesis we base on the work 
done by Ma & Bertschinger in 1995 [118], where they studied evolutionary equations 
in two different gauges (syncronous and Newtonian conformal). In particular, we focus 
on the syncronous gauge, briefly describing in the following sections the equations used 
for our work. 

Although the process of galaxy formation in recent epochs is well described by 
Newtonian gravity (and other physical processes such as hydrodynamics), a general 
relativistic treatment is required for perturbations on scales larger than the horizon 
size or before the horizon crossing time (see chapter 2). The use of general relativity 
brought in the issue of gauge freedom, which has caused some confusion over the 
years. Lifshitz (1946) [115] adopted the "synchronous gauge" for his coordinate system, 
which has since become the most commonly used gauge for cosmological perturbation 
theories.^ 

We consider only spatially flat background spacctimes with iscntropic scalar metric 
perturbations (sec §2.2 and §2.3). Here we give a complete discussion of CDM, baryons, 
photons and masslcss neutrinos in flat models and present the coupled, linearized 
Einstein, Boltzmann, and fluid equations for the metric and density perturbations. 

The CDM and the baryon components behave hke coUisionless and coUisional fluids, 
respectively, while the photons and the neutrinos require a phase-space description 
governed by the Boltzmann transport equation. We also derive analytically the time 
dependence of the perturbations on scales larger than the horizon. This information is 
needed in the initial conditions for the numerical integration of the evolution equations. 

The photon and neutrino distribution functions are expanded in Legendre polyno- 
mials, reducing the linearized Boltzmann equation to a set of coupled ordinary differ- 
ential equations for the expansion modes. 



^However, some complications associated with this gauge, such as the appearance of coordinate 
singularities and spurious gauge modes (see appendix A. 3), prompted Bardeen (1980) [11] and others 
(e.g. Kodama & Sasaki, 1984 [102]) to formulate alternative approaches that deal only with gauge- 
invariant quantities. 
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3.3.1 Phase space and the Boltzmann equation 

A phase space is described by six variables: three positions and their conjugate 
momenta Pj.^ The phase space distribution of the particles gives their number in a 
differential volume dx^dx^dx^dPidP2dP3 in phase space 

f{x\ Pj, r)dx^dx^dx^dPidP2dP3 = dN . (3.21) 

Importantly, / is a scalar and is invariant under canonical transformations. The zeroth- 
order phase space distribution is the Fermi-Dirac distribution for fermions (+ sign) and 
the Bose- Einstein distribution for bosons (— sign) 

where e = a(p^ + m^)^/^ = (P^ + a^m^)^/^ , Tq = aT denotes the temperature of the 
particles today, the factor gs is the number of spin degrees of freedom, and hp and /cb 
are the Planck and the Boltzmann constants. 

Following common practice (e.g.. Bond & Szalay 1983 [33]) we shall find it conve- 
nient to replace Pj by Qj = apj in order to eliminate the metric perturbations from the 
definition of the momenta. Moreover, we shall write the comoving 3-momentum qj in 
terms of its magnitude and direction: qj — quj where n^Ui = SijU^n^ = 1. Thus, we 
change our phase space variables, replacing f{x\Pj,T) by f{x'',q,nj,T). While this 
is not a canonical transformation (i.e., qi is not the momentum conjugate to x*), it is 
perfectly valid provided that we correctly transform the momenta in Hamilton's equa- 
tions. Note that we do not transform /. Because qj are not the conjugate momenta, 
d^xd^q is not the phase space volume element, and fd^xd^q is not the particle number. 

In the perturbed case we shall continue to define e as a(r) times the proper energy 
measured by a comoving observer, e = (g^-|-a^m^)^/^. For the models we are interested 
in, the photons and the massless neutrinos at the time of neutrino decoupling are 
all ultra-relativistic particles, so e in the unperturbed Fermi-Dirac and Bose-Einstein 
distributions can be simply replaced by the new variable q. 

The general expression for the energy-momentum tensor written in terms of the 
distribution function and the 4-momentum components is given by 

Tap^ J dP.dP^dP^i-g)-'/' ^f{x\Pj,T) , (3.23) 

^The conjugate momentuni is related to the proper momentum measured by an observer at a 
fixed spatial coordinate value by P, = a{Sij + ^hij)p' . 
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where g denotes the determinant of Qa/s- It is convenient to write the phase space 
distribution as a zeroth-order distribution plus a perturbed piece in the new variables 
q and rij 

f{x\ Pj, t) = fo{q) [l + q, nj, r)] . (3.24) 

In the synchronous gauge, from equation (3.23) we obtain to linear order in the per- 
turbations 



-t 

i 



-a-^ J q^dq ^ q^ + m?a'^ fo{q) (1 + ^) , 
a"^ J q^dq dQ q Ui fo{q) * , 

a-' [ q'dqdn -^4=T^ (1 + *) ■ 



(3.25) 



The phase space distribution evolves according to the Boltzmann equation, which can 
be written to first order in fc-space as 

1 fdf\ 



— + tl(]^.n)^ + -—^ 
OT 6 dmq 



V 



where 



Df df ^ dx^ df _^ dqdf ^ drii df 
dr dr dr dx'^ dr dq dr drii 



/o V^^. 



'df 



c 



dr 



(3.26) 



(3.27) 



c 



The terms in the Boltzmann equation depend on the direction of the momentum n 
only through its angle with k. Therefore, we shall assume that the initial momentum- 
dependence is axially symmetric, so that ^ depends on q = gn only through q and 
k • n. 



3.3.2 Cold dark matter 

CDM interacts with other particles only through gravity and can be treated as a 
pressureless perfect fluid. The CDM particles can be used to deflne the synchronous 
coordinates and therefore have zero peculiar velocities in this gauge. Setting 6 = a = 
and w = w = in equations (A. 32) leads to 

4 = , 4 = 0. (3.28) 

3.3.3 Massless neutrinos 

The energy density and the pressure for massless neutrinos are p,y = 3pi, — — T^o = T\. 
Prom equations (3.25) the unperturbed energy density pi, and pressure p^, are given by 

Piy = 3pt, = a^'^ J q^dqdVt qfo{q) , (3.29) 
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and the perturbations of energy density 5pi, and pressure 5pi, are 

5p^ = 35p^^ J q^dqdQqfoiq)-^ . (3.30) 

The Bohzmann equation simphfies for massless particles, for which e = q. To reduce 
the number of variables we integrate out the g-dependence in the neutrino distribution 
function and expand the angular dependence of the perturbation in a series of Legendre 
polynomials Pi{k ■ fi): 

F.(k, n, t) ^ ^ftr^ff"^ ^ f:(-0'(2i + T)Pi{k ■ n) . (3.31) 

Jq^dqqfoiq) 

In terms of the new variable -FV(k, n, r) and its harmonic expansion coefficients, the 
perturbations Si,, 9^, and (7,^, defined in eq. (A. 29), take the form 

^ ^ I ^^-^'^(^'"'''^) ^ -^'^o ' 



167r 



J dn(k- h)F„(k, h, r) = i , (3.32) 



F,(k,n,T) = iF,2 



Integrating equation (3.26) over (fdqqfQ{q) and dividing it by / q^dqqfo{q), the 
Boltzmann equation for massless neutrinos becomes 

OF 2 • 4 • 

+ ikpF,, = --h - -{h + 6fi)P2{iJ,) , (3.33) 

where ji = k ■ h and P2(a*) = (l/2)(3/x^ — 1) is the Legendre polynomial of degree 
2. Substituting the Legendre expansion for F^, and using the orthonormality of the 
Legendre polynomials and the recursion relation [l + l)Pi+i(//) = [21 + 1)ij,Pi{ijl) — 
lPi-i{li), we obtain the following equations: 

4 2 • 

Su — — Ou h , 

3 3 



i = - (7,) , 



8 3 4 8 

i^w = [/F.(z-i) - (/ + l)i^.a+i)] , ^ > 3 . (3.34) 

This set of equations governs the evolution of the phase space distribution of massless 
neutrinos. Note that a given mode Fi is coupled only to the (/ — 1) and (/ + 1) neigh- 
boring modes. In this way the Boltzmann equation (3.33) has been transformed into 
an infinite hierarchy of moment equations that must be truncated at some maximum 
multipole order Zmax (for more details see § 6.2). 
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3.3.4 Photons 

Photons evolve differently before and after recombination. Before recombination, pho- 
tons and baryons are tightly coupled, interacting mainly via Thomson scattering (and 
the electrostatic coupling of electrons and ions). After recombination, the Universe 
gradually becomes transparent to radiation and photons travel almost freely, although 
Thomson scattering continues to transfer energy and momentum between the photons 
and the matter. 

The evolution of the photon distribution function can be treated in a similar way as 
the massless neutrinos, with the exception that the collisional terms on the right-hand 
side of the Boltzmann equation are now present and they depend on polarization. We 
shall track both the sum (total intensity) and difference (Stokes parameter Q) of the 
phase space densities in the two polarization states, denoted respectively by F^(k, n, r), 
defined as in equation (3.31), and (^^(k, n, r). 

Expanding F^(k, n, r) and G^(k, n, r) in Legendre series as in equation (3.31), the 
collision operators can be rewritten as 



'dF\ 



Ai 



aUeCTT 



c 



1 1 \ °° 

(3.35)^ 
, (3.36) 



- {F,2 + G^o + G,2) (1 - P2) - E(-0'(2^ + 1)^7;^/ 

^ l>0 



where Ue is the proper mean density of the electrons. The left-hand-side of the Boltz- 
mann equation for F^ and G^ remain the same as for the massless neutrinos, so we 
obtain 
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1'Ft(i-i) — {I' + ^)Fj{i+i) — aneCTTF^i , />3, 



lG^{i_i) - (/ + l)G^(i+i) 
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-G^i + - {Fj2 + Gjo + G72) (sio + 



(3.37) 
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We truncate the photon Boltzmann equations in a manner similar to massless 
neutrinos, except that Thomson opacity terms must be added. For I — Zmax we replace 
equations (3.37) by 

Fji — /cF^(;_i) ~ aneCTrF-yi , 

Gji — —G^i — aneCTTG^i . (3.38) 

3.3.5 Baryons 

The baryons (and electrons) behave hke a non-relativistic fluid described, in the absence 
of coupling to radiation, by the energy-momentum conservation equation (A. 32) with 
5pi,/5pi) — — w <^ 1 and a — 0. Before recombination, however, the coupling of the 
baryons and the photons causes a transfer of momentum and energy between the two 
components. 

The momentum transfer into the photon component is represented by aneaT{Ob—0^) 
of equation (3.37). Momentum conservation in Thomson scattering then implies that 
a term (Ap^/Spb) aneCrriO^ — Ob) has to be added to the equation for 9i, (where we have 
used Pb <^ pb), so equation (A. 32) is modified to become 

Sb = -Ob - , 

Ob = -% + cySb+^aneaT{0^-9b) . (3.39) 
a Spb 

The square of the baryon sound speed is evaluated from 

Pb ji \ 6 ama J 

where is the mean molecular weight (including free electrons and all ions of H and 
He) and, in the second equality, we have neglected the slow time variation of iiJ 

3.3.6 Super-horizon-sized perturbations and initial conditions 

The evolution equations derived in the previous sections can be solved numerically once 
the initial perturbations are specified. We start the integration at early times when 
a given /c-mode is still outside the horizon, i.e., kr <^ 1, where kr is dimensionless. 



''This approximation is adequate because even during recombination, when jj, is largest, the baryons 
contribute very Uttle to the pressure of the photon-baryon fluid. 
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The behaviour of the density fluctuations on scales larger than the horizon is gauge- 
dependent. The fluctuations can appear as growing modes in one coordinate system 
and as constant modes in another. 

We are concerned only with the radiation-dominated era since the numerical in- 
tegration for all the /c-modes of interest will start in this era. At this early time the 
CDM and the baryons make a negligible contribution to the total energy density of 
the Universe: ptotai — Pv + P'y- The expansion rate is a/ a — t"^. We can analytically 
extract the time-dependence of the metric and density perturbations h, rj, 5, and 9 on 
super- horizon scales {kr <^ 1) from equations (3.34) and (3.37). The large Thomson 
damping terms in equations (3.37) drive the I > 2 moments of the photon distribution 
function F^i and the polarization function G^i to zero. Similarly, F,^; for Z > 3 can 
be ignored because they are smaller than F^2 by successive powers of kr. Equations 
(A.25), (A.27), (3.34), and (3.37) then give 

r^h + rh + 6[(1 - R^)5^ + RJ^] = , 

4 2- -1 

S-y + -9^ + -h^O, 9^- -k'^S^ = , 

O O 41: 

6u + ^9, + ^h = 0, 9, - ^k\5, - Aa,) = Q , (3.41) 

where we have defined R^, = Pu/{p-y + Pv)- For N^, flavors of neutrinos (A^^, = 3 in 
the standard model), after electron-positron pair annihilation and before the massive 
neutrinos become non relativistic, pu/ p-y — (7A^,^/8)(4/ll)'^/^ is a constant. 

To lowest order in kr, the terms oc k"^ in equations (3.41) can be dropped, and we 
have 9i, — 9^ — Then these equations can be combined into a single fourth-order 
equation for h 

whose four solutions are power laws: h oc t"' with n — 0, 1, 2, and —2. Prom equations 
(3.41) we also obtain 

h = A + B{kT)-^ + C{kTf + D{kT) , 
5 = (1 - R,)5^ + RJ, = -lB(kr)-^ -Icikrf -]:D{kT) , 

o o u 

9 = {1 - R,)9^ + RJ, = -^Dk, (3.43) 

8 

and A, B, C, and D are arbitrary dimensionless constants. The other metric pertur- 
bation T] can be found from equation (A.25) 

r] = 2C +^D{kT)-K (3.44) 
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Press & Vishniac (1980) [144] derived a general expression for the time dependence 
of the four eigenmodes. They showed that of these four modes, the first two (propor- 
tional to A and B) are gauge modes that can be eliminated by a suitable coordinate 
transformation. The latter two modes (proportional to C and D) correspond to phys- 
ical modes of density perturbations on scales larger than the Hubble distance in the 
radiation-dominated era. Both physical modes appear as growing modes in the syn- 
chronous gauge, but the C(A;t)^ mode dominates at later times. In fact, the mode 
proportional to D in the radiation-dominated era decays in the matter-dominated era 
(Ratra 1988 [146]). We choose our initial conditions so that only the fastest-growing 
physical mode is present (this is appropriate for perturbations created in the early 
Universe), in which case 9^ — 9^ — rj — to lowest order in kr. To get nonzero 
starting values we must use the full equations (3.41) to obtain higher order terms for 
these variables. To get the perturbations in the baryons we impose the condition of 
constant entropy per baryon. Using all of these inputs, we obtain the leading-order 
behaviour of super-horizon-sized perturbations in the synchronous gauge 
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Chapter 4 

The Mirror Universe 



4.1 Introduction to the mirror world 

In 1956 Li and Yang [112] proposed that the interactions of the fundamental particles 
could be non invariant under mirror reflection of the coordinate system (Parity trans- 
formations). Subsequently experiments confirmed that the weak interactions indeed 
violate the Parity. In particular they are left-chiral, i.e. they have V— A form. 

From the modern point of view, fundamental interactions (strong, weak and electro- 
magnetic) are described by the Standard Model (SM) based on the gauge symmetry 
SU{3) X SU{2) X U{1). The fermion fields quarks and leptons are ascribed to the 
certain representation of this symmetry. The electroweak symmetry SU{2) x U{1) is 
spontaneously broken down to electromagnetic U{l)em at the scale v ~ 100 GeV, and 
as a result, the gauge bosons and Z acquire masses ~ 100 GeV. Parity violation 
in weak interactions is related to the fact that the left-handed components of quarks 
Qi = {u,d)iL and leptons k = {i',e)iL transform as doublets of SU{2) x U{1), while 
the right-handed components Um, dm and em are the singlets {i = 1,2,3 is the family 
index). Thus, for particles one observes that their weak interactions have V— A form. 

However, the fact that the observable matter is constituted by quarks and leptons 
and not by their antiparticles is a consequence of the baryon asymmetry of our Universe. 
If instead the baryon asymmetry would have the opposite sign, then the observed 
Universe would be made of antiparticles, namely anti-quarks qi = Clji^ and anti- 

- — T 

leptons li = Cli , which are now right-handed, and so we would observe the V-l-A 
form in their weak interactions, i.e. our world would be the right-handed.^ 

^In fact, one cannot exclude the possibility that in some patch of the Universe, very far from us, 
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The fact that we have particle excess over antiparticles is probably related to the CP 
violating features of the gauge and particle sector of our world. So, in the context of the 
full theory (Grand Unification or seesaw model of neutrino masses or Supersymmetric 
Standard Model) the sign of the baryon asymmetry of the Universe is defined by the 
CP violation phases in the corresponding baryogenesis mechanism (respectively GUT 
baryogenesis, leptogenesis or electroweak baryogenesis). 

In short, for creating the matter in the Universe, one should have to violate not 
only C and P, but also the CP symmetry, which actually should be identified as the 
particle- antiparticle symmetry. 

Today, it is widely believed that mirror symmetry is in fact violated in nature. The 
weak nuclear interaction is the culprit, with the asymmetry being particularly striking 
for the weakly interacting neutrinos. For example, today we know that neutrinos only 
spin with one orientation. Nobody has ever seen a right-handed neutrino. 

If one wants to retain the mirror symmetry, the only possibility is a complete 
doubling of the number of particles. There is an old idea suggested by Li and Yang 
[112] that there can exist a hidden mirror sector of particles and interactions which is 
the exact duplicate of our visible world, i.e. with the Lagrangian exactly identical to 
the Lagrangian of the ordinary particles, however with right-handed interactions. In 
other words, for each type of particle, such as electron, proton and photon, there is 
a mirror twin, so that the ordinary particles favor the left hand, the mirror particles 
favor the right hand. If such particles exist in nature, then mirror symmetry would be 
exactly conserved. In fact one could have the Parity symmetry as an exact symmetry 
of exchange of the ordinary and mirror particles. 

The mirror particles can exist without violating any known experiment. Thus, the 
correct statement is that the experiments have only shown that the interactions of the 
known particles are not mirror symmetric, they have not demonstrated that mirror 
symmetry is broken in nature. 

What really needs to be done is to understand the experimental imphcations of the 
existence of mirror particles and find out whether such things could describe our Uni- 
verse. In particular, it is exactly the aim of this thesis to investigate the cosmological 
implications of the mirror world. 

The hypothesis of the mirror sector has attracted a significant interest over last 
years, in particular being motivated by the problems of neutrino physics [1, 23, 71, 
76, 178], gravitational microlensing [16, 21, 26, 67, 123], gamma ray bursts [27, 179], 



the baryon asymmetry has the opposite sign and thus antiparticles rather than particles are dominant. 



§4.1 Introduction to the mirror world 



51 



ultra-high energy cosmic rays [24], flavour and CP violation [17, 22, 152], etc. The 
basic concept is to have a theory given by the product G x G" of two identical gauge 
factors with the identical particle contents, which could naturally emerge, e.g., in the 
context of Eg x E'^ superstring. Two sectors communicate through gravity and per- 
haps also via some other messengers. A discrete symmetry P{G <-> G') interchanging 
corresponding fields of G and G", so called mirror parity, imphes that both particle 
sectors are described by the same Lagrangians.^ 

In particular, one can consider a minimal symmetry Gsm x G^sm; where Gsm — 
SU{3) X SU{2) X U{1) stands for the standard model of observable particles: three 
families of quarks and leptons qi, Ui, df, li, Ci {i — 1,2,3) and the Higgs doublet 0, 
while G'gj^ — [SU{3) x SU{2) x U{1)]' is its mirror gauge counterpart with analogous 
particle content: fermions q^, u'^, and the Higgs 0'. (Prom now on all fields 

and quantities of the mirror (M) sector will have an apex to distinguish from the ones 
belonging to the observable or ordinary (O) world.) The mirror parity implies that all 
coupling constants (gauge, Yukawa, Higgs) have the same pattern in both sectors and 
thus their microphysics is the same.^ 

According to this theory the mirror partners have the same mass as their ordinary 
counterparts, but the mirror particles interact with ordinary particles predominately 
by gravity only. The three non-gravitational forces act on ordinary and mirror parti- 
cles completely separately (and with opposite handedness: where the ordinary particles 
are left-handed, the mirror particles are right-handed). Por example, while ordinary 
photons interact with ordinary matter, they do not interact with mirror matter^. Sim- 
ilarly, the "mirror image" of this statement must also hold, that is, the mirror photon 
interacts with mirror matter but does not interact with ordinary matter. The upshot 
is that we cannot see mirror photons because we are made of ordinary matter. The 
mirror photons would simply pass right through us without interacting at all! 

The mirror symmetry requires that the mirror photons interact with mirror elec- 

^In the brane world picture, the M sector can be the same O world realized on a parallel brane, 
G' = G [6]. 

^The mirror parity could be spontaneously broken and the weak interaction scales {(f) = v and 
{(j)') = v' could be different, which leads to somewhat different particle physics in the mirror sector 
[1, 16, 21, 23, 123, 124, 125]. In this thesis we treat only the simplest case, v = v' , in which the M 
sector has exactly the same physics as the O one. A possibility for further studies is to add another 
free parameter related to the difference of the weak scales in the two sectors, but at the moment it 
would be a useless complication. 

^On the quantum level, small new fundamental interactions connecting ordinary and mirror par- 
ticles are possible. Various theoretical constraints suggest only a few possible types of interactions: 
neutrino-mirror neutrino mass mixing and photon-mirror photon kinetic mixing [70, 71, 101]. For the 
purposes of this thesis, this possibility is ignored. 
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trons and mirror protons in exactly the same way in which ordinary photons interact 
with ordinary electrons and ordinary protons. A direct consequence of this is that a 
mirror atom made of mirror electrons and a mirror nucleus, composed of mirror protons 
and mirror neutrons, can exist. In fact, mirror matter made of mirror atoms would 
also exist with exactly the same internal properties as ordinary matter, but would be 
completely invisible to us! 

Today the most important aspect of the mirror scenario is that it predicts the 
existence of dark matter in the Universe in a very natural manner. In fact, mirror 
matter would be invisible, making its presence felt by its gravitational effects, which 
is exactly the definition of "dark matter" ! 

One could naively think that due to mirror parity the O and M particles have not 
only the exactly identical particle physics (microphysics) , but also the same cosmology 
at all stages of the Universe evolution. However, identical Lagrangians of two particle 
sectors do not necessarily imply that the ordinary and mirror worlds should be realized 
in the same initial states.^ In fact, should the O and M particles have the same 
cosmological densities, this would be in the immediate conflict with the Big Bang 
nucleosynthesis (BBN) bounds on the effective number of extra hght neutrinos, AAT^ < 
1 [116]: the mirror photons, electrons and neutrinos would give a contribution to the 
Hubble expansion rate exactly equal to the contribution of the ordinary ones, equivalent 
to AA^,^ ~ 6.14. Therefore, to satisfy the nucleosynthesis limit, the M particles density 
in the early Universe should be appropriately reduced, i.e. at the BBN epoch they 
should have a temperature smaller than the ordinary one, T' < T. Namely, the 
conservative bound AA^,^ < 1 implies that T'/T < 0.64. 

This situation is plausible if two following conditions are satisfied: 

A. At the initial moment two systems are born with different densities. In partic- 
ular, the inflationary reheating temperature in the M sector should be lower than in 
the visible one, < Tr, which can be achieved in certain models [16, 21, 24, 103]. 

B. The M and O particles interact very weakly, so that two systems do not come 
into the thermal equilibrium with each other in the early Universe. This condition is 
automatically fulfilled if two worlds communicate only via the gravity.^ 

If two sectors have different reheating temperatures, and if they do not come into 
the thermal contact at later stages, then during the Universe expansion they evolve 

^For analogy two oscillators with the same Lagrangians could naturally have different amplitudes. 

^More generally, there could be other messengers like superheavy gauge singlet fields or light 
singlets of the moduli type. In either case, they should mediate the effective couplings between the 
O and M particles suppressed by a large mass factor M ~ Mp or so. 
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independently and approach the BBN epoch with different temperatures. 

Several possibilities were discussed in the literature when after inflation two particle 
sectors could be reheated at different temperatures. For example one can consider the 
chaotic inflation scenario with two inflaton fields, the ordinary one 77 and the mirror 
one 77', having identical harmonic potentials /i^r]^ + /i^r)'^ and interacting respectively 
with the O and M particles. In this case one can show that different initial amplitudes 
of rj and 77' would result in the different temperatures of the O and M thermal bathes 
born after inflaton decay [24, 90]. The thermal energies of two sectors simply reflect 
the difference of the initial energies of two oscillators rj and 77'. This naive scenario, 
however, is not valid for inflaton with unharmonic potentials A77^ + A77'^ [24] . 

An attractive realization of the inflationary paradigm is provided by supersym- 
metric models of hybrid inflation. The simplest model is based on the superpotential 
W — AyS'($^ — /i^) containing the inflaton field S and the additional "orthogonal" field 

where A is order 1 coupling constant and // is a dimensional parameter of the or- 
der of the GUT scale {fx ~ 10^^ GeV) [54]. The supersymmetric vacuum is located 
a,t S — 0, ^ — fjL, while for the field values ^ — 0, S > fx the tree level potential 
has a flat valley with an energy density V — A^//^. Since the supersymmetry is bro- 
ken by the non- vanishing F-term, Fg — A//^, the flat direction is lifted by radiative 
corrections and the potential of S gets a slope which is appropriate for the slow roll 
conditions. The COBE results on the CMB anisotropy (see previous chapter) imply 
that y^/^ ~ e^/^ X 7 • 10^^ GeV, where e <IC 1 is a slow-roll parameter. 

In ref. [19] there was suggested to modify the superpotential to the following form 

W^preheat = \S {^^ - /X^) + Q^^^ + ^'^) , (4.1) 

where and ip' are some additional superfields. Such a modification was motivated by 
the necessity to solve the initial conditions problem. Fast damping of $ is a pretext 
of inflationary stage which allows the inflaton energy density ~ fi^ to dominate and 
then S can slowly roll to the origin. This function is carried by the second term in 
(4.1) - the oscillating orthogonal field $ lastly decays into and particles which 
have practically no contact to inflaton S. In addition, with vanishing also effective 
mass terms of ^''s disappear and the latter fields start behaving as massless - they stop 
oscillating and freeze. In general, oscillations ^ and ^' freeze at different amplitudes, 
typically ~ /i, at which they are catched by the moment when their mass drops 
below the Hubble parameter (see fig. 3 of ref. [19]). When slow-roll ends up, all 
fields start oscillating around their vacuum values, 5 = 0, $ = ^, ^' = 0, and 
reheat the Universe e.g. through the decays of ^ ^ 0i02 and ^' 0i02) due to 
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superpotential terms ^^0102 and ^'0i02> where 0i^2 and 0^ 2 Higgs doublet 

superfields respectively for the O and M sectors. Then different magnitudes of ^ and 
^' at the end of slow-roll phase should reflect into difference of reheating temperatures 
Tr and in two systems, simply because of the energy difference stored into ^ and 
^' oscillations which decay respectively into O and M Higgses. ^ 

We have also to make sure that after reheating two sectors do not come into the 
thermal equilibrium to each other, or in other words, that interactions between O and 
M particles are properly suppressed. In particular, the quartic interaction between the 
O and M Higgs scalars A(0V)(0'V') are dangerous, since they would bring to sectors in 
thermal equilibrium unless the coupling constant is very small, A < 10~^ [21]. Another 
dangerous coupling can be presented by the kinetic mixing of the O and M photons, 
aF'^'^F^^, which can be safe only if a < 10~^. 

The operator which can hnk O and M sectors, has a higher dimension 

^{k4>)m + h.c. , (4.2) 

and is cutoff by a large mass factor M ~ Mp or so, and thus are safe. This operator 
creates mixing between the O and M neutrinos and so they could oscillate in each 
other. From the point of view of the ordinary observer, the mirror neutrinos should 
appear as the sterile ones. 

In a supersymmetric version of the theory the contact between two sectors is nat- 
urally weak. In particular, the dangerous mixed term A(0V)(0'V') is not allowed 
anymore, and the lowest order operator between the Higgses of two sectors in the 
superpotential has a dimension 5 

^{M2)m'2) , (4.3) 

which are suppressed by a large mass factor M as the operators (4.2), and 
thus are safe. The same holds true for soft supersymmetry breaking terms like 
("^3/2/-^)(0i02)(0i02), etc., where 771,3/2 ~ 100 GeV is the gravitino mass. 

As for the kynetic mixing term F^'^F'^^^ between the fleld-strength tensors of the 
gauge factors U{1) and C/(l)', it can be forbidden by embedding Gsm x G'^m the 
grand unifled group like SU(5) x SU{5y or 50(10) x 50(10)'. 



For other scenarios of the asymmetric reheating see [16, 21, 24, 103, 123, 124, 125]. 
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4.2 Thermodynamics of the Mirror Universe 

Once the O and M systems arc decoupled aheady after reheating, at later times t they 
will have different temperatures T{t) and T'{t), and so different energy and entropy 
densities 

Pit) = ^9*iT)T' , At) = , (4.4) 

<t) = -^9s{T)T^ , s'{t) = —gmr' . (4.5) 

The factors g^, gg and g'^, g'^ account for the effective number of the degrees of freedom 
in two systems, and in general can be different from each other. During the Universe 
expansion, the two sectors evolve with separately conserved entropies. Therefore, the 
ratio X = {s'/sY^^ is time independent^, while the ratio of the temperatures in two 
sectors is simply given by 



T'jt) _ 



g'siT') 



(4.6) 



The Hubble expansion rate is determined by the total energy density p — p -\- p' , 
nP- Therefore, at a given time i in a radiation dominated epoch we 

have 

H(t) = i = 1.66^|1 = 1.66VS(n£ (4.7) 

in terms of O and M temperatures T{t) and T'{t), where 

UT) = 9.{T){1 + ax') , -giiT') = g^T') (l + A) • (4.8) 

Here the factor a(T,T') = {g'^T') I g^{T)\ ■ {g^ij)! g'XT')fl^ takes into account that for 
T' the relativistic particle contents of the two worlds can be different. However, 
except for very small values of x, we have a ~ 1. So hereafter we always take g^^iT) — 
g^{T){l + x^) and g'^{T') = gl{T'){l + x"'). In particular, in the modern Universe we 
have a{To, T^) = 1, ^^(ro) = g'^{T^) = 3.91, and x = Tj^/Tq, where Tq, are the present 
temperatures of the O and M rehc photons.^ 



^We assume that expansion goes adiabatically in both sectors and neglect the additional entropy 
production due to the possible weakly first order electroweak or QCD phase transitions. 

''The frozen ratio of the neutrino and photon temperatures in the M sector Tq = T^^q/Tq depends 
on the v' decoupling temperature from the mirror plasma, which scales approximatively as ~ 
x-^/^Td, where = 2 — 3 MeV is the decoupling temperature of the usual neutrinos. Therefore, 
unless X < 10~^, has a standard value ro = Ti,a/To = (4/11)^/'^. For x < lO"'', becomes 
larger than the QCD scale A ~ 200 MeV, so that due to the mirror gluons and light quarks u',d',s' 
contribution we would obtain = (4/53)^3, g'^[T^) = 2.39 and so T^/Tq ~ 1.2 a;. However, in the 
following such small values of x are not of our interest. 
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It is useful to note that, due to the difference in the initial temperature conditions 
in the two sectors, reactions at the same temperature T* = Tl occur at different 
times t'^ — xH^, which implies different rates of the Hubble expansion (in particular 
H{t'^) > H{t^), given x < 1), while reactions at the same time = t'^ occur at different 
temperatures Tl — xT^. 

In fact, X is the only free parameter in our model, and in general it could be 
constrained by the BBN bounds. The observed abundances of hght elements are in 
good agreement with the standard nucleosynthesis predictions, when at T ~ 1 MeV 
we have — 10.75 as it is saturated by photons 7, electrons e and three neutrino 
species i'cij^t- The contribution of mirror particles (7', e' and i^e./^.r) would change it to 
g* — g*{l + x^). Deviations from g^ — 10.75 are usually parametrized in terms of the 
effective number of extra neutrino species, Ag — g* — 10.75 = 1.75AA^,^. Thus we have 

AN^ = 6.14 • x^ . (4.9) 

In view of the present observational situation, a rehable bound is AA^,^ < 1 [19], which 
translates as x < 0.64. This hmit very weakly depends on AA^,^; e.g., AA^,^ < 1.5 
imphes x < 0.70. 

As far as ^ 1, in a relativistic epoch the Hubble expansion rate (4.7) is dom- 
inated by the O matter density and the presence of M sector practically does not 
affect the standard cosmology of the early ordinary Universe. However, even if the 
two sectors have the same microphysics, the cosmology of the early mirror world can 
be very different from the standard one as far as the crucial epochs like baryogenesis, 
nucleosynthesis, etc. are concerned. Any of these epochs is related to an instant when 
the rate of the relevant particle process r(T), which is generically a function of the 
temperature, becomes equal to the Hubble expansion rate H{T) (see chapter 1). Ob- 
viously, in the M sector these events take place earlier than in the O sector, and as a 
rule, the relevant processes in the former freeze out at larger temperatures than in the 
latter. 

In the matter domination epoch the situation becomes different. In particular, we 
know that ordinary baryons can provide only a minor fraction of the present cosmolog- 
ical density, whereas the observational data indicate the presence of dark matter. So, 
it is interesting to question whether the missing matter density of the Universe could 
be due to mirror baryons. 
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4.3 Baryogenesis features in ordinary and mirror 
sectors 

As we expressed in §1.6, in general there are two kinds of baryogenesis scenarios: the 
electroweak baryogenesis and the GUT baryogenesis or leptogenesis. 

Whatever the mechanism responsible for the observed baryon asymmetry (BA) 
in the O sector 77 = n^/n^, it is most hkely to think that the BA in the M sector 
rj' — n'^/n'^ is produced by the same mechanism, and moreover the rates of the B and 
CP violation processes are parametrically the same in both cases. However, the out of 
equilibrium conditions should be different since at relevant temperatures the Universe 
expansion is faster for the M sector. Below we show that by this reason 77' typically 
emerges larger than 77 for either type (a) or (b) scenarios. 

The M baryons can be of the cosmological relevance if exceeds — 3.66 x 
10^?7/i~^ = 0.04 — 0.05, whereas > 1 would overdose the Universe. So we are 
interested in a situation when the ratio 

falls in the range from 1 to few tens. Since n'^ = x^n^, we obtain f3 = x^rj'/r]. Therefore, 
f]' > 1] does not a priori mean that (3 > 1, and in fact there is a lower limit x > 10~^ 
or so for the relevant parameter space. Indeed, it arises from x^ = firj/rj' by recalling 
that T] ~ 10^^, while rj' can be taken at most ~ 10^"^, the biggest value which can be 
principally realized in any baryogenesis scheme under the realistic assumptions. 

GUT Baryogenesis 

The GUT baryogenesis mechanism is typically based on a superheavy boson X under- 
going the B and CP violating decays into quarks and leptons. The following reaction 
rates are of relevance: 

Decay: Td ~ axMx for T < Mx oiVd ^ axM^jT for T > Mx, where olx is the 
coupling strength of X to fermions and Mx is its mass; 

Inverse decay: Fjo ~ ToiMx/Tf^^ ex.-p{-Mx/T) for T < or Fid ~ Pd for 

T>Mx; 

The X boson mediated 2^2 processes: P5 ~ rixcr ~ ^Oi'x'^^/ i-^x + T^V, where the 
factor A amounts for the possible reaction channels. 

The final BA depends on a temperature at which X bosons go out from equilibrium. 
One can introduce a parameter which measures the effectiveness of the decay at the 
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epoch T Mx [104]: k = {Vd/'2H)t=Mx = 0.3^* ' {axMpi/Mx). The larger is k the 
longer equilibrium is maintained and the freeze-out abundance of X boson becomes 
smaller. Hence, the resulting baryon number to entropy ratio, B — ub/s ~ O.I477 is a 
decreasing function of k. It is approximately given as S ~ {e/gs)F{k,kc), where e is 
the CP violating factor and 

r 1 if A; < 1 

F{k, = 0.3A;-i(log k)-^-^ iil<k<k, (4.11) 

Here kc is a critical value defined by equation A;c(log/cc)~^'^ = 'iQQ/{Aax)- It distin- 
guishes between the regimes k < kc, in which inverse decay is relevant, and k > kc, in 
which instead 2 -^-^ 2 processes are the dominant reason for baryon damping. 

In a general context, without referring to a particular model, it is difficult to decide 
which range of parameters k and kc can be relevant for baryogenesis. One can impose 
only the most reasonable constraints gs{T — Mx) > 100 and e < 10"^, and thus e/ Qs < 
10~^ or so. For a given mechanism responsible for the observed baryon asymmetry 
B ~ 10~^°, this translates into a lower bound F{k, kc) > 10"^. 




Figure 4.1: Panel A. The combination F{kx^ ,kc) as a function of k for kc = 10^ and 
X = 0.6, 0.1, 0.01 (dash line). The solid curve corresponding to a; = 1 in fact measures the 
possible BA in the ordinary world, F{k, kc) = {qs/ ^)B{k). Panel B. The curves confining the 
parameter region in which (3 = Q.'^/Vt\y varies from 1 to 100, for x = 0.6 (solid line) and for 
X = 0.01 (dash). The parameter area above thick solid curve corresponds to F{k, kc) < 10~^ 
and it is excluded by the observable value of 77. 

The presence of the mirror sector practically does not alter the ordinary baryogen- 
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esis. The effective particle number is g*{T) — g^{T){l + x^) and thus the contribution 
of M particles to the Hubble constant at T ~ Mx is suppressed by a factor (which 
is very small, if we consider the BBN limits explained in § 4.2). 

In the mirror sector everything should occur in a similar way, apart from the fact 
that now at T' ~ Mx the Hubble constant is not dominated by the mirror species 
but by ordinary ones: g'^{T') ~ g'^{T'){l + x~^). As a consequence, we have k' — 
{T'j^/2H)\T'=Mx — k^^- Since the value of kc is the same in the two sectors, the mirror 
baryon asymmetry can be simply obtained by replacing k ^ k' — kx^ in eq. (4.11), 
i.e. B' — n's/s' ~ {elg'^F{k' — kx'^,k^. Since F is a decreasing function of k, then 
for X < 1 we have F{kx^, kc) > F{k, kc) and thus we conclude that the mirror world 
always gets a larger BA than the visible one, B' > B. 

However, this does not a priori mean that is always larger than Qb- Since the 
entropy densities are related as s'/s — x^, for the ratio /3 = ^{^/^h we have 

^^^^"ns ~ Bs F{k,kc) ■ ^ ^ 

The behaviour of the factor x^F{kx'^. kc) as a function of k for different values of the 
parameter x is given in the fig. 4.1 A. Clearly, in order to have fl'^ > fib the function 
F{k, kc) have to decrease faster than between k' — kx"^ and k. Closer inspection 
of the function (4.11) reveals that the M baryons can be overproduced only if k is order 
kc or larger. In other words, the relevant interactions in the observable sector maintain 
equilibrium longer than in the mirror one, and thus ordinary BA can be suppressed 
by an exponential Boltzmann factor, while the mirror BA could be produced still in 
non-exponential regime k' < kc- 

In fig. 4.1 B we show the parameter region in which (3 = Q'^/Q]^ falls in the range 
1 — 100, in confront to the parameter area excluded by condition F{k, kc) > 10~^. We 
see that for x = 0.6 there is an allowed parameter space in which f3 can reach values 
up to 10, but P — 100 is excluded. For a limiting case x = 10~^, as it was expected, 
the parameter space for > 1 becomes incompatible with F{k, kc) > 10~^. For 
intermediate values of x, say a; ~ 0.1 — 0.3, also the values (3 ~ 100 can be compatible. 

The above considerations can be applied also in the context of leptogenesis. One 
should remark, however, that potentially both the GUT baryogenesis or leptogenesis 
scenarios are in conflict with the supersymmetric inflation scenarios, because of the 
upper limit on the reheating temperatures about Tr < 10^ GeV from the thermal 
production of gravitinos [59]. Moreover, it was shown recently that the non-thermal 
gravitino production can impose much stronger limits, Tr < 10^ GeV or so [64, 84]. 
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This problem can be fully avoided in the electroweak baryogenesis scenario, which 
instead is definitely favoured by the supersymmetry. 

Electroweak Baryogenesis 

The electroweak (EW) baryogenesis mechanism is based on the anomalous B- violating 
processes induced by the sphalerons, which are quite rapid at high temperatures, but 
become much slower when temperature drops below 100 GeV. A successfuU scenario 
needs the first order EW phase transition and sufficient amount of CP violation, which 
conditions can be satisfied in the frames of the supersymmetric standard model, for 
certain parameter ranges [149]. 

The characteristic temperature scales of the electroweak phase transition are fixed 
entirely by the form of the finite temperature effective potential 



where all parameters can be expressed in terms of the fundamental couplings in the 
Lagrangian. For large temperatures, T » 100 GeV, the electroweak symmetry is 
restored and V{(f),T) has a minimum at = 0. With the Universe expansion the 
temperature drops, approaching the specific values which define the sequence of the 
phase transition. These are all in the 100 GeV range and ordered as Ti > Tc > Tt, > ^o- 
Namely, below T — Ti the potential gets a second local minimum (f)+{T). At the 
critical temperature T — Tc the latter becomes degenerate with the false vacuum 
= 0. At temperatures T < Tc to the true vacuum state (f) — (t)+{T) becomes 
energetically favoured, and transition to this state can occur via thermal quantum 
tunneling, through the nucleation of the bubbles which then expand lastly, percolate 
and finally fill the whole space within a horizon by the true vacuum. 

The bubble production starts when the free energy barrier separating the two min- 
ima becomes small enough. The bubble nucleation temperature Tj, is defined as a 
temperature at which the probability for a single bubble to be nucleated within a 
horizon volume becomes order 1 



where Fc{T) is the free energy and uj is an order 1 coefficient [149]. In particular, in 
the limit of thin wall approximation we have 



V{(l), T) = D(T^- T^)(P^ - ETcj)^ + At0^ , 



(4.13) 




(4.14) 




(4.15) 
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The condition (4.14) results into large values of Fc(Tb)/Tb, typically order 10^. 

Once the bubble nucleation and expansion rate is larger than the Hubble parameter, 
the out of equilibrium condition for anomalous B violating processes is provided by 
the fast phase transition itself. The BA can be produced inside the bubbles due to the 
CP violation since the quarks and antiquarks have different reflection coefficients on 
the bubble wall. The baryogenesis rate is completely independent from the Universe 
expansion and it occurs practically in one instant as compared to the cosmological 
time scale of this epoch. 

As for the mirror sector, it is described by the same Lagrangian as the ordinary one 
and so the effective thermal potential of the mirror Higgs V{(j)', T') has the same form 
as (4.13). Then the temperature scales, which are defined entirely by the form of the 
effective potential, should be exactly the same for O and M sectors. Namely, T[ — Ti 
and T;' = Tc- 

The equation (4.14) is the same for the M sector apart from the fact that the 
corresponding Hubble constant is different: H{T' — — x~'^H{T — T^. Therefore, 
we obtain 



+ 81ogx, (4.16) 



which in turn tells that the bubble nucleation temperatures in two sectors are practi- 
cally equal: — Tb( 1 + 0.01 logx). (Clearly, the phase transition in M sector occurs at 
earlier time than in O sector: t'^ ~ a;~^tb-) The reason is that between the temperature 
scales Tc and Tq the free energy FdT) is a rapidly changing function but the change in 
temperature itself is insignificant. Hence, we expect that the initial BA's produced right 
at the phase transition to be the same in O and M sectors: B{T — T^) — B'{T' — Tb). 

However, the instantly produced baryon number can be still washed out by the 
sphaleron interactions. The anomalous B violation rate T{T) ~ exp[— F(T')/T'], where 
F{T) is the sphaleron free energy at finite T, may be large enough inside the bubble 
as far as the temperature is large. But it quickly falls as the temperature decreases, 
and baryon number freezes out as soon as r{t) drops below H{t), given by eq. (4.7). 
The wash-out equation dB/dt — —r{t)B can be rewritten as 

= (4.17) 

B HT ^ ^ 

and integrated. Then the final BA in the O and M sectors can be expressed respectively 
as 

B = B{T^)D^'+^"^~"'" , B' = B{T^)D^"^'+^"^~"'" , (4.18) 
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where D < 1 is the baryon number depletion factor 



D = exp 



PI 1^ dT^ 



(4.19) 



and g^, ~ O(IOO) in the supersymmetric standard model. Thus we always have B' > B, 
while the M baryon mass density relative to O baryons reads 



O' R' 1 — 

(3(x) = )^ = x'^ = x'D-^^^\ K{x) 



B 



(4.20) 




Figure 4.2: The contours of /? = ri[,/r2b in the plane of the parameters x and D, correspond- 
ing to /3 = 1, 10 and 100 from top to bottom. 

Lacking a precise theory for non-perturbative sphaleron transitions in the broken 
phase, the exact value of D cannot be calculated even in the context of concrete 
models. If D ~ 1, the wash-out is ineffective and practically all BA produced right 
at the bubble nucleation is conserved. In this case Q'^ should be smaller than fib- 
However, if D is enough small, one can achieve sufficiently large Q'^. The contour plot 
for the parameters x and D for which (3 falls in the range 1 — 100 is given in fig. (4.2). 
For small x we have essentially f3 ^ x^D~^ and thus 100 > f3 > 1 requires a depletion 
factor in the interval D = (10^^ — l)x^. Once again, for x ~ 10~^ one needs the 
marginal values D ~ 10^*^ — 10^^ below which the observable BA B ~ 10~^° cannot 
be produced at all. 
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4.4 Baryogenesis via particle exchange between or- 
dinary and mirror sectors 

Long time ago Sakharov [153] has suggested that a non-zero baryon asymmetry can 
be produced in the initially baryon symmetric Universe if three conditions are fulfilled: 
B-violation, C- and CP- violation and departure from thermal equilibrium. These con- 
ditions can be satisfied in the decays of heavy particles of grand unified theories. On 
the other hand, the non perturbative sphaleron processes, which violate B + L but 
conserve B — L, are effective at temperatures from about 10^^ GeV down to 100 GeV 
[108]. Thus, one actually needs to produce a non- zero B — L rather than just B, a fact 
that disfavors the simplest baryogenesis picture based on grand unification models like 
SU (5). When sphalerons are in equilibrium, the baryon number and B — L are related 
as B = a{B — L), where a is a model dependent order one coefficient [149]. Hence, 
the observed baryon to entropy density ratio, B — ub/s — (0.6 — 1) x 10~^°, needs to 
produce B - L 10"^°. 

The seesaw mechanism for neutrino masses offers an elegant possibility of generating 
non-zero B — L in CP- violating decays of heavy Majorana neutrinos N into leptons and 
Higgses, the so called leptogenesis scenario [81, 82]. Namely, due to complex Yukawa 
constants, the decay rates r{N ^ Icj)) and r(iV — Icj)) can be different from each other, 
so that leptons I and anti-leptons I are produced in different amounts. 

An alternative mechanism of leptogenesis based on the scattering processes rather 
than on decay was suggested in ref. [14]. The main idea consists in the following. 
There exists some hidden (shadow) sector of new particles which are not in thermal 
equilibrium with the ordinary particle world as far as the two systems interact very 
weakly, e.g. if they only communicate via gravity. However, other messengers may well 
exist, namely, superheavy gauge singlets like right-handed neutrinos which can mediate 
very weak effective interactions between the ordinary and hidden leptons. Then, a net 
B — L could emerge in the Universe as a result of CP- violating effects in the unbalanced 
production of hidden particles from ordinary particle collisions. 

Here we consider the case when the hidden sector is a mirror one. As far as the 
leptogenesis is concerned, we concentrate only on the lepton sector of both O and M 
worlds. Therefore we consider the standard model, containing among other particles 
species, the lepton doublets li — {p, e), (i = 1,2,3 is the family index) and the Higgs 
doublet 0. Then the mirror standard model contains the lepton doublets /■ = {u', e')i 
and the Higgs doublet (p', so that the products Zj0 and l^cf)' are gauge invariant in both 
sectors. 
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In the spirit of the seesaw mechanism, one can introduce some number of heavy 
fermions A^o(a = l,2,..), called heavy neutrinos, which are gauge singlets and thus can 
couple to Z, as well as to l',(f)'. In this way, they play the role of messengers between 
ordinary and mirror particles. The relevant Yukawa couplings have the form 

h,ahNa(t> + KjiNa(l)' + ^Ma^NaN, + h.C. (4.21) 

(charge-conjugation matrix C is omitted); all fermion states N, I' are taken to be left- 
handed while their C-conjugate, right-handed anti-particles are denoted as /, N, I'. It 
is convenient to present the heavy neutrino mass matrix as Mab = gabM, M being the 
overall mass scale and gab some typical Yukawa constants. (Without loss of generality, 
Qab can be taken diagonal and real.) 

Due to mirror Parity, both particle sectors are described by identical Lagrangians, 
that is, all coupling constants (gauge, Yukawa, Higgs) have the same pattern in both 
sectors and thus their microphysics is the same. Namely, the M Parity as a discrete 
symmetry under the exchange d> ^ (t)'\ I ^ I', etc., implies h'^^ = h*^^. 

Integrating out the heavy states in the couplings (4.21), we get the effective 
operators 

^hlj<P<P + ^kW + ^W<t>' + h.c. , (4.22) 

with couphng constant matrices of the form A — hg~^hF , A' — h'g'^h!^ and D — 
hg~^h'^ . Thus, the first operator in eq. (4.22), due to the ordinary Higgs vacuum 
expectation value (VEV) {(j)) — v ^ 100 GeV, induces the small Majorana masses of 
the ordinary (active) neutrinos. In addition, if the mirror Higgs 0' also has a non- 
zero VEV {(f)') — v' <^ M, then the third operator provides the masses of the mirror 
neutrinos contained in I' (which in fact are sterile for the ordinary observer). And 
finally, as we see, the operator (4.2) is also induced by the heavy neutrino exchanges 
(the second term in (4.22), which gives rise to the mixing mass terms between the 
active and sterile neutrinos). The total mass matrix of neutrinos i/ C I and u' C I' 
reads as [23] 

Thus, this model provides a simple explanation of why sterile neutrinos could be light 
(on the same grounds as the active neutrinos) and could have significant mixing with 
the ordinary neutrinos. For example, if v' ~ 10^ v, then the mirror neutrinos i^' with 
masses of keV order could provide the warm dark matter component in the Universe 
[16, 21, 24]. Instead, if (0') = 0, the u' are massless and unmixed with the ordinary 
neutrinos. For our considerations with exact mirror Parity, we have v' — v. 
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Let us discuss now the leptogenesis mechanism in our scenario. A crucial role in 
our considerations is played by the reheating temperature Tr, at which the inflaton 
decay and entropy production of the Universe is over, and after which the Universe 
is dominated by a relativistic plasma of ordinary particle species. As we discussed 
above, we assume that after the postinflationary reheating, different temperatures are 
estabhshed in the two sectors: < Tr, i.e. the mirror sector is cooler than the 
visible one, or ultimately, even completely "empty". This situation is motivated by 
the BBN constraints, and it can be achieved in the context of certain inflationary 
models [16, 20, 21, 24, 103]. In addition, the two particle systems should interact very 
weakly so that they do not come in thermal equilibrium with each other after reheating. 
The heavy neutrino masses are much larger than the reheating temperature Tr and 
thus cannot be thermally produced. As a result, the usual leptogenesis mechanism via 
N ^ Icf) decays is ineffective. 

Now, the important role is played by lepton number violating scatterings mediated 
by the heavy neutrinos N. The "cooler" mirror world starts being "slowly" occupied 
due to the entropy transfer from the ordinary sector through the AL — 1 reactions 
li(f) — > l'^4>': k4> ~^ I'k'P'- III general these processes violate CP due to complex Yukawa 
couphngs in eq. (4.21), and so the cross-sections with leptons and anti-leptons in the 
initial state are different from each other. As a result, leptons leak to the mirror sector 
more (or less) effectively than antileptons and a non-zero B — L is produced in the 
Universe. It is essential that these processes stay out of equilibrium. In other words, 
their rate should be less than the Hubble parameter H — 1.66 gl^^T'^/Mpi (g^ being 
the effective number of particle degrees of freedom) for temperatures T <Tr. 

For the rate of AL = 1 reactions we have Fi = aiUgq, where Ueq ~ (1.2/7r^)T^ is an 
equilibrium density per degree of freedom and ai is the total cross section of Z0 — > l'^)' 
scatterings 



The sum is taken over all flavor and isospin indices of initial and flnal states, and 
Qi = Tr{DW) = Tr[{h'^h')g-\h^h)*g-^]. Hence, the out-of-equilibrium condition for 
this process reads as 



However, there are also scattering processes with AL = 2 like l(p l(p etc., which 
can wash out the produced B — L unless they are out of equilibrium [81, 82]. Their total 



rate is given as F2 ~ {2,Q2/^TiM^)n^ where Q2 = Ti{A^A) = Ti[{h%)g-^{h%)*g-\ 



Qi 



(4.24) 



87rM2 ■ 




(4.25) 
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Figure 4.3: Trcc-lcvcl and one-loop diagrams contributing to the CP-asymmetries in Icf) 
l'4>' (left column) and Icf) — > I' (J)' (right column). 



Therefore, for a given reheat temperature Tr, the eq. (4.25) and the analogous con- 
dition K2 — {r2/2H)ji < 1 translate into the lower limit on the heavy neutrino mass 
scale M 

MuQi '^^ > A.2g:^I^Tl'\ MisQa > l^Ag^^l^Tl'^ . (4.26) 

where M12 = (M/IO^^ geV), Tg = (T^/IO^ GeV) and g., ^ 100 in the standard 
model. Clearly, if the Yukawa constants hia and /i-^ are of the same order, the out-of- 
equilibrium conditions for AL = 1 and AL = 2 processes are nearly equivalent to each 
other. 

Let us turn now to CP-violation. In AL = 1 processes the CP-odd lepton number 
asymmetry emerges from the interference between the tree-level and one-loop diagrams 
of fig. 4.3. The tree-level amplitude for the dominant channel l(p goes as 1/M and 

the radiative corrections as For the channel /0 instead, both tree-level 

and one-loop amplitudes go as As a result, the cross section CP asymmetries 

are the same for both 1(f) I'cf)' and l(j) channels (on the contrary, the diagrams 

with I'cj)' inside the loops, not shown in fig. 4.3, yield asymmetries, iAcr', symmetric 
to each other). However, CP- violation takes also place in AL — 2 processes (see fig. 
4.4). This is a consequence of the very existence of the mirror sector, namely the 
contribution of the mirror particles to the one-loop diagrams of fig. 4.4. The direct 
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Figure 4.4: Trcc-lcvcl and one-loop diagrams contributing to the CP-asymmetry of — > 1(f). 
The external leg labels identify the initial and final state particles. 

calculation gives 

a{l(l) l'4>') - am l'(f)') = (-A(7 - Aa')/2 , 

a{l(j) - am = {-Aa + Aa')/2 , 

a{l(t) U) - a(l4> 14)) = Aa ; (4.27) 

where J = lmTT[{h''^h')g^'^{h'^h)g^^mh)*g^^] is the CP- violation parameter and 5* is 
the cm. energy square (Act' is obtained from Act by exchanging h with h'). 

This is in perfect agreement with CPT invariance that requires that the total cross 
sections for particle and anti-particle scatterings arc equal to each other: (rij^cf) X) = 
a(l(j) — > X). Indeed, taking into account that a{l(f) l(p) = ailcp l(p) by CPT, we 
see that CP asymmetries in the AL — 1 and AL — 2 processes should be related as 

a{l(t) X') - am X') = -[a{l(f) 14>) - a^ 1(f))] = -Act , (4.29) 

where X' are the mirror sector final states, l'(f)' and I'ci)' . That is. the AL = 1 and 
AL — 2 reactions have CP asymmetries with equal intensities but opposite signs. 
But, as L varies in each case by a different amount, a net lepton number decrease is 
produced, or better, a net increase ol B — L oc Act (recall that the lepton number 
L is violated by the sphaleron processes, while B — L is changed solely by the above 
processes). 

As far as we assume that the mirror sector is cooler and thus depleted of particles, 
the only relevant reactions are the ones with ordinary particles in the initial state. 
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Hence, the evolution of the B — L number density is determined by the CP asymmetries 
shown in eqs. (4.27) and obeys the equation 



dUB-L . ^rj 3 2 



+ 3HnB-L = 7 A(7 . (4.30) 



Since the CP asymmetric cross section Aa is proportional to the thermal average cm. 
energy square S ~ 17T^ and H = 1/2 1 oc T^, one integrates the above equation from 
T = Tr to the low temperature limit and obtains the final B — L asymmetry of the 
Universe as 



B-L^^ 



Aan, 



AHs 



(4.31) 

R 



where s = {27i'^ /A5)g^,T^ is the entropy density. 

The following remark is in order. In fact, the lepton number production starts as 
soon as the inflaton starts decaying and the particle thermal bath is produced before 
the reheating temperature is established. (Recall that the maximal temperature at the 
reheating period is usually larger than Tr.) In this epoch the Universe is still dominated 
by the inflaton oscillations and therefore it expands as t^/'^ while the entropy of the 
Universe grows as t^^^. The integration of eq. (4.30) from some higher temperatures 
down to Tr gives an asymmetry 1.5 times larger than the estimation (4.31). Taking 
all these into account, the flnal result can be recasted as follows ^° 

where we have taken again 100. Taking also into account the lower hmits (4.26), 
we obtain the upper limit on the produced B — L 

B-L< 10-' ^ ; g = max{gi, 6 Q2} . (4.33) 

This shows that for Yukawa constants spread e.g. in the range 0.1 — 1, one can achieve 
B — L — 0(10"^^) for a reheating temperature as low as Tr ~ 10^ GeV. Interestingly, 
this coincidence with the upper bound from the thermal gravitino production, Tr < 
4 X 10^ GeV or so [59], indicates that our scenario could also work in the context of 
supersymmetric theories. 



Observe that the magnitude of the produced B — L strongly depends on the temperature, namely, 
larger B — L should be produced in the patches where the plasma is hotter. In the cosmological context, 
this would lead to a situation where, apart from the adiabatic density /temperature perturbations, 
there also emerge correlated isocurvature fluctuations with variable B and L which could be tested 
with the future data on the CMB anisotropics and large scale structure. 
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As far as the mirror sector includes the gauge couphng constants, which are the 
same as the standard ones, the mirror particles should be thermalized at a temperature 
T'. Once Ki < 1, T' will be smaller than the parallel temperature of the ordinary 
system T. Obviously, the presence of the out-of-equihbrium hidden sector does not 
affect much the Big Bang nucleosynthesis epoch. Indeed, if the two sectors do not come 
into full thermal equilibrium at temperatures T ~ Tr, then they evolve independently 
during the Universe expansion and approach the nucleosynthesis era with different 
temperatures. For Ki < 1, the energy density transferred to the mirror sector will be 
crudely p' {SKi/g^)p, where 100) is attained to the leptogenesis epoch. Thus, 
assuming that at the BBN epoch the mirror sector is dominated by relativistic degrees 
of freedom, we obtain an effective number of extra light neutrinos AA^,^ ^ Ki/2. 

Now it is important to stress that this mechanism would generate the baryon asym- 
metry not only in the observable sector, but also in the mirror sector. In fact, two 
sectors are completely similar, and have similar CP-violating properties. We have 
scattering processes which transform the ordinary particles into their mirror partners, 
and CP-violation effects in this scattering owing to the complex coupling constants. 
These exchange processes are active at some early epoch of the Universe, perhaps even 
enough close to equilibrium. In this hypothetical O observer could detect dur- 

ing the contact epoch that (i) matter slowly (in comparison to the Universe expansion 
rate) disappears from the thermal bath of our world, and, in addition, (ii) particles 
and antiparticles disappear with different rates, so that after the contact epoch ends 
up, he observes that his world is left with non-zero baryon number even if initially it 
was baryon symmetric. 

On the other hand, his mirror analogue, M observer, should see that (i) the mat- 
ter creation takes place in his world, and (ii) particles and antiparticles emerge with 
different rates. Therefore, after the contact epoch, he also would observe the non-zero 
baryon number in his world. 

One would naively expect that in this case the baryon asymmetries in the O and 
M sectors should be literally equal, given that the CP-violating factors are the same 
for both sectors. However, we show in that in reality, the BA in the M sector, since 
it is colder, can be about an order of magnitude bigger than in O sector, as far as 
washing out effects are taken into account. Indeed, this effects should be more efficient 
for the hotter O sector while they can be negligible for colder M sector, which could 
provide reasonable differences between two worlds in case the exchange process is 
not too far from equilibrium. The possible marriage between dark matter and the 
leptobaryogenesis mechanism is certainly an attractive feature of our scheme. 
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Now, let us discuss how the mechanism considered above produces not only the lep- 
ton number in the ordinary sector, but also the lepton prime asymmetry in the mirror 
sector. The amount of this asymmetry will depend on the CP- violation parameter that 
replaces J in eqs. (4.27) and Act', namely J' = ImTr[(h''h)g-'^(h"'h')g-\h'^h')*g-'^]. 
The M parity under the exchange (f) ^ I ^ f , etc., implies that the Yukawa cou- 
phngs are essentially the same in both sectors, = Therefore, in this case also 
the CP-violation parameters are the same, J' — J. 

Therefore, one naively expects that ub-l — "it^'b-l the mirror baryon number 
density should be equal to the ordinary baryon density, = f^b-^^ However, now we 
show that if the AL — 1 and AL — 2 processes are not very far from equilibrium, i.e. 
-f^i,2 ~ 1, the mirror baryon density should be bigger than the ordinary one. 

Here we should notice that eq. (4.30) for ub-l was vahd if these processes were 
very far from equilibrium, Ki^2 <^ 1. If instead Ki^2 ~ 1 then eq. (4.30) should be 
modified as 

+ SHuB-L + PriB-L = ^Aa nl^ , (4.34) 

where F = Fi Fa = (Qi ^Q2)n^/%'KM'^ is the total rate of the AL = 1 and AL = 2 
reactions. Now, solving this equation we obtain the expression for S — L as 

B-L^{B-L)o- D{2K) , (4.35) 

where {B — L)o is the solution of eq. (4.30), given by expressions (4.31) or (4.32), and 
the depletion factor D(k) is given by 

DkD{k)^le-''F{k) + lG{k) (4.36) 

where 

F{k) = f(2A; -If + 6k -5 + Ge'^'' 



G{k) = ^[2-{k' + 2k + 2)e-''\ . (4.37) 

These two terms in D[k) correspond to the integration of (4.34) respectively in the 

epochs before and after reheating (T > and T < Tr). Obviously, for A; <^ 1 

the depletion factor D{k) 1 and thus we recover the result as in (4.31) or (4.32): 

B — L = {B — L)o. However, for large k the depletion can be reasonable, e.g. for 

A; = 1, 2 we have respectively D{k) — 0.34, 0.1. 

^^Thc mirror parity could be also spontaneously broken by the difference in weak scales {<p) = v and 
((/>') = v' , which would lead to somewhat different particle physics in the mirror sector [16, 21, 23, 24], 
e.g. the mirror leptons and baryons could be heavier than the ordinary ones. But, as the mechanism 
only depends on the Yukawa constant pattern in (4.21), one still has ns-L = J^b-d while > Ob- 
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As far as the mirror sector is concerned, the evolution of the B — L number density, 
n'{B — L), obeys the equation 



^ + 3Hn'^_, + r'n^_L = | Aa < , (4.38) 



where now F' = T[ + T2 = {Qi + 6Q2)n'^q/87rM'^ is the total reaction rate of the 
AL' = 1 and AL' = 2 processes in the mirror sector, and n'^^ = (1.2/7r^)T'^ = x^Ueq is 
an equilibrium number density per a degree of freedom in a mirror sector. Therefore 
we have K' = r'/2H = x^K, and so for X ~ 1 and x < 0.64 (BBN hmit) K' <^ 1. 

For the mirror sector we have (B — L)' = — L)q- D{2K'), where the depletion can 
be irrelevant. Now taking into the account that in both sectors the B — L densities are 
reprocessed in the baryon number densities by the same sphaleron processes, we have 
B — a{B — L) and B' — a{B — L)', with coefficients a 1/3 equal for both sectors. 
Therefore, we see that the cosmological densities of the ordinary and mirror baryons 
should be related as 

- 'D(2Kj ■ ^^-^^^ 
Therefore, for K <C 1 and thus K' — x^K <^ 1, depletion factors in both sectors are 
D D' Ki 1 and thus we have that the mirror and ordinary baryons have the same 
densities, ~ Jib- In this case mirror baryons are not enough to explain all dark 
matter and one has to invoke also some other kind of dark matter, presumably cold 
dark matter. 

If instead X ~ 1, then we would have n(, > Qb, Sind thus all dark matter of the 
Universe could be in form of the mirror baryons. E.g., for K = 0.5 we have from 
eq. (4.39) that il'^ 5 fib, and hence for 0.05 we have ~ 0-25, which is 

exactly about the best fit value of the dark matter density. 



4.5 Primordial nucleosynthesis and mirror helium 
abundance 

The time scales relevant for standard BBN are defined by the "freeze-out" temperature 
of weak interactions Tw ~ 0.8 MeV {tw ~ 1 s) and by the "deuterium bottleneck" 
temperature Tjv ^ 0.07 MeV {In ~ 200 s) [104]. When T > Tw, weak interactions 
transform neutrons into protons and viceversa and keep them in chemical equilib- 
rium. The neutron abundance X„ = Un/uB, defined as the ratio of neutron to baryon 
densities, is given by Xn{T) = [1 -|- exp(Am/r)]~^, where Am ~ 1.29 MeV is the 
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neutron-proton mass difference. For T < Tw the weak reaction rate ~ GpT^ 
drops below the Hubble expansion rate H{T) ~ 5.5T'^/Mpi, the neutron abundance 
freezes out at the equilibrium value X„(T^) and it then evolves only due to the neutron 
decay: Xn{t) — exp(— i/r), where r = 886.7 s is the neutron hfetime. 

At temperatures T > Tn, the process p + n-^d + ^ is faster than the Universe 
expansion, and free nucleons and deuterium are in chemical equilibrium. The light 
element nucleosynthesis essentially begins when the system cools down to the temper- 
ature 

where Bd = 2.22 MeV is the deuterium binding energy, and uin is the nucleon mass. 
Thus, the primordial ^He mass fraction is 

n - 2XM = ^^"P(-;^/;) ^ ^ 0.24 . (4.41) 

As we have already discussed (see § 4.2), the presence of the mirror sector with a 
temperature T' <T has practically no impact the standard BBN, in the limit x < 0.64, 
which in fact has been set by uncertainties of the present observational situation. In 
the mirror sector nucleosynthesis proceeds along the same lines. However, the impact 
of the O world for the mirror BBN is dramatic! 

For any given temperature T', using eqs. (4.7) and (4.8) now we have 

H{T') ~ 5.5(1 + x~y^^T'^/Mpi (4.42) 

for the Hubble expansion rate. Therefore, comparing H{T') with the reaction rate 
r(r') oc T'^ (see eq. (1.55)), we find a freeze-out temperature = (1 + x'^^fl^Tw, 
which is larger than Tw, whereas the time scales as t'y^ — tw/{'^ + x~^)^/^ < tw 
(obtained using eq. (4.7) and the relation t oc H"^). In addition, 77' is different from 
77 ~ 5 X 10~^°. However, since Tn depends on baryon density only logarithmically (see 
eq. (4.40)), the temperature TJ^ remains essentially the same as TV, while the time t'j^ 
scales diS t'j^ — / {1 -\- x""^)^/^. Thus, for the mirror ^He mass fraction we obtain: 

V'^OY' (f \ 2cxp[-t;v/r(l+x-^)V^] 

- ^XnitN) - i + exp[Am/Tv.(l + a;-)V6] " ^^'^^^ 

We see that is an increasing function of x~^. In particular, for a; ^ one has 

In reality, eq. (4.43) is not valid for small x, since in this case deuterium production 
through reaction n -\- p <-> d -|- 7 can become ineffective. By a simple calculation 
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Figure 4.5: The primordial mirror ^Hc mass fraction as a function of x. The dashed curve 
represents the approximate result of eq. (4.41). The solid curves obtained via exact numerical 
calculation correspond, from bottom to top, to r]' varying from 10"^° to 10"^. 

one can make sure that for x < 0.3 • {rj' x 10^°)"^/^, the rate at which neutrons are 
captured to form the deuterium nuclei, T'^ = n'^aN ~ 'q'n'^aN, where n'^ ~ T'^ is 
the M photons density and — 4.5 ■ 10~^° cm^ s~^ is the thermal averaged cross 
section, becomes smaller than the Hubble rate H{T') for temperatures T' > T'j^. In 
this case M nucleosynthesis is inhibited, because the neutron capture processes become 
ineffective before deuterium abundance grows enough to initiate the synthesis of the 
heavier elements. Therefore, for any given Yl first increases with increasing l/x, 
reaches a maximum and then starts decreasing. The true dependence of on the x 
computed for different values of rj' by standard BBN code [98], is presented in fig. 4.5. 
The Hubble expansion rate of the mirror world was implemented, for each value of x, 
by taking an effective number of extra neutrinos as AiVj, = 6.14 ■ 

We have to remark, however, that in the most interesting situation when /5 — 
Vt'^,/Vt^, = x^rj'/rj > 1, the condition a; < 0.3 • [r]' x 10^°)"^/^ is never fulfilled and the 
behaviour of Y^ is well described by the approximate formula (4.43). Hence, in this 
case Yl is always bigger than I4. In other words, if dark matter of the Universe is 
represented by the baryons of the mirror sector, it should contain considerably bigger 
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fraction of primordial ^He than the ordinary world. In particular, the helium fraction 
of mirror matter is comprised between 20% and 80%, depending on the values of x 
and r]'. This is a very interesting feature, because it means that mirror sector can 
be a helium dominated world, with important consequences on star formation and 
evolution, and other related astrophysical aspects, as we will see in chapter 7. 



4.6 Mirror baryons as dark matter: general discus- 



As we explained in § 4.1, mirror matter is an ideal candidate for the inferred dark matter 
in the Universe, being naturally stable and "dark" , interacting only gravitationally with 
the ordinary matter. 

At present time various observations indicate that the amount of dark matter is 
order 10 times the amount of ordinary matter in the Universe (for a picture of the 
present observational status see § 1.7). This is not a problem for the mirror matter 
theory, because a mirror symmetric microscopic theory does not actually imply equal 
numbers of ordinary and mirror atoms in the Universe, as someone could erroneously 
think. The point is that, according to what we said in the previous sections, the initial 
conditions need not be mirror symmetric, and the Universe could have been created 
with more mirror matter than ordinary matter. 

In the most general context, the present energy density contains relativistic (radia- 
tion) component $1^, non-relativistic (matter) component and the vacuum energy 
density (cosmological term). According to the inflationary paradigm the Universe 
should be almost flat, Qq — + 0,^ + 0,^. ~ 1, which well agrees with the recent 
results on the CMB anisotropy [9, 109]. For redshifts of the cosmological relevance, 
1 + z — T/Tq 1, the Hubble parameter is expressed by (using eqs. in appendix A. 2) 



where now Qr and flm represent the total amount of radiation and matter of both 
ordinary and mirror sectors. In the context of our model, the relativistic fraction is 
represented by the ordinary and mirror photons and neutrinos, and, using eq. (4.8) 
and the value of the observable radiation energy density flrh^ ~ 4.2x10"^, it is given 



sion 



H{z) = Ho Qril + zf + Q^(l + z) 



3^1/2 



(4.44) 



by 



= 4.2 X 10"^ (1 + x^) ~ 4.2 X 10"^ , 



(4.45) 



§4.6 Mirror baryons as dark matter: general discussion 



75 



where the contribution of the mirror species is neghgible in view of the BBN constraint 
X < 0.64. As for the non-relativistic component, it contains the O baryon fraction Q;, 
and the M baryon fraction Q'j, = /^Q^, while the other types of dark matter could also 
present. Obviously, since mirror parity doubles all the ordinary particles, even if they 
are "dark" (i.e., we are not able to detect them now), whatever the form of dark matter 
made by some exotic ordinary particles, there will exist a mirror partner made by the 
mirror counterpart of these particles. In the context of supersymmetry, the CDM 
component could exist in the form of the lightest supersymmetric particle (LSP). It 
is interesting to remark that the mass fractions of the ordinary and mirror LSP are 
related as JIlsp — ^^lsp- Iii addition, a significant HDM component could be due 
to neutrinos with order eV mass. The contribution of the mirror massive neutrinos 
scales as Q'^ — x^Qi, and thus it is irrelevant. In any case, considering the only CDM 
component, which is now the preferred candidate, we can combine both the ordinary 
and mirror components, given that their physical effects are exactly the same. So, in 
our work we consider a matter composition of the Universe made in general by 

= Ob + + QcDM ■ (4.46) 

The important moments for the structure formation are related to the matter- 
radiation equality (MRE) epoch and to the plasma recombination and matter-radiation 
decoupling (MRD) epochs (see § 1.5). The MRE occurs at the redshift 

l + .eq = ^-2.4.10^i^, (4.47) 

which is always smaller than the value obtained for an ordinary Universe, but approx- 
imates it for low X (see fig. 4.6). If we consider only ordinary and mirror baryons and 
photons, we find 

^ P^+Pk ^ P^i^ + (^) ii^,)^ P\ + (^'l (1 + ,) , (4.48) 

^ P^o pto + p^o P7(i + ^^) p;(i + ^ > 

where pb and p-y indicate respectively the sums of baryons and photons of the two 
sectors. So, in the presence of a mirror sector the matter - radiation equality epoch 
shifts as 

The MRD, instead, takes place in every sector only after the most of electrons 
and protons recombine into neutral hydrogen and the free electron number density 
diminishes, so that the photon scattering rate = rieCr = Xerjn-yaT drops below 
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the Hubble expansion rate H{T), where ax = 6.65 ■ 10~^^ cm^ is the Thomson cross 
section. In condition of chemical equilibrium, the fractional ionization Xg = rie/nb is 
given by the Saha equation, which for <^ 1 reads 

- (1 - Y,f^ ^ f^y'^\^^/^^ , (4.50) 

where B = 13.6 eV is the hydrogen binding energy. Thus we obtain the familiar result 
that in our Universe the MRD takes place in the matter domination period, at the 
temperature Tdcc — 0.26 eV which corresponds to redshift 1 + Zdec = T^ec/To ~ 1100. 



10^ 
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Fi gure 4.6: The M photon decoupling redshift 1 + z'^^^ as a function of x (solid line). 
The long-dash line marks the ordinary decoupling redshift 1 + z^ec = 1100. We also show 
the matter-radiation equality redshift 1 + Zeq for the cases n^h'^ = 0.2 and Vt^h^ = 0.6 
(respectively lower and upper dash). 

The MRD temperature in the M sector T'^^^ can be calculated following the same 
lines as in the ordinary one. Due to the fact that in either case the photon decoupling 
occurs when the exponential factor in eq. (4.50) becomes very small, we have T'^^^ ~ 
Tdec, up to small corrections related to r^', Yl different from 77, Ki. Hence, considering 
that T' = X ■ T (see § 4.2), we obtain 

1 + 4,, ~ x-\l + Zdec) ^ 1.1 ■ lO'x^^ , (4.51) 

so that the MRD in the M sector occurs earher than in the ordinary one. Moreover, 
comparing eqs. (4.47) and (4.51), which have different trends with x, we find that, for 
X less than a value Xeq given by 

Xeq « 0.046(fi„/i2)-i , (4.52) 



§4.6 Mirror baryons as dark matter: general discussion 



77 



the mirror photons would decouple yet during the radiation dominated period (see 
fig. 4.6). Assuming, e.g., the "standard" values Qrn — 0.3 and h — 0.65, we obtain 
the approximate value Xeq — 0.36, which indicates that below about this value the 
mirror decoupling happens in the radiation dominated period, whith consequences on 
structure formation (see next chapter). 

We have shown that mirror baryons could provide a significant contribution to the 
energy density of the Universe and thus they could constitute a relevant component 
of dark matter. Immediate question arises: how the mirror baryonic dark matter 
(MBDM) behaves and what are the differences from the more familiar dark matter 
candidates as the cold dark matter (CDM), the hot dark matter (HDM), etc. In next 
chapters we discuss the problem of the cosmological structure formation in the presence 
of M baryons as a dark matter component. Namely, M baryons being a sort of self- 
interacting dark matter could provide interesting signatures on the CMB anisotropy, 
the large scale structure of the Universe, the form of the galactic halos, microlensing, 
etc. 



Chapter 4. The Mirror Universe 



Chapter 5 



Structure formation for a Mirror 
Universe 



5.1 Introduction 

In chapter 2 wc presented the hnear structure formation theory in a non-baryonic dark 
matter scenario (which is now considered the "standard" paradigm). We introduced 
some fundamental concepts (hke the Jeans length and the coUisionless dissipation) 
and gave both the Newtonian and the general relativistic approaches (briefly reviewed 
in appendix B). Here, we extend the theory to the case of dark matter with a non- 
negligible mirror baryon component. 

In a Mirror Universe we assume that a mirror sector is present , so that the matter is 
made of ordinary baryons (the only certain component), non-baryonic (dark) matter, 
and mirror baryons. Thus, it is necessary to study the structure formation in all 
these three components. This requires essentially the study of the baryonic structure 
formation (given also that we already described the non-baryonic case in chapter 2). 
There are two reasons for that: first, the baryons are the only certain component 
of the universal matter (even if the lowest); second, as seen in the previous chapter, 
the physics of mirror baryons is the same as ordinary ones but with different initial 
conditions and time-shifted key epochs. 

As for the non-baryonic dark matter, we continue to focus only on adiabatic per- 
turbations, that we introduced in § 2.2. Here we remember only that an adiabatic 
perturbation satisfies the condition for adiabaticity 

3 

(^m=|(^r, (5.1) 
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which relates perturbations in matter (^m) and radiation (Sj.) components. 



5.2 Baryonic structure formation (ordinary and 
mirror) 



We will now consider cosmological models where baryons, ordinary or mirror, are the 
dominant form of matter, but it should be made clear at the outset that purely ordi- 
nary baryonic models cannot successfully explain the origin of the observed structure. 
Nevertheless, it's very important to look at the details of these scenarios for two rea- 
sons. First, even if they are very few, ordinary baryons do exist in the Universe, and 
it's necessary to know as they behave. Second, mirror baryons, which follow the same 
physics of the ordinary ones, could be the dominant form of matter in the Universe. 
Thus, it is crucial to study the interaction between baryonic matter and radiation dur- 
ing the plasma epoch in both sectors, and the simplest way of doing it is by looking at 
models containing only these two matter components. 

As explained in appendix B, the relevant length scale for the gravitational insta- 
bilities is characterized by the Jeans scale (see § B.l), which now needs to be defined 
in both the ordinary and mirror sectors. 

5.2.1 Evolution of the adiabatic sound speed 

Given the expression of the Jeans length (B.16), it is clear that the key issue is the 
evolution of the sound speed, since it determines the scale of gravitational instability. 
If we remember eq. (B.ll), we obtain for the two sectors 



where S denotes the total entropy and w' is relative to the mirror equation of state 
p' — w'p'. 

The ordinciry sound speed 

In a mixture of radiation and baryonic matter the total density and pressure are p = 
P-y + Pb and p pr^ — p-y/3 respectively (recall that pb — 0). Hence, the adiabatic 





(5.2) 
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sound speed is given by 

^^V^4fl + ^V". (5.3) 



' \dpj V3\ 4pJ 
where we have used the adiabatic condition (5.1). In particular, using eqs. (1.22) and 
(1.23) together with the definition of matter - radiation equahty (1.59) (where now we 
consider that matter is only baryonic), we obtain 



1 + 



3 /l + ^eqM-^/' 



(5.4) 



4 V 1 + ^ 

In fact, the relation above is valid only in an ordinary Universe, and it is an approxima- 
tion, for small values of x and the mirror baryon density (remember that (3 — Q^flb), of 
the more general equation for a Universe made of two sectors of baryons and photons, 
obtained using eqs. (5.3) and (4.48) and given by 



^ Q /I I ™4\ /I I ^ \ ^1/2 

v.(z) ^ ^ 




x'^\ /I + z, 



eq 



(5.5) 



1 + P J \l + z 

In the most general case, the matter is made not only of ordinary and mirror baryons, 
but also of some other form of dark matter, so the factor is replaced by I+^+^dm, 
where jSoM = (^m — — ^b) l^b- From eqs. (5.4) and (5.5) wc obtain that, given the 
conditions x < 1 (from the BBN bound presented in § 4.2) and f5 > 1 (cosmologically 
interesting situation, i.e., significant mirror baryonic contribution to the dark matter), 
the ordinary sound speed in a Universe made of two sectors is always higher than that 
in an ordinary Universe. This is due to the presence of the term [(1 + x*)/{l + /?)] 
linked to the shift of matter-radiation equality epoch (4.49). 

Now we define the scale factor corresponding to the redshift 

(1 + z^yy) = (ab^)-^ = (Oft/O^) = 3.9 • 10^(06/1^) . (5.6) 

Given that 1 -|- Zj-ec — 1100, baryon-photon equipartition occurs before recombination 
only if Qb^^ > 0.026 (which seems unlikely, given the current estimates reported in 
§ 1.7). On the other hand, a^^ is always higher than Og^, given that Jim > f^b- In the 
radiation era p^,, ensuring that Vg — l/v^. In the interval between equipartition 
and decoupling, when pb ^ P-y, eq. (5.3) gives Vs — ^^Ap^/Spi, oc a~^/^. After decou- 
pling there is no more pressure equilibrium between baryons and photons, and Vs is 
just the velocity dispersion of a gas of hydrogen and helium, Vs oc a"^. The situation 
whit Qb^^ > 0.026 is resumed below: 

const a < a}yy , 
Vs{a) oc \ a~^/^ aiyy <a < a^ec , (5.7) 

a > Qdec ■ 
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If Qy^h"^ < 0.026, which is a more probable and simple case, at-y > ddec-, and the 
intermediate situation does not arise. It's very important to observe (see Padmanabhan 
(1993) [132]) that at decoupling drops from (p-y/pb) to (pb/Pb)- Since p^ oc n^T while 
Pb oc UbT with {n^/ub) — 10^ ^ 1, this is a large drop in t>s and consequently in Aj. 
More precisely, v"^ drops from the value {l/3){py/ pb) = (l/3)(i7^/i7f,)(l + Zdec) to the 
value {5/3){Tdec/mb) = (5/3)(To/mfe)(l + Zdec), with a reduction factor 



Fi(fib/i^ > 0.026) 



(^S )dGC 



6.63 ■ lO-^iQbh^) 



(5.8) 



where (fs)dcc [v'sY^^ indicate the sound speed respectively just before and after 
decoupling. In the case fib/i^ < 0.026, drops directly from (1/3) to {5/3){Tdec/mb) = 
(5/3)(To/mb)(l + Zdec) with a suppression 



F2(fib/i^ < 0.026) 



1.9 ■ 10"^ . 



(5.9) 



In fig. 5.1 we plot the previously discussed trends of the ordinary sound speed for a 
typical model with fib^^ > 0.026. If we reduce the value fib^^, cih-y goes toward higher 
values, while Odec remains fixed, so that for Q\,h'^ < 0.026 decouphng happens before 
equipartition and the intermediate regime, where oc a^^/^, disappears. 



1/3' 



const. 



a 



Figure 5.1: The trends of the ordinary sound speed as discussed in the text for the case 
J^b/i^ = 0.08; in this case there is also the intermediate regime Vs oc a ' . 
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The mirror sound speed 

The mirror plasma contains more baryons and less photons than the ordinary one, 



= Pph and p'^ = X py. Then, using eqs. (5.3) and (4.48), we have 



V's(z) 




-1/2 ^ 




l+x-'^\ fl + z. 



cq 



.-1/2 



(5.10) 



Let us consider for simplicity the case when dark matter of the Universe is entirely 
due to M baryons, ~ Q'^ (i.e., /3 » 1). Hence, for the redshifts of cosmological 
relevance, z ~ ^^eq, we have v'^ ~ 2x^/3, which is always less than Vg ~ l/-\/3 (some 
example: if x = 0.7, Vg ^ 0.5 • Vg] ii x — 0.3, Vg ^ 0.1 • Vg). In expression (5.10) 
it is crucial the presence of the factor [(1 + x~^)/{l + which is always greater 

than 1, so that v'g < Vg during all the history of the Universe, and only in the limit 
a <S Oeq we obtain ~ ~ 1/^3. As we will see in the following, this has important 
consequences on structure formation scales. 

According to what found in § 4.6, in the mirror sector the scale of baryon - photon 
equality a'^^ is dependent on x; if we remember the definition of the quantity Xeq ~ 
0.046(fi„i/i^)~^, we find that for x > x^q the decoupling occurs after the equipartition 
(as in the ordinary sector for fl]^h'^ > 0.026), but for x < x^q it occurs before (as for 
Qb^^ < 0.026). It follows that, by taking care to interchange a\yy with and a^ec 
with a^gp, we have for the sound speed the same trends with the scale factor in both 
sectors, as expressed in eq. (5.7), though with the aforementioned differences in the 
values. 

In fact, the matter-radiation equality for a single sector (ordinary) Universe, 
(aeq)ord) IS always bigger than that for a two sectors (mirror) one, (aeq)mir, accord- 
ing to 

_ (1 + x^) 

(Qeq)mir — _|_ (Qeq)ord < (Qeq)ord ) (^■l-'-) 

while the baryon-photon equipartition transforms as 

a^. = ^ - ^ = a., J < a,, . (5.12) 
Recalling our hypothesis x < 1 and /5 > 1, it is always verified that 

a'b^ < Oeq < . (5.13) 

If we consider now the drop in (Wg)^ at decoupling and call F[, F2 the factors of 
this drop in the cases when, respectively, a'^^ < a'^^^ or > a'^^^, we find 

F[ = /3x-^Fi and F^ ^ F2 . (5.14) 
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Some example: for x = 0.7, Fl ^ 2.9f3Fi; for x = 0.5, F[ = 8f3Fi; for x = 0.3, 
F[ ^ 37/9 -Fl. We note that, if /? > 1, F[ is at least an order of magnitude larger than 
Fl. In fact, after decoupling {vi,)^ = (5/3)(r^^^/mb) = (5/3)(Tdec/mb) = {v^y (given 
T^gj, = Tdcc), and between equipartition and recombination (t>g)^ < (fs)^. The relation 
above means that the drop is smaller in the mirror sector than in the ordinary one. 
Obviously, before equipartition (fg)^ = (fg)^ = 1/3, and for this reason the parameter 
F2 is the same in both sectors. 

In figure 5.2 we show the trends with scale factor of the mirror sound speed, in 
comparison with the ordinary one. The ordinary model is the same as in figure 5.1, 
while the mirror model has x = 0.6 and (3 = 2 (this means that mirror baryons are 
twice the ordinary ones, but in these models we chose the latter four times their current 
estimation to better show the general behaviour). In the same figure are also indicated 
the aforementioned relative positions of the key epochs (photon-baryon equipartition 
and decoupling) for both sectors, together with the matter-radiation equahty. 




Figure 5.2: The trends of the mirror sound speed (blue line) as a function of the scale factor, 
compared with the ordinary sound speed (red line). The ordinary model has Q.\yh? = 0.08 
(as in figure 5.1), while the mirror model has x = 0.6 and (3 = 2. Are also reported all the 
key epochs: photon-baryon equipartition and decoupling in both sectors, and the matter- 
radiation equality. 

In figure 5.3 the same ordinary and mirror sound speeds are plotted together with 
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the velocity dispersion of a typical non baryonic cold dark matter candidate of mass 
~ IGeV. Note that the horizontal scale is expanded by some decade, because the key 
epochs for the CDM velocity evolution (when the particles become non relativistic, 
ttnr, and decoupling, Od) occur at a much lower scale factor (for a brief review of the 
dark matter see chapter 2). 




Figure 5.3: The trends of the mirror (blue line) and ordinary (red) sound speed compared 
with the velocity dispersion of a typical non baryonic cold dark matter candidate of mass 
~ 1 GeV (green) ; Onr and indicate the scale factors at which the dark matter particles 
become non relativistic or decouple. The models are the same as in figure 5.2, but the 
horizontal scale is expanded by some decade to show the CDM velocity. 



5.2.2 Evolution of the Jeans length and the Jeans mass 
Ordinary sector 

Recalling what stated in § B.l, we have the Jeans length 



X, = VsJ-;^, (5.15) 

V tjPdom 

where pdom is the density of the dominant species, and the Jeans mass 

Mj = ^7^pJ^y = ^pU\jf , (5.16) 
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where the density is now that of the perturbed component (baryons). 

Using the results of § 5.2.1 relative to the sound speed, we find for the evolution of 
the adiabatic Jeans length and mass in the case Q^^h'^ > 0.026 



Aj oc 



a 
a 



Mj oc ^ oc 



a 
a 



a < a, 



3/2 



'eq 



a. 



a 



1/2 



const. 



(5.17) 



In figures 5.4 and 5.5 we plot with the same horizontal scale the trends of the 
ordinary Jeans length and mass for Qi^h"^ = 0.08, a case in which the intermediate 
regime Ob-y < a < Odec is present. 




Figure 5.4: The trends of the ordinary Jeans length as discussed in the text for ^l^^h"^ = 0.08, 
a case in which the intermediate regime a^-y < a < Odec is present. The horizontal scale is 
the same as in figure 5.1. 



The greatest value of the Jeans mass is just before decoupling (see Padmanabhan 
(1993) [132]), in the interval a^^ < a < a^cc , where 



Mj(a < Ode 



1.47 ■ lO^^Mo ( ) (n^h'' 



-2 



that for = — O.l/i is ~ 10 Mq. Just after decoupling we have 



Mj(a> adcc) = 2.5-1O3M0 



(5.1^ 



(5.19) 
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that for r^jn = f2b — O.l/i ^is~ IO^Mq. This drop is very sudden and large, changing 
the Jeans mass by F^^^ ~ 1.7 ■ IQ-^^Q^h'^Y/'^ . 

The maximum possible value of the Jeans mass is obtained for Ob-y = a^ec, keeping 
constant aeq (i.e., substituting a fraction of baryons with the same amount of dark 
matter), and is given by 




(5.20) 



Figure 5.5: The trends of the ordinary Jeans mass as discussed in the text for 17^^^ = 0.08, 
a case in which the intermediate regime Ob-y < a < a^ec is present. The horizontal scale is 
the same used for plotting A j in figure 5.4. 

Otherwise, if flhh"^ < 0.026, ab7 > adec, there is no intermediate phase a^'y < a < 
Odec, and Mj(a < Odec) is larger 

Mj(a < adec) ^ 3.1 ■ lO^^M^ (l^m/^')"'^' , (5.21) 

while after decoupling it takes the value in eq. (5.19), so that the drop is larger, 
F^^^ ~ 8.3 ■ 10-^\ 

We note that, with the assumption = fib = 1, ^j,max (which is the first scale 
to become gravitationally unstable and collapse soon after decoupling) has the size of 
a supercluster of galaxies. 
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Mirror sector 

In the mirror sector tlie Jeans lengtli and mass are 




(5.22) 



In this case it's no more sufficient to interchange a^-y with a^^ and a^ec with a'^^^, as 
made for the sound speed, because from relation (5.13) we note that in the mirror 
sector the photon-baryon equipartition happens before the matter-radiation equality 
(due to the fact that we are considering a mirror sector with more baryons and less 
photons than the ordinary one). It follows that, due to the shifts of the key epochs, 
the intervals of scale factor for the various trends are different. As usual, there are two 
different possibilities, x > x^q and x < x^q (which correspond roughly to Q]^h'^ > 0.026 
and f2b^^ < 0.026 in an ordinary Universe), where, as discussed in § 5.2.1, for the 
second one the intermediate situation is absent. 




Figure 5.6: The trends of the mirror Jeans length (blue line) as a function of the scale factor, 
compared with the ordinary Jeans length (red line). The ordinary model has f^b^^ = 0.08 
(as in figure 5.4), while the mirror model has x = 0.6 and (3 = 2. The horizontal scale is the 
same as in figure 5.2. We note that the same behaviours reported for the ordinary sector in 
figure 5.4 are present in the mirror sector for different intervals of scale factor. 
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Using the results of § 5.2.1 for the sound speed, we find the evolution of the adiabatic 
Jeans length and mass in the case x > Xeq 

M> oc , < " < > (5.23) 

a a-^ \ const. Ogq < a < a^^^ 

We plot in figures 5.6 and 5.8 the trends of the mirror Jeans length and mass com- 
pared with those for the ordinary sector; the parameters of both mirror and ordinary 
models are the ones previously used, i.e. flhh"^ = 0.08, x = 0.6 > Xeq and /? = 2. If we 
remember eq. (5.14), for the mirror model the drop in the Jeans mass at decoupling 
is (Fj')'^/^ = /^^/^^"^/^(Fi)'^/^, which, given our bounds on x and /3, is greater than 
(Fi)^/^. We give here some numerical example: for x = 0.7, (F[)'^^'^ ^ 5/?^''^(Fi)^''^; for 
X = 0.5, (F{)3/2 ^ 23/?3/2(ir^)3/2. for ^ ^ 0.3, (F;)3/2 ^ 225/?3/2(ir^)3/2. for ^ ^ 0.6 and 

p = 2 (the case of figures 5.6 and 5.8), {F{f/^ ^ 28(Fi)3/2. 




Figure 5.7: The trends of the mirror (blue line) and ordinary (red) Jeans length compared 
with that of a typical non baryonic cold dark matter candidate of mass ~ 1 GeV (green) ; Onr 
and Od indicate the scale factors at which the dark matter particles become non relativistic 
or decouple, respectively. The models are the same as in figure 5.6, but the horizontal scale 
is expanded by some decade to show the CDM Jeans length, as in figure 5.3. 
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Figure 5.8: The trends of the mirror Jeans mass (blue line) as a function of the scale factor, 
compared with the ordinary Jeans mass (red line). The models and the horizontal scale are 
the same as in figure 5.6. 

If we now consider the expression (5.12), we have 



which can be used to express the value of the mirror Jeans mass in the interval Ogq < 
a < Ojgj. (where Mj takes the maximum value) in terms of the ordinary Jeans mass in 
the corresponding ordinary interval a^^-y < a < a^ec- We obtain 



which, for (3 > 1 and x < 1, means that the Jeans mass for the M baryons is lower 
than for the O ones over almost the entire (/3-x) parameter space, with implications 
for the structure formation process. If, e.g., x = 0.6 and (3 = 2, then Mj ~ 0.03 Mj. 
We can also express the same quantity in terms of Qh, x and (3 only, in the case that 
all the dark matter is in the form of mirror baryons, as 
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It's important to stress that these quantities are strongly dependent on the values 
of the free parameters x and (3, which shift the key epochs and change their relative 
position. We can describe some case useful to understand the general behaviour, but 
if we want an accurate solution of a particular model, we must unambiguously identify 
the different regimes and solve in detail the appropriate equations. 




Figure 5.9: The trends of the mirror (blue line) and ordinary (red) Jeans mass compared 
with those of a typical non baryonic cold dark matter candidate of mass ~ 1 GeV (green) ; Onr 
and ad indicate the scale factors at which the dark matter particles become non relativistic 
or decouple, respectively. The model parameters and the horizontal scale are the same as in 
figure 5.7. 

In figures 5.7 and 5.9 we plot the trends of the mirror and ordinary Jeans length 
and mass compared with those of a typical non baryonic cold dark matter candidate of 
mass ~ IGeV. Apart from the usual expansion of the horizontal scale, due to the much 
lower values of the CDM key epochs as compared to the baryonic ones, a comparison 
of the mirror scenario with the cold dark matter one shows that the maximal value 
of the CDM Jeans mass is 10^^ times lower than that for mirror baryons, which is a 
very big value. This implies that a very large range of mass scales, which in a mirror 
baryonic scenario oscillate before decoupling, in a cold dark matter scenario would 
grow unperturbed during all the time (for more details see § 5.4). 
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For the case x < Xeq, both a[^^ and a'^^^ are smaller than the previous case x > Xeq, 
while the matter-radiation equality remains practically the same; as explained in § 4.6, 
the mirror decoupling (with the related drop in the associated quantities) happens 
before the matter-radiation equality, and the trends of the mirror Jeans length and 
mass are the following 





a < a^^ 






const. 






acq < a 



xWi""'^' Mjoc^oc<;;:_^^ (5-27) 

In figure 5.10 we plot the mirror Jeans mass for the three different possibilities: 
X = Xeq (the transition between the two regimes). 




Figure 5.10: The trends of the mirror Jeans mass for the cases x < x^q (sohd hne), x = Xgq 
(dotted) and x > Xcq (dashed). The model with x > Xeq is the same as in figure 5.8, the others 
are obtained changing only the value of x and keeping constant all the other parameters. 
As clearly shown in the figure, the only key epoch which remains almost constant in the 
three models is the matter-radiation equality; the mirror baryon-photon equipartition and 
decoupling indicated are relative to the model with x < x^q. It's also evident the change in 
the trends when x becomes lower than Xeq, due to the fact that a'^^^ becomes lower than Ceq. 
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5.2.3 Evolution of the Hubble mass 

The trends of the Hubble length and mass for the baryonic component are the same 
as for for the non-baryonic ones (defined in § 2.4.3) 



A„oc|«3/. MHOC^ocj;,/. (5.28) 

In this case, it should be emphasized that during the period of domination of 
photons (a < a^^) the baryonic Jeans mass is of the same order of the Hubble mass. 
In fact, following definitions (2.4), (2.5) and (2.12), given = (87r/3)Gp^ ~ 
we find 



2 (\„\3 ( ^3 



(Ah) j a < tte 



3 



\ ^ I - 26 . (5.29) 

We plot the trends of the Hubble mass in figures 5.12 and 5.13 together with other 
fundamental mass scales. 

5.2.4 Dissipative effects: coUisional damping 

We now turn our attention to physical dissipative processes that can modify the 
purely gravitational evolution of perturbations. In baryonic models the most im- 
portant physical phenomenon is the interaction between baryons and photons in the 
pre-recombination era, and the consequent dissipation due to viscosity and heat con- 
duction. 

Ordinciry baryons 

Adiabatic perturbations in the photon-baryon plasma suffer from coUisional damp- 
ing around the time of recombination because the perfect fluid approximation breaks 
down. As we approach decoupling, the photon mean free path increases and photons 
can diffuse from the overdense into the underdense regions, thereby smoothing out any 
inhomogeneities in the primordial plasma. The effect is known as "coUisional dissipa- 
tion" or "Silk damping" (Silk (1967) [164]). To obtain an estimate of the effect, we 
consider the physical (proper) distance associated with the photon mean free path 

A^ = — ^ ~ lO'^^a^X-^ (^bh'^Y^ cm , (5.30) 

where is the electron ionization factor, ne oc is the number density of the 
free electrons and ctt is the cross section for Thomson scattering. Clearly, all bary- 
onic perturbations with wavelengths smaller than A-y will be smoothed out by photon 
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free streaming. The perfect fluid assumption breaks down completely when A <ti A-y. 
Damping, however, occurs on scales much larger than A-y as the photons slowly diffuse 
from the overdense into the underdense regions, dragging along the still tightly coupled 
baryons. Integrating up to recombination time we obtain the total physical distance 
traveled by a typical photon (see Kolb & Turner (1990) [104]) 

As = ^^(A-y)rec^rec - 3.5 (^fe/l')"'''^ Mpc , (5.31) 

and the associated mass scale, which is known as the "Silk mass" , is given by 

Ms^^TTpbj^yj ~ 6.2 X 10^^ (n^/iy Mo , (5.32) 

which, assuming fl^h^ ~ 0.02 (as actually estimated, see § 1.7), gives Ms ~ 8x 10^^ Mq. 
The dissipative process we considered above causes that fluctuations on scales below 
the Silk mass are completely obliterated by the time of recombination and no structure 
can form on these scales. Alternatively, one might say that adiabatic perturbations 
have very little power left on small scales. 

Mirror bctryons 

In the mirror sector too, obviously, the photon diffusion from the overdense to under- 
dense regions induces a dragging of charged particles and washes out the perturbations 
at scales smaller than the mirror Silk scale 

A5 ~ 3 X /(x)(n^/i^)-^/^ Mpc , (5.33) 

where f{x) = x^l'^ for x > x^q and f{x) = {x/xeq)^^^xl^^ for x < x^q, and we considered 
the initial hypothesis Qm — ^i- 

Thus, the density perturbation scales which can run the linear growth after the 
matter-radiation equality epoch are limited by the length A^. The smallest perturba- 
tions that survive the Silk damping will have the mass 

M's ~ [f{x)/2]%n,^h')-'^H0'' Mo , (5.34) 

which should be less than 2 x 10^^ M© in view of the BBN bound x < 0.64. Interestingly, 
for X ~ Xea we obtain 



M'six = Xeq) ~ 10^ (Om/i')-' M, 



(5.35) 
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which, for the current estimate of flmh'^ (see § 1.7), gives M'g ~ 3 x 10^° Mq, a typical 
galaxy mass. 

At this point it is very interesting a comparison between different damping scales, 
collisional (ordinary and mirror baryons) and collisionless (non-baryonic dark matter) . 
Recalling what stated in § 2.4.2, we have that for hot dark matter (as a neutrino of 
mass ~10 eV) Mp^ ~ 10^^ Mq, while for a typical warm dark matter candidate of mass 
~1 keV, Mfs^^'^ ~ 10^ - 10^° Mq. From eq. (5.35) it is evident that the dissipative 
scale for mirror Silk damping is analogous to that for WDM free streaming. Conse- 
quently, the cutoff effects on the corresponding large scale structure power spectra are 
similar, though with important differences due to the presence of oscillatory features 
in the mirror baryons spectra, which makes them distinguishable one from the other 
(for a detailed presentation of the mirror power spectra see next chapter). In figure 
5.11 we show this comparison together with the trend of the mirror Silk mass over a 
cosmologically interesting range of x. 




0.1 X.. 1 X 



Figure 5.11: The trend of the mirror Silk mass (blue line) over a cosmologically interesting 
range of x, which contains Xeq (we considered Q^h? ~ 0.15, so Xeq — 0.3). The axis are both 
logarithmic. We show for comparison also the values of the ordinary Silk mass (red) and of 
the free streaming mass (green) for typical HDM and WDM candidates. 
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5.3 Scenarios 

After the description of the fundamental scales for structure formation, let us now put 
together all the informations and discuss the mirror scenarios. They are essentially 
two, according to the value of x. which can be higher or lower than x^q, and are shown 
respectively in figures 5.12 and 5.13, which will be our references during the present 
section. 

Typically, adiabatic perturbations with sizes larger than the maximum value of 
the Jeans mass, which is Mj(aeq) for x > Xeq and Mj{a'^^^) for x < Xgq, experience 
uninterrupted growth. In particular, as discussed in appendix B and summarized in 
table 2.1, they grow as Si, oc before matter-radiation equality and as 5b oc a after 
equality. Fluctuations on scales in the mass interval Mg < M < Mj max grow as 
5b oc while they are still outside the Hubble radius. After entering the horizon and 
until recombination these modes oscillate like acoustic waves. The amplitude of the 
oscillation is constant before equilibrium but decreases as a~^/^ between equipartition 
and recombination. After decoupling the modes become unstable again and grow as 
5b oc a. Finally all perturbations on scales smaller than the value of the Silk mass are 
dissipated by photon diffusion. 

Given this general behaviour, the schematic evolution of an adiabatic mode with 
a reference mass scale Mpdt, with Af^ < Mpgrt < Mj(aeq), is depicted in figure 5.12 
for X > Xeq- We distinguish between three evolutionary stages, called A, B and C, 
depending on the size of the perturbation and on the cosmological parameters flbh'^, 
X and j3, which determine the behaviour of the mass scales, and in particular the 
key moments (time of horizon crossing and decoupling) and the dissipative Silk scale. 
During stage A, i.e. before the horizon crossing (a < acnt < (icq), the mode grows as 
5b oc a^; throughout stage B (aent < a < a'^^c) ^^e perturbation enters the horizon, 
baryons and photons feel each other, and it oscillates; finally, in stage C (a > a^g^), the 
photons and baryons decouple and the mode becomes unstable again and grows as 5b oc 
a. Note that fiuctuations with size greater than Mj(acq) grow uninterruptedly (because 
after horizon crossing the photon pressure cannot balance the gravity), changing the 
trend from before MRE to a after it, while those with sizes smaller than Mg are 
completely washed out by photon diffusion. 

After decoupling all surviving perturbations (those with Mpcrt > Mg) grow steadily 
until their amplitude becomes of order unity or larger. At that point the linear theory 
breaks down and one needs to employ a different type of analysis. 

If we look, instead, at the schematic evolution of an adiabatic mode with the same 
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Figure 5.12: Typical evolution of a perturbed scale Mpert (black line) in adiabatic mirror 
baryonic dark matter scenario with x > Xeq- The figure shows the Jeans mass Mj (blue), the 
Silk mass Mg (red) and the Hubble mass Mh (green). The time of horizon crossing of the 
perturbation is indicated by Cent- Are also indicated the three evolutionary stages: during 
stage A (a < Oent < Oeq) the mode grows as 6^ oc a^; throughout stage B (oent < a < a^ec) 
the perturbation oscillates; finally, in stage C (a > a^ec) mode becomes unstable again 
and grows as 6^ oc a. Note that fluctuations with size smaller than Mg are wiped out by 
photon diffusion. 

reference mass scale Mpert but for x < Xeq, as reported in figure 5.13, we immediately 
notice the lower values of the maximum Jeans mass and the Silk mass, which are now 
similar. So, for the plotted perturbative scale there are now only the two stages A and 
C. In general, depending on its size, the perturbation mass can be higher or lower than 
the Silk mass (and approximatively also than the maximum Jeans mass), so modes 
with Mport > Ml. grow continuously before and after their horizon entry, while modes 
with Mpert < A^s completely washed out. 

Notice that Mj becomes smaller than the Hubble horizon mass Mh starting from 
a redshift 

= 3750 x-^ (n^h^) , (5.36) 

which is about z^q for x = 0.64, but it sharply increases for smaller values of x, as shown 
in figure 5.14. We can recognize this behaviour also watching at the intersections of the 
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Figure 5.13: Typical evolution of a perturbed scale Mpert (black line) in adiabatic mirror 
baryonic dark matter scenario with x < x^q. The value of Mpert is the same as in figure 5.12. 
The time of horizon crossing of the perturbation is indicated by Oent- The figure shows the 
Jeans mass M'j (blue), the Silk mass Mg (red) and the Hubble mass Mh (green). Unlike the 
case X > Xeq (shown in the previous figure), now there are only the two evolutionary stages 
A (a < Oent) and C (a > Oent)- Fluctuations with size smaller than Mg are wiped out by 
photon diffusion, but now the Silk mass is near to the maximum of the Jeans mass. 

lines for and Mh in figures 5.12 and 5.13. Thus, density perturbation scales which 
enter horizon at z ~ Zeq have mass larger than Mj and thus undergo uninterrupted 
linear growth immediately after t^q. Smaller scales for which > Mh would instead 
first oscillate. Therefore, the large scale structure formation is not delayed even if the 
mirror decoupling did not occur yet, i.e. even if a; > Xeq- 

When compared with the non baryonic dark matter scenarios (see § 2.4), the main 
feature of the mirror baryonic dark matter scenario is that the M baryon density 
fluctuations should undergo the strong collisional damping around the time of M re- 
combination, which washes out the perturbations at scales smaller than the mirror 
Silk scale. It follows that density perturbation scales which undergo the linear growth 
after the MRE epoch are limited by the length A^. This could help in avoiding the 
excess of small scales (of few Mpc) in the CDM power spectrum mentioned in § 2.4.4 
without tilting the spectral index. To some extent, the cutoff effect is analogous to 
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Figure 5.14: The M photon decouphng redshift 1 + z'^^^ as a function of x (solid line). The 
long-dash line marks the ordinary decoupling redshift 1 + Zdec = 1100. We also show the 
matter-radiation equality redshift 1 -|- Zeq and the mirror Jeans- horizon mass equality redshift 
1 -|- z^, for the cases Q.mh'^ = 0.2 (respectively lower dash and lower dot) and ^mh'^ = 0.6 
(upper dash and dot). 

the free streaming damping in the case of warm dark matter (WDM), but there are 
important differences. The point is that, alike usual baryons, the MBDM shows acous- 
tic oscillations with an impact on the large scale power spectrum. In particular, it is 
tempting to imagine that the M baryon oscillation effects are related to the anoma- 
lous features observed in large scale structure power spectra (see next chapter for a 
complete discussion). 

In addition, the MBDM oscillations transmitted via gravity to the ordinary baryons, 
could cause observable anomalies in the CMB angular power spectrum for /'s larger 
than 200. This effect can be observed only if the M baryon Jeans scale Aj is larger than 
the Silk scale of ordinary baryons, which sets a principal cutoff for CMB oscillations 
around / ~ 1200. As we will see in the next chapter, this would require enough large 
values of x, near the upper bound fixed by the BBN constraints, and, together with 
the possible effects on the large scale power spectrum, it can provide a direct test for 
the MBDM (verifiable by the higher sensitivity of next CMB and LSS experiments). 

Clearly, for small x the M matter recombines before the MRE moment, and thus it 
behaves as the CDM as far as the large scale structure is concerned. However, there 
still can be crucial differences at smaller scales which already went non-linear, like 
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galaxies. In our scenario, dark matter in galaxies and clusters can contain both CDM 
and MBDM components, or can be even constituted entirely by the mirror baryons. 

One can question whether the MBDM distribution in halos can be different from 
that of the CDM. Simulations show that the CDM forms triaxial halos with a den- 
sity profile too clumped towards the center, and overproduces the small substructures 
within the halo. Since MBDM constitutes a kind of coUisional dark matter, it may 
potentially avoid these problems, at least the one related with the excess of small 
substructures. 

It's also worth noting that, throughout the above discussion, we have assumed that 
the matter density of the Universe is close to unity. If, instead, the matter density is 
small and a vacuum density contribution is present, there is an additional complication 
due to the fact that the Universe may become curvature dominated starting from some 
redshift Zcurv Given the current estimate ft^ — 0.7 (see § 1.7), this transition has yet 
occurred and the growth of perturbations has stopped around Zcurv, when the expansion 
became too rapid for it. 

At the end, we spend few words to mention that the main difficulty with the 
ordinary baryonic adiabatic scenario is the excess of angular fluctuations in the CMB 
temperature respect to the observational limits. More specifically, one needs S^, ~ 10~^ 
at recombination, but in the adiabatic picture matter inhomogeneities are accompanied 
by perturbations in the radiation field, and this will inevitably lead to temperature 
fluctuations of order 5T/T ~ 5^ ~ 10~^ at decouphng, which is in direct disagreement 
with observations. In this sense the mirror baryonic adiabatic scenario has various ways 
to overcome this problem, first of all the fact that mirror recombination takes place 
x"^ times before the ordinary one (see § 4.6), and structures have more time to grow 
after decouphng, in a way compatible with the CMB observations of ST/T ~ 10~^. 

5.4 Evolution of perturbations 

As a result of the studies done in previous sections, in this section we finally consider 
the temporal evolution of perturbations, as function of the scale factor a. All the 
plots are the results of numerical computations obtained making use of a Fortran code 
originally written for the ordinary Universe and modified to account for the mirror 
sector (for more details see the beginning of next chapter). 

We used the synchronous gauge, described in § 3.3 and appendix A. 3. The difference 
in the use of other gauges is limited to the gauge-dependent behaviour of the density 
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fluctuations on scales larger than the horizon. The fluctuations can appear as growing 
modes in one coordinate system and as constant mode in another, that is exactly what 
occurs in the synchronous and the conformal Newtonian gauges. 

In the flgures we plot the evolution of the components of a mirror Universe, namely 
the cold dark matter^, the ordinary baryons and photons, and the mirror baryons and 
photons, changing some parameter to evaluate their influence. Note that in flgures are 
plotted the density contrasts not normalized to the average density. 

First of all, we comment flgure 5.15b, which is the most useful to recognize the 
general features of the evolution of perturbations. Starting from the smallest scale 
factor, we see that all three matter components and the two radiative components 
grow with the same trend (as a^), but the radiative ones have a slightly higher den- 
sity contrast (with a constant rate until they are tightly coupled); this is simply the 
consequence of considering adiabatic perturbations, which are linked in their matter 
and radiation components by the adiabatic condition (5.1). This is the situation when 
the perturbation is out of horizon, but, when it crosses the horizon, around a ~ 10~^, 
things drastically change. Baryons and photons, in each sector separately, become 
causally connected, feel each other, and begin to oscillate for the competitive effects 
of gravity and pressure. Meanwhile, the CDM density perturbation continues to grow 
uninterruptedly, at flrst reducing his rate from to In a (due to the rapid expansion 
during the radiation era), and later, as soon as MRE occurs (at a ~ 3 x 10~^ for 
the considered model), increasing proportionally to a. The oscillations of baryons and 
photons continue until their decoupling, which in the mirror sector occurs before than 
in the ordinary one (scaled by the factor x) . This moment is marked in the plot as the 
point where the lines for the two components move away one from the other. Prom 
this point, the photons continue the oscillations until they are completely damped, 
while the baryons rapidly fall into the potential wells created by the dark matter and 
begin to grow as a. Note that it's important the way in which the oscillation reaches 
the decoupling; if it is compressing, flrst it slows down (as if it continues to oscillate, 
but disconnected from the photons), and then it expands driven by the gravity; if, 
otherwise, it is expanding, it directly continues the expansion and we see in the plot 
that it immediately stops to oscillate. In this flgure we have the flrst behaviour in the 
mirror sector, the second one in the ordinary sector. 

In flgures 5.15, 5.16 and 5.17 we compare the behaviour of different scales for the 
same model. The scales are given by log A; = —0.5, —1.0, —1.5, —2.0, —2.5, —3.0, where 

^As non baryonic dark matter we consider only the cold dark matter, which is now the standard 
choice in cosmology. 
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/c = 27r/A is the wave number and is measured in Mpc~^. The effect of early entering 
the horizon of small scales (those with higher k) is evident. Going toward bigger 
scales the superhorizon growth continue for a longer time, delaying more and more 
the beginning of acoustic oscillations, until it occurs out of the coordinate box for the 
bigger plotted scale (log A; = —3.0). Starting from the scale log/c = —1.5, the mirror 
decoupling occurs before the horizon entry of the perturbations, and the evolution 
of the mirror baryons density is similar to that of the CDM. The same happens to 
the ordinary baryons too, but for log A; ^ —2.0 (since they decouple later), while the 
evolution of mirror baryons is yet indistinguishable from that of the CDM. For the 
bigger scales (log A; <^ —2.5) the evolution of all three matter components is identical. 

As previously seen, the decoupling is a crucial point for structure formation, and it 
assumes a fundamental role specially in the mirror sector, where it occurs before than 
in the ordinary one: mirror baryons can begin before to grow perturbations in their 
density distribution. For this reason it's important to analyze the effect of changing 
the mirror decoupling time, obtained changing the value of x and leaving unchanged 
all other parameters, as it is possible to do using figures 5.15b, 5.18 and 5.19 for 
X — 0.6,0.5,0.4,0.3,0.2 and the same scale log A; = —1.0. It is evident the shift of 
the mirror decoupling toward lower values of a when reducing x, according to the law 
(4.51), which states a direct proportionality between the two. In particular, for x < 0.3 
mirror decoupling occurs before the horizon crossing of the perturbation, and mirror 
baryons mimic more and more the CDM, so that for x ~ 0.2 the perturbations in 
the two components are indistinguishable. For the ordinary sector apparently there 
are no changes, but at a more careful inspection we note some difference due to the 
different amount of relativistic mirror species (proportional to x^), which slightly shifts 
the matter-radiation equality. This effect is more clear in figure 5.20, where we plot 
only the CDM and the ordinary baryons for the cases x = 0.2 and 0.6: for the lower 
value of x there are less mirror photons, the MRE occurs before and the perturbation 
in the collisionless component starts growing proportionally to the scale factor before; 
thus, when the baryons decouple, their perturbation rapidly grows to equalize that in 
the CDM, which meanwhile has raised more for the lower x. 

Obviously, these are cases where the CDM continues to be the dominant form of 
dark matter, and drives the growth of perturbations, given its continuous increase. In 
any case, if the dominant form of dark matter is made by mirror baryons the situation 
is practically the same, as visible comparing figures 5.18b and 5.21a (where we see only 
slight differences on the CDM and mirror baryons behaviours in the central region of 
the plots), since mirror baryons decouple before than ordinary ones and fall into the 
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potential wells of the CDM, reinforcing them. 

Finally, in the interesting case where mirror baryons constitute all the dark matter, 
they drive the evolution of perturbations. In fact, in figure 5.21b we clearly see that 
the density fluctuations start growing in the mirror matter and the visible baryons 
are involved later, after being recombined, when they rewrite the spectrum of already 
developed mirror structures. This is another effect of a mirror decoupling occurring 
earlier than the ordinary one: the mirror matter can drive the growth of perturbations 
in ordinary matter and provide the rapid growth soon after recombination necessary 
to take into account of the evolved structures that we see today. 

Given all the considerations made in this chapter, it is evident that the case of 
mirror baryons is very interesting for structure formation, because they are coUisional 
between themselves but coUisionless for the ordinary sector, or, in other words, they 
are self-coUisional. In this situation baryons and photons in the mirror sector are 
tightly coupled until decoupling, and structures cannot grow before this time, but the 
mirror decoupling happens before the ordinary one, thus structures have enough time 
to grow according to the hmits imposed by CMB and LSS (something not possible in 
a purely ordinary baryonic scenario). Another important feature of the mirror dark 
matter scenario is that, if we consider small values of x, the perturbation evolution 
is very similar to the CDM case, but with a fundamental difference: there exist a 
cutoff scale due to the mirror Silk damping, which kills the small scales, overcoming 
the problems of the CDM scenario with the excessive number of small sateUites. These 
are important motivations to go further in the work and investigate the effects of the 
mirror sector on the CMB and LSS power spectra, as we will do in the next chapter. 
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Figure 5.15: Evolution of perturbations for the components of a Mirror Universe: cold 
dark matter (dark line), ordinary baryons and photons (solid and dotted red lines) and 
mirror baryons and photons (solid and dotted blue lines). The model is a flat Universe with 
= 0.3, Qb^^ = 0.02, n'^h"^ = 0.04 {p = 2), h = 0.7, x = 0.6, and plotted scales are 
logA:(Mpc~i) = -0.5 (a) and -1.0 (6). 
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Figure 5.16: The same as in figure 5.15, but for scales logA;(Mpc ^) = 
{b). 



-1.5 (a) and -2.0 
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Figure 5.18: Evolution of perturbations for the components of a Mirror Universe: cold 
dark matter (dark line), ordinary baryons and photons (solid and dotted red lines) and 
mirror baryons and photons (solid and dotted blue lines). The model is a flat Universe with 
n^a = 0.3, 17b/i^ = 0.02, = 0.04 (/3 = 2), h = 0.7, x = 0.5 (a) or 0.4 (6), and plotted 

scale is logfc(Mpc~^) = —1.0. 
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Figure 5.19: The same as in figure 5.18, but for x = 0.3 (a) and 0.2 (6). 
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Figure 5.20: Evolution of perturbations in a Mirror Universe for cold dark matter (dark 
lines) and ordinary baryons (red and green lines). The models are flat, with Jim = 0.3, 
Ob/i^ = 0.02, n'y^h'^ = 0.04 {p = 2), h = 0.7, x = 0.2 (solid black and red lines) and 0.6 
(dashed black and green lines), and log ^(Mpc^^) = —1.0. 
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Figure 5.21: Evolution of perturbations for the components of a Mirror Universe: cold dark 
matter (dark line), ordinary baryons and photons (solid and dotted red lines) and mirror 
baryons and photons (solid and dotted blue lines). The models are flat with = 0.3, 
^b^^ = 0.02, n'^ = (2/3)(17m-Ob) (/3 = 4) (a) or (O^-Ob) (no CDM) (6), /i = 0.7, x = 0.4, 
and logA;(Mpc"^) = -1.0. 



Chapter 6 

Cosmic microwave background and 
large scale structure for a Mirror 
Universe 



6.1 Introduction 

Let us now consider the consequences of the existence of a mirror sector in terms of 
signatures on the Cosmic Microwave Background (CMB) and the Large Scale Structure 
(LSS). In the last decade their study is providing a great amount of observational data, 
and the continuous improvement is such that we are in the condition to verify different 
models of Universe, with the possibility to fit the cosmological parameters (for recent 
works see, e.g., Wang et al. (2002) [180] or Efstathiou ct al. (2002) [57]). These 
powerful cosmological instruments could help us to understand the nature of the dark 
matter of the Universe analyzing the implications of different types of dark matter on 
their cosmological observables, and comparing them with the available experimental 
results. In this view we must study the consequences of the Mirror Universe on the 
CMB and LSS in order to compare our models with data, study their compatibility 
and possibly reduce the available parameter space. 

To do this we need to numerically compute mirror models and their CMB and LSS 
power spectra (see § 3.1.1 and § 3.2). We used the evolutionary equations written 
by Ma and Bertschinger in 1995 [118] and described in § 3.3, inserted in a Fortran 
code developed by the group of the prof. N. Vittorio of the University of Rome "Tor 
Vergata" and kindly provided to us. Obviously, this program was written for a standard 
Universe made only by the ordinary sector, and the first step of our work was to modify 
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the code in order to simulate a Mirror Universe (which, as we know now, has two 
sectors instead). This required to modify and add several subroutines to the code. In 
fact, according to what said in previous chapters, we have to treat two self-interacting 
sectors, and thus double all the equations governing the evolution of the considered 
components, relativistic (photons and massless neutrinos) and non relativistic (baryons 
and cold dark matter). In addition, the two sectors communicate via gravity, so that 
they are coupled and influence each other through this interaction. Therefore, there 
are more regimes respect to the standard case, because now, instead of being simply 
the ordinary ones, they are made up of the couplings between the different ordinary 
and mirror regimes, the latter being time-shifted from the former according to the laws 
exposed in § 4.6 and § 5.2.1. 



We computed many models for Mirror Universe, assuming adiabatic scalar pertur- 
bations, a flat space-time geometry, and different mixtures of ordinary and mirror 
baryons, photons and massless neutrinos, cold dark matter, and cosmological constant. 

In Fourier space, all the k modes in the linearized Einstein, Boltzmann, and fluid 
equations evolve independently; thus the equations can be solved for one value of k 
at a time. Moreover, all modes with the same k (the magnitude of the comoving 
wavevector,) obey the same evolutionary equations. We integrated the equations of 
motion numerically over the range —5.0 < log A; < —0.5 (where k is measured in Mpc~^) 
using points evenly spaced with an interval of A log A; = 0.01. The full integration was 
carried to z — 0. 

As shown in § 3.3, the Boltzmann equation for massless neutrinos (3.33) has been 
transformed into an inflnite hierarchy of moment equations that must be truncated 
at some maximum multipole order Zmax- One simple but inaccurate method is to set 
Fi,i — ioT I > Imax-^ According to what suggested by Ma and Bertschinger (1995) 
[118] instead, an improved truncation scheme is based on extrapolating the behaviour 

of F^i to I = Zmax + 1 as 



The problem with this scheme is that the couphng of multipoles in equations leads to the propaga- 
tion of errors from /max to smaller I. Indeed, these errors can propagate to Z = in a time r w Zmax/fc 
and then reflect back to increasing I, leading to amplification of errors in the interval < I < /max- 
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However, time-variations of the potentials during the radiation-dominated era make 
even equation (6.1) a poor approximation if Zmax is chosen too small. 

Thus, in order to have a relative accuracy better than 10~^ in our final results, in 
the computation of the potential and the density fields the photons and the massless 
neutrinos phase space distributions for both the ordinary and mirror sectors were 
expanded in Legendre series truncating the Boltzmann hierarchies at Zmax = 2000, 
using the truncation schemes given by equations (6.1) and (3.38). 

At last, we normahze our results to the COBE data, following the procedure de- 
scribed by Bunn and White in 1997 [37]. 



Parameter 


min. value 


max. value 




0.1 


1.0 




0.010 


0.030 


< 


0.0 


— 


X 


0.1 


1.0 


h 


0.5 


0.9 


n 


0.90 


1.10 



Table 6.1: Parameters and their ranges used in mirror models. The values are not evenly 
spaced, but arbitrarily chosen in the parameter space. Not listed there are the total and 
vacuum densities, but, being fiat models, they are Q,q = 1 and Ha = 1 — Om- 

We considered different values of the cosmological parameters, where now we add 
to the usual ones two new mirror parameters: the ratio of the temperatures in the two 
sectors x (see § 4.2) and the mirror baryons density 0!^^ (also expressed via the ratio of 
the baryonic densities in the two sectors (3, defined in § 4.3). Starting from an ordinary 
reference model (we choose the so-called "concordance model" of Wang et al. (2002) 
[180]), we study the influence of the mirror sector varying the two parameters that 
describe it for a given ordinary sector. Furthermore, we evaluate the influence of all 
the parameters for a Mirror Universe, changing all of them. The values used for the 
parameters are not on a regular grid, but arbitrarily chosen for the only purpose to 
better understand the CMB and LSS for a Mirror Universe. At the moment we are not 
able to make a grid thin enough to perform a fit of all free parameters (even if we fix 
some of them) for two reasons: first, we have two further parameters, which lengthen 
a lot the computational time (by a factor 10^), and second our program is much slower 
than others commonly used for a standard Universe, because our models are more 
complicated in terms of calculus and for our choice to privilege precision instead of 
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performance (at least in this moment^). We list the parameters used and their value 
ranges in table 6.1 (remember the comment on the grid choice); the total and vacuum 
densities (not listed in the table) are fixed by our choice of a flat geometry: — 1 
and D,\ — 1 — In addition, in order to study the parameter dependence, we did 
computations also for different numbers of extra-neutrino species (AA^,^ = 0.5, 1.0, 1.5). 

6.3 The cosmic microwave background for a Mirror 
Universe 

As anticipated when we studied the structure formation, we expect that the existence 
of a mirror sector influences the cosmic microwave background radiation observable 
today; now we want to evaluate this effect. 

We choose a starting standard model and add a mirror sector simply removing cold 
dark matter and adding mirror baryons. At present, we think that a good reference 
model is the so-called "concordance model" of Wang et al. (2002) [180] (in this case 
they performed an 11-parameters flt of all the available experimental data). We are 
aware that for the reader this model will be surely obsolete, given the current rate of 
new experimental data. However, this is not a shortcoming, because here we want only 
to put in evidence the differences respect to a representative reference model, and this 
is a good model for this purpose, unless new (very unlikely) revolutionary observations 
will change the scenario. The parameter values for this reference model are: — 1, 
= 0.34, Qa = 0.66, lo^ ^ Vt^h? ^ 0.02, = 0.91, h = 0.64, with only cold dark 
matter (no massive neutrinos) and scalar adiabatic perturbations. 

Prom this starting point, we flrst substitute all the cold dark matter with mirror 
baryonic dark matter (MBDM) and evaluate the CMB angular power spectrum varying 
X from 0.3 to 0.7. This is shown in top panel of flgure 6.1, where mirror models are 
plotted together with the concordance model. The flrst evidence is that the deviation 
from the standard model is not linear in x: it grows more for bigger x and for a; < 0.3 
the power spectra are practically coincident. This is important, because it means that a 
Universe where all the dark matter is made of mirror baryons could be indistinguishable 
from a CDM model if we analyse the CMB only. We see the greatest separation from 

^One of the future steps of the work will be just to write a new program (probably based on the 
commonly used Seljak and Zaldarriaga's CMBFAST code [160, 186]) much faster than the one now 
used, which should allow us to fit the parameters and compare the results with the other cosmological 
models. 
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Figure 6.1: CMB angular power spectrum for different values of x and a;(, = compared 
with a standard model (solid line) . Top panel. Mirror models with the same parameters as 
the ordinary one, and with x = 0.3,0.5,0.7, and a;(, = Om^^ — ^^b (no CDM) for all models. 
Bottom panel. Mirror models with same parameters as the ordinary one, and with x = 0.7 

and Lo[^ = <^b) 2c<;b! 3wb) '^^h- 
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the standard model for x — 0.7, but it will increase for larger values of x. The height 
of the first acoustic peak grows for x > 0.3, while the position remains nearly constant. 
For the second peak occurs the opposite, i.e. the height remains practically constant, 
while the position shifts toward higher I; for the third peak, instead, we have a shift 
both in height and position (the absolute shifts are similar to the ones for the first two 
peaks, but the height now decreases instead of increasing). Observing also the other 
peaks, we recognize a general pattern: except for the first one, odd peaks change both 
height and location, even ones change location only. 

In bottom panel of figure 6.1 we show the intermediate case of a mixture of CDM 
and MBDM. We consider x — 0.7, a high value which permits us to see well the 
differences, and change u;'^ from uy^ (20% of Jldm) to Acoy^ (80% of fidm)- The dependence 
on the amount of mirror baryons is lower than on the ratio of temperatures x. In fact, 
the position of the first peak is nearly stable for all the mirror models (except for a very 
low increase of height for growing a;(,), while differences appear for the other peaks. 
In the second peak the position is shifted as in the case without CDM independently 
from while the height is inversely proportional to cu'^ with a separation appreciable 
for < Scuy^. For the third peak the behaviour is the same as for the case without 
CDM, with a shghtly stronger dependence on lo^, while for the other peaks there is 
a weaker dependence on a;(,. A common feature is that the heights of the peaks are 
not linearly dependent on cu'^, while their positions are practically insensitive to cu'^ but 
depend only on x. 

We will analyse in more detail the x and cu'^ dependence of the peaks, together with 
other parameters, in § 6.5. 

6.3.1 The mirror cosmic microwave background radiation 

In the same way as ordinary photons at decoupling from baryons formed the CMB 
we observe today, also mirror photons at their decoupling formed a mirror cosmic mi- 
crowave background radiation, which, on the contrary, we cannot observe because they 
don't couple with the ordinary baryons of which we are made^ (it would be possible 
for an hypothetical mirror observer, instead). Nevertheless, its study is not only spec- 
ulative, being a way to better understand the cosmology of the Mirror Universe and 
our observable CMB. 

•^Indeed, there is in principle the possibihty of an influence of the mirror CMB on our photons 
in case of existence of a photon-mirror photon mixing (as supposed, for example, by Foot (2002) 
[69] and references therein), but its detection is not possible with present and probably even future 
experiments. 
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Figure 6.2: Angular power spectra for ordinary (solid line) and mirror (dashed line) CMB 
photons. The models have = !> = 0.3, = = 0.02, h = 0.7, n = 1.0, and x = 0.7 
(top panel) and x = 0.5 (bottom panel). 
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We computed two models of mirror CMB, in order to have elements to compare with 
the corresponding observable CMBs. The chosen parameter values are those currently 
estimated (see § 1.7), and the amount of mirror baryons is the same as the ordinary 
ones, while x is taken as 0.7 or 0.5 to explore different scenarios. Thus the parameters 
of the models are: = 1, — 0.3, u;b — u;{j — 0.02, x — 0.7 or 0.5, h — 0.7, n — 1.0. 
In figure 6.2 we plot the ordinary and mirror CMB spectra corresponding to the same 
Mirror Universe. 

The first evidence is that, being scaled by the factor x the temperatures in the 
two sectors, also their temperature fluctuations will be scaled by the same amount, as 
evident if we look at the lowest i values (the fluctuations seeds are the same for both 
sectors) . Starting from the top flgure, we see that the first mirror CMB peak is much 
higher and shifted to higher multipoles than the ordinary one, while other peaks are 
both lower and at higher i values, with a shift growing with the order of the peak. 

Observing the bottom plot, we note the effect of a change of the parameter x on the 
mirror CMB: (i) for lower x-values the first peak is higher (for x — 0.5 it is nearly one 
and half the ordinary one) ; (ii) the position shifts to much higher multipoles (with the 
same horizontal scale we can no more see some peaks) . The reason is that a change of 
X corresponds to a change of the mirror decoupling time. The mirror photons, which 
decouple before the ordinary ones, see a smaller sound horizon, scaled approximately 
by the factor x; since the first peak occurs at a multipole i oc (sound horizon)"^, we 
expect it to shifts to higher ^-values by a factor x~^, that is exactly what we see in the 
figure. 

We have verified (even if not shown in the figures) that increasing x the mirror 
CMB is more and more similar to the ordinary one, until for x — 1 the two power 
spectra are perfectly coincident (as expected, since in this case the two sectors have 
exactly the same temperatures, the same particle contents, and then their photons 
power spectra are necessarily the same). 

If we were able to detect both the ordinary and mirror CMB photons, we had two 
snapshots of the Universe at two different epochs, which were a powerful cosmolog- 
ical instrument, but unfortunately this is impossible, because mirror photons are by 
definition completely invisible for us. 



§6.4 The large scale structure for a Mirror Universe 



119 



6.4 The large scale structure for a Mirror Universe 

Given the oscillatory behaviour of the mirror baryons (different from the smooth one 
of cold dark matter), we expect that MBDM induces specific signatures also on the 
large scale structure power spectrum (see § 3.2). 

In order to evaluate this effect, we computed LSS power spectra using the same 
reference and mirror models used in § 6.3 for the CMB analysis. So the two panels 
of figure 6.3 are the LSS corresponding to those of figure 6.1. In order to remove the 
dependences of units on the Hubble constant, we plot on the x-axis the wave number 
in units of h and on the y-axis the power spectrum in units of h^. The minimum scale 
(the maximum k) plotted depends on the limit of the linear regime, placed between 
k/h — 0.3 and 0.4 Mpc~^, according to what described by Hamilton & Tegmark (2000) 
[88]. 

In top panel of the figure we show the dependence on x for different mirror models 
without CDM; in this case, where all the dark matter is made of mirror baryons, the 
oscillatory effect is obviously maximum. The first evidence is the strong dependence 
on X of the beginning of oscillations: it goes to higher scales for higher x and, starting 
from X ~ 0.3, the power spectrum for a Mirror Universe approaches more and more 
the CDM one. This behaviour is a consequence of the x-depcndence of the mirror Silk 
scale (see § 5.2.4 and figure 5.11): this dissipative scale induces a cutoff in the power 
spectrum, which is damped with an oscillatory behaviour (it will be more evident in 
figures 6.4 and 6.5, where we extend our models to smaller scales inside the non linear 
region). Oscillations begin at the same time of the damping, and they are so deep 
(because there are many mirror baryons) to go outside the coordinate box. In any case 
the mirror spectra are always below the ordinary one for every value of x. 

The dependence on the amount of mirror baryons is instead shown in the bottom 
panel of the figure, where only a fraction of the dark matter is made of mirror baryons, 
while the rest is CDM. Contrary to the CMB case, the matter power spectrum strongly 
depends on a;^. The oscillations are deeper for increasing mirror baryons densities and 
the spectrum goes more and more away from the pure CDM one. We note also that 
the damping begins always at the same scale, and thus it depends only on x and not 
on a;[,, as we know from expression (5.33). The same considerations are valid for the 
oscillation minima, which become much deeper for higher mirror baryon densities, but 
shift very slightly to lower scales, so that their positions remain practically constant. 
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Figure 6.3: LSS power spectrum in the linear regime for different values of x and oo'^ = ^'^h'^, 
compared with a standard model (solid line). In order to remove the dependences of units on 
the Hubble constant, we plot on the x-axis the wave number in units of h and on the y-axis 
the power spectrum in units of h^. Top panel. Mirror models with the same parameters as 
the ordinary one, and with x = 0.3, 0.5, 0.7 and uj[^ = ^mh'^ — oj^ (no CDM) for all models. 
Bottom panel. Mirror models with the same parameters as the ordinary one, and with x = 0.7 
and Lu[^ = 2a;b, Swb, 4a;b- 
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6.4.1 Extension to smaller (non linear) scales 

Let us now extend the behaviour of the matter power spectrum to lower scales, which 
already became non linear. Obviously, since our treatment is based on the linear 
theory, it is no longer valid in non linear regime. Nevertheless, even if it is not useful 
for comparison with observations, the extension of our models to these scales is very 
useful to understand the behaviour of the power spectrum in a mirror baryonic dark 
matter scenario, in particular concerning the position of the cutoff (remember the 
discussion in § 2.4.4 about the failure of the CDM scenario). 

Therefore, in figures 6.4 and 6.5 we extend the power spectra up to k/h = 10 
Mpc^^ (corresponding to galactic scales), well beyond the limit of the linear regime, 
given approximately by k/h < 0.4 Mpc"-"^. 

In figure 6.4 we plot in both panels the same models as in figure 6.3, except for the 
spectral index now set to 1.0, a value chosen to eliminate the effect of a power spectrum 
tilt on the cutoff.^ For comparison we show also a standard model characterized by 
a matter density made almost completely of CDM, with only a small contamination 
of baryons {fl^ — 0.2% instead of ~ 4% of other models). In the top panel, the x- 
dependence of the mirror power spectra is considered: the vertical scale extends to 
much lower values compared to figure 6.3, and we can clearly see the deep oscillations, 
but in particular it is evident the presence of the previously cited cutoff. For larger 
values of x oscillations begin earlier and cutoff moves to higher scales. Moreover, 
note that the model with almost all CDM has more power than the same standard 
model with baryons, which in turn has more power than all mirror models for any x 
and for all the scales. In the bottom panel we show the dependence on the baryon 
content. It is remarkable that all mirror models stop to oscillate at some low scale 
and then continue with a smooth CDM-like trend. This means that, after the cutoff 
due to mirror baryons, the dominant behaviour is the one characteristic of cold dark 
matter models (due to the lack of a cutoff for CDM) . Clearly, for higher mirror baryon 
densities the oscillations continue down to smaller scales, but, contrary to the previous 
case, where all the dark matter was mirror baryonic, there will always be a scale below 
which the spectrum is CDM-like. 

An interesting point of the mirror baryonic scenario is his capability to mimic a 
CDM scenario under certain circumstances and for certain measurements. To explain 
this point, in figure 6.5 we show models with low x- values (0.2 or 0.1) and all dark 

^Indeed, here we want to simply compare the standard and mirror power spectra, and a common 
tilt for both models does not influence our conclusions; nevertheless, we prefer to use a scale invariant 
spectrum. 
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Figure 6.4: LSS power spectrum beyond the linear regime for different values of x and uo'^ = 
^'^h"^, compared with a standard model (solid line). The models have the same parameters 
as in figure 6.3, except for the spectral index, which is now set to 1.0. For comparison we 
also show a standard CDM model with a negligible amount of baryons (ilb ^ 0.2%). 
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matter made of MBDM; we see that for x — 0.2 the standard and mirror power spectra 
are aheady practically coincident in the linear region. If we go down to x — 0.1 the 
coincidence is extended up to k/h ~ 1 Mpc"^. In principle, we could still decrease 
X and lengthen this region of equivalence between the different CDM and MBDM 
models, but we have to remember that we are dealing with linear models extended 
to non linear scales, so neglecting all the non linear phenomena (such as merging or 
stellar feedback), that are very different for the CDM and the MBDM scenarios. In the 
same plot we also considered a model with x — 0.2 and dark matter composed equally 
by mirror baryons and by CDM. This model shows that in principle it's possible a 
tuning of the cutoff effect reducing the amount of mirror matter, in order to better 
reproduce the cutoff needed to explain, for example, the low number of small satellites 
in galaxies. 

This work provided for us the linear transfer functions, which constitute the prin- 
cipal ingredient for the computation of the power spectrum at non linear scales. This 
calculation is out of the aim of this thesis, but is one of the next steps in the study of 
the Mirror Universe. 




Figure 6.5: LSS power spectrum beyond the linear regime for two low values of x (0.1 and 
0.2) and different amounts of mirror baryons (W^^ = — or = ^Icbm), compared with 
a standard model (solid line). The other parameters are the same as in figure 6.3, except for 
the spectral index, which is now set to 1.0. For comparison we show also a standard CDM 
model with a negligible amount of baryons (J^b ~ 0.2%). 
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6.5 Parameter trends 

The shapes, heights and locations of peaks and oscillations in the photons and matter 
power spectra are predicted by all models based on the inflationary scenario. Further- 
more the details of the features in these power spectra depend critically on the chosen 
cosmological parameters, which in turn can be accurately determined by precise mea- 
surements of these patterns. In this section we briefly discuss the sensitivity of the 
CiS and P(A;)'s on the values of some fundamental parameters in a mirror baryonic 
scenario. 

In particular, the exact form of the CMB and LSS power spectra is greatly depen- 
dent on assumptions about the matter content of the Universe. Apart from the total 
density parameter Qq (our models are flat, so Qq is always 1), the composition of the 
Universe can be parametrized by its components and and the components of 
the matter density f2b, i^b, f^cDM- Further parameters are the tilt of scalar fluctua- 
tions n, the Hubble parameter h, and the ratio of temperatures in the two sectors x. In 
addition, we consider also the dependence on the number of massless neutrino species 
A*";^, in order to compare it with the x-dependence. This is important if we remember 
that the relativistic mirror particles can be parametrized in terms of effective number 
of extra-neutrino species (see § 4.2). 

Starting from the reference model of parameters = 0.3, U]^ = uj{^ = 0.02, 
X = 0.2, h = 0.7 and n = 1.0, we change one parameter each time, compute the 
respective models, and plot the CMB and LSS power spectra in order to show the 
dependence on it (flgures 6.6 - 6.13). Then, we compute the relative locations and 
heights of the flrst three acoustic peaks of the CMB angular power spectrum and plot 
them in flgures 6.14 and 6.15 to compare their sensitivities to the parameters. 

In the following we briefly analyse the dependences on every parameter, referring 
to the figure where the respective models are plotted. 

• Matter density (fig. 6.6): varies from 0.1 to 0.5. In flat models a decrease in 
Jim implies two things: an increase in Qa (with the consequent delay in matter- 
radiation equality) and a decrease in JIcdm (if we leave unchanged the O and M 
baryon densities). Both these things correspond to boosting and shifting effects 
on the acoustic peaks, while the matter power spectrum goes down, given the 
decreasing in the coUisionless species (CDM) and the progressive relative growth 
of the baryon densities, which are responsible for the oscillatory features. 

• O baryon density (fig. 6.7): a;b = u!{^ varies from 0.01 to 0.03. An increase 
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of the baryon fraction increases odd peaks (compression phase of the baryon- 
photon fluid) due to extra-gravity from baryons with respect to the even peaks 
(rarefaction phase of the fluid oscillation) in the CMB, and generate deeper 
oscillations in the LSS. In particular, the relative magnitudes of the first and 
second acoustic peaks are sensitive to a;b, as we see in figure 6.15. These effects 
are completely due to O baryons. In fact, even if not shown, we have verified that 
an increase of cuy^ with a constant cu'^ has exactly the same consequences for this 
value of X, while to see the effect of M baryons we have to raise the temperature 
of the mirror sector. 

Hubble constant (figs. 6.8 and 6.9): h varies from 0.50 to 0.90, but now we can 
leave constant either fib = 6.8) or co'b = uj{^ (fig. 6.9). In both cases a 

decrease in h corresponds to a delay in the epoch of matter-radiation equality and 
to a different expansion rate. This boosts the CMB peaks and slightly changes 
their location toward higher ts (similar to the effect of an increase in Qa), and 
induces a decrease in the LSS spectrum. There are slight differences between 
the two situations of Q or a; constant, evident in particular on the first acoustic 
peak and on the matter oscillations. If fact, when we consider a;b,b' constant, the 
baryon densities flh,h' = '^b,b'/^^ grow for decreasing h, then favouring the raise 
of the first peak in the CMB and the onset of oscillations in the LSS. 

Spectral index (fig. 6.10): n varies from 0.90 to 1.10. Increasing n will raise the 
power spectra at large ts with respect to the low ^'s and at large values of k 
with respect to low values. This is not so evident in figure (except before the 
first acoustic peak) , where the curves seem nearly parallel as if they were simply 
vertically shifted; this means a low sensitivity to the spectral index in this range, 
as also evident in figures 6.14 and 6.15 for the CMB. 

Extra-neutrino species (fig. 6.11): AN^y varies from 0.0 to 1.5. The effect of 
increasing the number of massless neutrino species is a slow raise of the first 
acoustic peak and a shift to higher i values for next peaks, together with a shght 
lowering of the matter power spectrum; all these changes are nearly proportional 
to AN,y, as shown also in figures 6.14 and 6.15. 

Ratio of temperatures (fig. 6.12): x varies from 0.2 to 0.7. Concerning the CMB, 
the effect of raising x is qualitatively the same as an increase in AN„, but more 
pronounced (for these "cosmologically compatible" ranges) and with a non-linear 
dependence. In the LSS spectrum, instead, the situation is different from the case 
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of extra-neutrino species, as now a growth of x induces the onset of the oscillatory 
features at lower values of k. This behaviour has been studied in more details in 
§ 6.3 and § 6.4. 

• M baryon density (fig. 6.13): uj'^ varies from 0.01 (cJb / 2) to 0.08 (4cJb)- The 
value of X is now raised to 0.7, because for 0.2 there aren't differences between 
models with different uj{^ values (we start observing small deviations only for the 
higher A;- values in the matter power spectrum). Also this behaviour has been 
studied in more details in § 6.3 and § 6.4; here we want to emphasize the low 
sensibility of the CMB on cj(, (with a slightly stronger dependence starting from 
the third peak) and, on the contrary, the high sensitivity of the LSS. For the first 
one, an increase in u'^ causes a very low increase of the height of the first peak 
and a progressive more pronounced decrease for the next peaks, while for the 
second one there is a fast deepening of the oscillations, slightly changing their 
locations. 

In figures 6.14 and 6.15 we focus our attention on the CMB first three peaks, choos- 
ing some indicator which could quantify the sensitivity to the parameters previously 
discussed. In figure 6.14 we analyse the locations of the peaks, plotting the differences 
of the locations between the various models and the reference model for the three 
peaks; in figure 6.15 we plot the deviations of the differences between the heights of 
the peaks from the same quantities obtained for the reference model. In this way we 
obtain a clear picture of the trends of these indicators varying the parameters. These 
plots provide a useful reference in order to evaluate the influence of each parameter on 
the CMB and LSS power spectra, and they contain a number of informations; we can 
extract some of them particularly worth of noting. 

Looking at the locations, we see a great sensitivity on the matter density and the 
Hubble constant, and a negligible one on the spectral index and the amount of mirror 
baryons. Concerning the extra-neutrino species and the temperature of the sectors, 
the sensitivities are comparable, but the trends are different: they are respectively a 
constant slope for AA^,^ and an increasing one for x. 

As regards the peak temperatures, the most sensitive parameter, besides and h, 
is cuy^; the dependence on n and is a bit greater than what it is for the locations, and 
the differences between the trends with A^,^ and x are slightly more evident, specially 
for values x > 0.6. 
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Figure 6.6: Dependence of the shape of photon and matter power spectra on the matter 
density ilm- The reference model (solid line) has: Q^a = 0.3, = = 0.02, x = 0.2, 
h = 0.7 and n = 1.0. For other models, all the parameters are unchanged except for the one 
indicated: = 0.1 (dot line), = 0.2 (dash line), Jim = 0.4 (long dash line), = 0.5 
(dot-dash line). 
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Figure 6.7: Dependence of the shape of photon and matter power spectra on the ordinary 
baryon density with = ^'h- The reference model (solid line) has: Om = 0.3, = wj^ = 
0.02, X = 0.2, h = 0.7 and n = 1.0. For other models, all the parameters are unchanged 
except for the one indicated: co^ = 0.010 (dot line), lo\) = 0.015 (dash line), uj\, = 0.025 (long 
dash line), = 0.03 (dot-dash line). 
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Figure 6.8: Dependence of the shape of photon and matter power spectra on the Hubble 
parameter h with = Q'-^ = const. The reference model (solid line) has: ilm = 0.3, 
= n'^ = 0.0408 (the value obtained for cob = 0.02 and h = 0.7), x = 0.2, h = 0.7 and 
n = 1.0. For other models, all the parameters are unchanged except for the one indicated: 
h = 0.5 (dot line), h = 0.6 (dash line), h = 0.8 (long dash line), h = 0.9 (dot-dash line). 
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Figure 6.9: Dependence of the shape of photon and matter power spectra on the Hubble 
parameter h with lo^ = io[^ = const. The reference model (solid line) has: = 0.3, 
= cj{, = 0.02, X = 0.2, h = 0.7 and n = 1.0. For other models, all the parameters are 
unchanged except for the one indicated: h = 0.5 (dot line), h = 0.6 (dash line), h = 0.8 (long 
dash line), h = 0.9 (dot-dash line). 
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Figure 6.10: Dependence of the shape of photon and matter power spectra on the scalar 
spectral index n. The reference model (solid line) has: rim = 0.3, = = 0.02, x = 0.2, 
h = 0.7 and n = 1.0. For other models, all the parameters are unchanged except for the 
one indicated: n = 0.90 (dot line), n = 0.95 (dash line), n = 1.05 (long dash line), n = 1.10 
(dot-dash line). 
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Figure 6.11: Dependence of the shape of photon and matter power spectra on the number of 
extra-neutrino species AN^,. The reference model (sohd line) has: Q^a = 0.3, lox, = io'^ = 0.02, 
X = 0.2, h = 0.7, n = 1.0, and AiVj^ = 0. For other models, all the parameters are unchanged 
except for the one indicated: ANi, = 0.5 (dot line), ANi, = 1.0 (dash line), ANi, = 1.5 (long 
dash line). 



§6.5 Parameter trends 



133 



80 

70 : 



0.01 



X 





0.10 



k (hMpc ') 



Figure 6.12: Dependence of the shape of photon and matter power spectra on the ratio 
of the temperatures of two sectors x. The reference model (sohd hne) has: il^ = 0.3, 
Wb = = 0.02, X = 0.2, h = 0.7 and n = 1.0. For other models, all the parameters are 
unchanged except for the one indicated: x = 0.3 (dot line), x = 0.4 (dash line), x = 0.5 (long 
dash line), x = 0.6 (dot-dash line), x = 0.7 (dot-dot-dot-dash line). 
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Figure 6.13: Dependence of the shape of photon and matter power spectra on the mirror 
baryon density u'^ keeping constant Wb. The reference model (soUd hne) has: rim = 0.3, 
= = 0.02, X = 0.7 (not 0.2, as previous figures), h = 0.7 and n = 1.0. For other 
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Figure 6.14: Dependences of the locations of the CMB acoustic peaks on the values of the 
cosmological parameters: Om, wj, = a;{,, h with $7^ = $7^^ = const., h with ujy, = Li>[^ = const., 
n, N^, X, and lo'^ with Wb constant and x = 0.7. The three indicators used here are the 
differences of the positions of the first three peaks of the models from the ones of the reference 
model. The reference model has: Qm = 0.3, ojh = oj{^ = 0.02, x = 0.2, h = 0.7 and n = 1.0. 
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Figure 6.15: Dependences of the temperatures of the CMB acoustic pealcs on the values of the 
cosmological parameters: ^l^a, = <^b! h with Q^, = Q'^ = const., h with lo]^ = = const., 
n, N^, X, and lo'^ with constant and x = 0.7. The three indicators used here are the 
deviations of the differences between the temperatures of the peaks from the same quantities 
obtained for a reference model. The reference model has: Jim = 0.3, = uj[^ = 0.02, x = 0.2, 
h = 0.7 and n = 1.0. 
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6.6 Comparison with observations 

So far we have studied the behaviour of the photon and matter power spectra varying 
many parameters with special attention to the two mirror parameters, i.e. the ratio of 
the temperatures of two sectors x and the amount of mirror baryons a;(,. 

Now we want to compare these models with some experimental data, in order to 
evaluate the compatibility of the mirror scenario with observations and possibly restrict 
the parameter ranges. 

As written in § 6.2, we are not able to fit the parameters now. It is due to the 
slowness of our present version of the numerical code, but our game will be to choose 
some representative model and compare it with observations. 

In the last decade the anisotropics observed in the CMB temperature became the 
most important source of information on the cosmological parameters: a lot of exper- 
iments (ground-based, balloon and satellite) were dedicated to its measurement. At 
the same time, many authors (see, e.g., Percival et al. [140, 141] or Wang et al.[180]) 
proved that its joint analysis with the fluctuations in the matter distribution (they 
have both the same primordial origin) are a powerful instrument to determine the pa- 
rameters of the Universe. As in § 6.3 and § 6.4, we analyse separately the variation of 
X and a;(, in the mirror models, using now both the CMB and LSS informations at the 
same time. 

In order to compare with observations, we use the best available data: for the 
CMB the COBE-DMR [167], MAXIMA [111], BOOMERANG [131], DASI [87, 145] 
and CBI Mosaic [136] data, and for the LSS the IRAS PSCz survey [155] (in particular 
the decorrelated power spectrum provided by Hamilton & Tegmark (2002) [88]). 

We start from figure 6.16, where we plot the models of the upper part of figures 6.1 
and 6.3 together with the observations. This is useful to analyse models with different 
values of x and without CDM (i.e. all the dark matter is made of mirror baryons). In 
top panel, we see that with the accuracy of the current anisotropy measurements the 
CMB power spectra for mirror models are perfectly compatible with data. Indeed, the 
deviations from the standard model are weak, even in a Universe full of mirror baryons 
(see § 6.3). In lower panel, instead, the situation is completely different: oscillations 
due to mirror baryons are too deep to be in agreement with data, and only models 
with low values of x (namely x < 0.3) are acceptable. This is an example of the great 
advantage of a joint analysis of CMB and LSS power spectra, being this conclusion 
impossible looking at the CMB only. Thus, we see the first strong constraint on the 
mirror parameters space: models with high mirror sector temperatures and all the dark 
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matter made of mirror baryons have to be excluded. 

In figure 6.17 we sfiow the models of the bottom panels of figures 6.1 and 6.3, 
compared with the same observations, in order to analyse models with the same x, 
but different mirror baryon contents. The above mentioned low sensitivity of the 
CMB power spectra on doesn't give us indications for this parameter (even for high 
values of x) , but the LSS power spectrum helps us again, confirming a sensitivity to the 
mirror parameters greater than the CMB one. This plot tells us that also high values 
of X can be compatible with observations if we decrease the amount of mirror baryons 
in the Universe. It is a second useful indication: in case of high mirror temperatures 
we have to change the mirror baryon density in order to reproduce the oscillations 
present in the LSS data. 

Therefore, after the comparison with experimental data, we are left with three 
possibilities for the Mirror Universe parameters: 

• high X and low a;(, (differences from the CDM in the CMB, and oscillations in 
the LSS with a depth modulated by the baryon density) ; 

• low X and high a;(, (completely equivalent to the CDM for the CMB, and few 
differences for the LSS in the hnear region) ; 

• low X and low c<j(, (completely equivalent to the CDM for the CMB, and nearly 
equivalent for the LSS in the linear region and beyond, according to the mirror 
baryon density). 

Thus, with the current experimental accuracy, we can exclude only models with high 
X and high a;^. 

Observing the figures, we are tempted to do some (surely premature, given the low 
experimental accuracy and the lack of enough mirror models to fit the parameters) 
guesses. In fact, it's worth noting that in the CMB spectrum the mirror models seem 
to better fit the high first peak, while the mirror baryons could hopefully reproduce the 
oscillations present in the LSS power spectra. We are waiting for more accurate data 
(as expected for example from the MAP satellite) and preparing a faster numerical 
code, in order to fit the new data with a mirror model and establish if there are 
indications for a standard Universe or a Mirror Universe. 

Our next step will be to consider some interesting mirror models and compute their 
power spectra. 

In figure 6.18 we plot models with equal amount of ordinary and mirror baryons 
and a large range of temperatures. This is an interesting situation, because the case 
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= Jib could be favoured in some baryogenesis scenario, as the one proposed by 
Bento and Berezhiani (2001) [14], which considered a lepto-baryogenesis mechanism 
that naturally lead to equal baryon number densities in both visible and hidden sectors. 
These models are even more interesting when we consider both their CMB and LSS 
power spectra. In top panel of figure we see that the temperature anisotropy spectra 
are fully compatible with observations, without large deviations from the standard 
case. In bottom panel, we have a similar situation for the matter power spectra, with 
some oscillations and a shghtly greater slope, that could be useful to better fit the 
oscillations present in the data and to solve the discussed problem of the desired cutoff 
at low scales. Let us note that we are deliberately neglecting the biasing problem, 
given that an indication on its value can come only from a fit of the parameters; so, 
we have in fact a small freedom to vertically shift the curves in order to better fit the 
experimental data. 

Models of a Mirror Universe where the dark matter is composed in equal parts by 
CDM and mirror baryons are plotted in figure 6.19. Now we concentrate on x- values 
lower than the previous figure, because the greater mirror baryonic density would 
generate too many oscillations in the linear region of the matter power spectrum. In 
top panel we show that, apart from little deviations for the model with higher x, 
all other models are practically the same. In bottom panel, instead, deviations are 
big, and we can still use LSS as a test for models. Indeed, models with x > 0.4 are 
probably to exclude, even taking into account a possible bias. Models with lower x are 
all consistent with observations. 

Extremal models 

At the end of this chapter we want to present some so-called "extremal" model. This 
name is due to the characteristic that some of their parameters arc probably outside 
the present estimated limits, but nevertheless they show interesting features. 

Let we consider a flat Mirror Universe where all the energy density is due to matter 
without cosmological constant contribution: Qq — Qm = 1. We can immediately place 
two objections to this scenario. The first is that supernovae measurements [143] seem to 
indicate an accelerated expansion of the Universe, something easily understandable as 
a vacuum energy density effect. The second is related to the constraints on the matter 
density arising from the baryon fraction in galaxy clusters (combined with primordial 
nucleosynthesis data) [55] and from dynamical measurements of fl^a [34]: they both 
predict low matter density values. Against these arguments, we can say that there 
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are still doubts on the reliability of the SNIa measurements; the accelerated expansion 
could be due to some other mechanism (see, e.g., Schwarz (2002) [158]); cosmological 
constant has been always opposed by scientific community and there aren't commonly 
accepted theories on it; the cited estimates on the matter densities could not be valid 
in a Mirror Universe, given the differences between the mirror baryonic and the cold 
dark matter structure formation scenarios. 

Our main point here is to show that a flat Mirror Universe all made of matter could 
be compatible with observations of CMB and LSS power spectra alone, forgetting 
the possible constraints coming from other measurements. To do this we have two 
possible kinds of models: those with low values and those with high x- values. The 
first models have the problem that the value of h is at the limit of the 2a confidence 
interval estimated by the HST Key Project [77] (but well inside other estimations cited 
in § 1.7); in addition, many CMB and LSS joint analyses (included that containing the 
concordance model of Wang ct al. (2002) [180]) obtain low values of h, if they don't 
impose priors on its value. The second ones have the problem that high values of x 
are in contrast with BEN limits {x < 0.64), as discussed in § 4.2. However, not very 
plausible though, one could consider the case with larger x as well, up to x = 1. This 
could be achieved, without contradiction to the BBN limits, if after the BBN epoch 
(but before the matter-radiation equality) there is some additional entropy production 
in the mirror sector (or both in the O and M sectors), due to decay of some light 
metastable particles in both O and M photons. 

Just for demonstration, one can imagine some axion-like boson g with mass Mg ~ 
100 KeV which was in equilibrium sometimes in the early Universe, but has decoupled 
already before electroweak epoch, T ~ 100 GeV. In this case their temperature Tg 
at the nucleosynthesis epoch is related to the temperature of the rest of the thermal 
bath, T ~ 1 MeV, as Tg/T < [g^{l MeV)/5*(100 GeV)^/^ ^ 0.5, where g^{l MeV) = 
10.75 and 5'*(100 GeV) = 106.5 stand for the total number of the particle degrees of 
freedom in the standard model (see appendix A. 5). Therefore their contribution in the 
nucleosynthesis is negligible, namely it corresponds to about AA*",^ ^ 0.03. However, 
when the temperature drops down below 100 keV, the contribution of these particles in 
the energy density of the Universe becomes more relevant, and at some moment it can 
approach the contribution of the photon energy density (because after these particles 
become non relativistic, at Tg < Mg, their energy density rescales as oc MgTg, whereas 
the photon one as oc T^). If at this moment these particles decay into both ordinary 
and mirror photons, clearly with the same rate due to mirror symmetry, this would 
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heat up both photon sectors.^ Therefore, if before the decay of g the temperature of 
the M photons was smaller than the O ones, T' < 0.64 T (as required by the BBN 
hmits) so that energy densities were related as p'^/ < (T'/T)^ Ri 0.17, after decay 
the same amount of energy will be deposited to both sectors, Ap'^ — Ap^. Thus, e.g. if 
the energy density of the ordinary photons is nearly doubled, i.e. Ap-y p^, then the 
O and M photons energy densities after decay should be related as p'^/pj ^ 0.5, which 
corresponds now to their temperature ratio x — T'/T 0.5^/^ 0.85. In other words, 
in the epoch relevant for CMB formation, due to additional entropy production, the 
temperature of mirror sector relative to the ordinary one can be much larger than at 
the nucleosynthesis epoch.^ 

In figure 6.20 we plot all these extremal models comparing them with the concor- 
dance model. For low /i- values, the mirror model is coincident with the concordance 
one for the first two peaks in the temperature anisotropy spectrum and higher for 
other peaks (still in agreement with observations). In the matter power spectrum, 
apart from the change in the slope at higher scales, there are some deviations at lower 
scales (still compatible with experimental data) . Note that the ordinary model for the 
same h and fim values is absolutely incompatible with observations, specially those for 
the LSS. As regards models with high x, they are in good agreement not only with the 
CMB data, but also with the LSS ones, with some possible problem for the latter only 
at high scales. 

We want to stress again that these models have the unique purpose of being spec- 
ulative exercises or curiosities, without any ambition at this stage. 



^One has to remember that this decay should take place before t ^ lO^s, otherwise it would affect 
the thermal Planck spectrum of the CMB photons. 

^Of course, the ratio of the baryon to photon number densities, ij = rib/n-y, would be different 
before and after decay of g, however one has to recall that the value of r] determined by the BBN is 
compatible within a factor of two with the one provided by the CMB analysis. 
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Figure 6.16: CMB and LSS power spectra for various mirror models with different values 
of X, compared with observations and with the concordance model of Wang et al. (2002) 
[180] (solid line) of parameters Qq = I, = 0.34, Qa = 0.66, uy, = Qy,h^ = 0.02, = 0.91, 
h = 0.64. The mirror models have the same parameters as the standard one, but with 
X = 0.3,0.5,0.7 and uj'^ = — Wb (no CDM) for all models. Top panel. Comparison 

of the photon power spectrum with the MAXIMA, BOOMERANG, DASI and CBI Mosaic 
data. Bottom panel. Comparison of the matter power spectrum with the IRAS PSCz data. 
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Figure 6.17: CMB and LSS power spectra for various mirror models with different values of 
mirror baryon density, compared with observations and with the concordance model of Wang 
et al. (2002) [180] (solid line) of parameters Qq = 1, = 0.34, Qa = 0.66, uy, = ^kh^ = 
0.02, Us = 0.91, h = 0.64. The mirror models have the same parameters as the standard 
one, but with x = 0.7 and for uj'^ = lo^, 2u)\^, 3c<Jb, 4LiJb- Top panel. Comparison of the photon 
power spectrum with the MAXIMA, BOOMERANG, DASI and CBI Mosaic data. Bottom 
panel. Comparison of the matter power spectrum with the IRAS PSCz data. 
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Figure 6.18: CMB and LSS power spectra for various mirror models with different values 
of X and equal amounts of ordinary and mirror baryons, compared with observations and 
with a standard reference model (solid line) of parameters = 1, = 0.30, I^a = 0.70, 
Wb = r^b/i^ = 0.02, Us = 1.0, h = 0.70. The mirror models have the same parameters as the 
standard one, but with u'^ = uj^ and x = 0.3,0.5,0.6,0.7. Top panel. Comparison of the 
photon power spectrum with the MAXIMA, BOOMERANG, DASI and CBI Mosaic data. 
Bottom panel. Comparison of the matter power spectrum with the IRAS PSCz data. 
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Figure 6.19: CMB and LSS power spectra for various mirror models with different values 
of X and equal amounts of CDM and mirror baryons, compared with observations and with 
a standard reference model (solid line) of parameters JIq = !> ^^m = 0.30, JIa = 0.70, 
= r^b^i^ = 0.02, Hg = 1.0, h = 0.70. The mirror models have the same parameters as the 
standard one, but with uj'^ = wcdm and x = 0.2, 0.3, 0.4, 0.5. Top panel. Comparison of the 
photon power spectrum with the MAXIMA, BOOMERANG, DASI and CBI Mosaic data. 
Bottom panel. Comparison of the matter power spectrum with the IRAS PSCz data. 
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Figure 6.20: CMB and LSS power spectra for standard and mirror models with = 1 
compared with observations and with the concordance model of Wang et al. (2002) [180] 
(solid line) of parameters Qq = 1, Q^i = 0.34, = 0.66, = ^hh"^ = 0.02, = 0.91, 
h = 0.64. The models can be divided into two groups: those with low /i-values (dotted and 
dashed lines), and those with high x- values (dot-dash and dot-dot-dot-dash). Top panel. 
Comparison of the photon power spectrum with the MAXIMA, BOOMERANG, DASI and 
CBI Mosaic data. Bottom panel. Comparison of the matter power spectrum with the IRAS 
PSCz data. 



Chapter 7 

Mirror stars and other mirror 
astrophysical consequences 



7.1 The mirror astrophysical picture inside the 
galaxy 

Until now we studied the so-called linear structure formation for a Mirror Universe, 
reaching the conviction that mirror matter could exist in the Universe, since it is 
compatible with all the available observations (BBN, CMB, LSS); in addition, we 
obtained some useful constraints on the parameters describing the mirror sector. Thus, 
we are pushed to continue our analysis of the Mirror Universe, and study also other 
astrophysical consequences of the mirror matter. It is well known that some topics (as 
the galaxy and star formations) present many difficulties and obscure points also for a 
standard Universe, where we have many direct observables, so it is even more difficult 
to treat them in a mirror scenario. Nevertheless, we can expose our ideas on various 
arguments, and go into more details where we well know the physics responsible for 
the processes, as for the mirror stellar evolution. 

In previous chapters we described the evolution of a Mirror Universe from the in- 
flation to the structure formation in linear scales, i.e., until the rich galaxy clusters. 
We also argued that for these scales the mirror matter should manifest as the CDM 
for some parameter choices {x < 0.3-0.2). Thus, it is crucial to extend our study to 
smaller scales which already went non-linear, like galaxies, and beyond until the small- 
est astrophysical structures, in order to understand the expected crucial differences 
between the mirror baryonic dark matter (MBDM) and CDM scenarios. 
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The galcLxy 

According to our picture, dark matter in galaxies and clusters is made of a mixture 
of mirror baryonic dark matter and cold dark matter, or can be even constituted 
entirely by the mirror baryons (only for low values of x, given the bounds obtained in 
previous chapter). The presence of mirror matter has obviously consequences on the 
distribution of matter in galactic halos, given the different physics of mirror baryons 
and CDM: the first one is a self-collisional dark matter with exactly the same physical 
laws of our visible sector, while the second one is a collisionless matter component. In 
particular, simulations show that the CDM forms triaxial halos with a density profile 
too clumped towards the center, and overproduces the small substructures within the 
halo. In principle, MBDM could avoid both these problems. In chapter 6 the study of 
the large scale structure demonstrated that the power spectrum for a Mirror Universe 
presents a cutoff at scales which depend on the temperature of the mirror sector and 
its baryon content, so that small substructures have much less power than in CDM 
scenarios. Regarding the halo profiles, the self-collisionality of mirror baryons could 
avoid the high central concentration typical of CDM simulations. Clearly we need 
a numerical study of galaxy formation (based on N-body simulations) to test these 
suggestions, and this will be just one of the next works after this thesis. 

Meanwhile, given the complexity of the physics of galaxy formation (actually this 
process in still to be understood), we can introduce some general consideration on 
it. At a stage during the process of gravitational collapse of the protogalaxy. the 
opacity of the system becomes so high that the gas prefers to fragment into protostars. 
This complex phenomenon lead a part of the protogalactic gas to form the first stars 
(probably very massive, with M ~ lO^-lO^M©). The details of this process are clearly 
dependent on the boundary conditions, which are different for the two sectors (as 
explained in chapter 4): (i) the temperature of the mirror sector is lower than that 
of the ordinary one by a factor x; (ii) the mirror baryonic density is higher or equal 
than the ordinary one; (iii) the baryonic chemical composition is very different in the 
two sectors, since we know (see § 4.5) that mirror sector is a helium-dominated world 
(with a concentration Y' dependent on x according to the formula (4.43), and in any 
case greater than Y). 

Anyway, whatever the details of the scenario dependent on the exact composition 
of matter in the Mirror Universe, one has to take into account the occurrence that 
during the galaxy evolution the bulk of the M baryons could lastly fragment into 
stars. A difficult question to address here is related to the star formation in the M 
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sector, also taking into account that its temperature/density conditions and chemical 
contents are much different from the ordinary ones. In any case, the fast star formation 
would extinguish the mirror gas and thus could avoid the M baryons to form disk 
galaxies as ordinary baryons do. The M protogalaxy, which at certain moment before 
disk formation essentially becomes the coUisionless system of the mirror stars, could 
maintain a typical elliptical structure.^ Certainly, in this consideration also the galaxy 
merging process should be taken into account. As for the O matter, within the dark M 
matter halo it should typically show up as an observable elliptic or spiral galaxy, but 
some anomalous cases can be also possible, like certain types of irregular galaxies or 
even dark galaxies dominantly made out of M baryons. The central part of halo can 
nevertheless contain a large amount of ionized mirror gas and it is not excluded that 
it can have a quasi-spherical form. Even if the stellar formation is very efficient, the 
massive mirror stars in the dense central region fastly evolve (see § 7.3) and explode 
as super novae, leaving behind compact objects like neutron stars or black holes, and 
reproducing the mirror gas and dust. 

We note that, although the hydrogen cross section an is large, it does not necessarily 
implies that the galaxy core will collapse within a dynamical time, since the inner halo 
should be opaque for M particles. They undergo many scatterings and escape from 
the system via diffusion, so the energy drain can be small enough and the instability 
time can substantially exceed the age of the Universe [89] . 

Mirror stars and MACHOs 

In order to understand the details of the process of galaxy formation and evolution, 
it is crucial to study the mirror star formation (beginning and speed) and evolution. 
Stars play an important role: the fraction of baryonic gas involved in their formation 
becomes coUisionless on galactic scales, and supernovae explosions enrich the galaxy 
of processed collisional gas (stellar feedback). By this reason the following sections are 
devoted to the study of the mirror stellar evolution. 

The existence of mirror stars is guaranteed by the existence of ordinary stars: given 
that two sectors have the same microphysics, stars necessarily form in both of them 
[29]. Clearly, being different the boundary conditions, stars have some differences 
(see § 7.3). The fact that dark matter made of mirror baryons has the property of 

"'^In other words, wc speculate on the possibility that the mirror baryons form mainly the elliptical 
galaxies. For a comparison, in the ordinary world the observed bright galaxies are mainly spiral while 
the elliptical galaxies account about 20 % of them. Remarkably, the latter contains old stars, very 
little dust and shows no sign of active star formation. 
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clumping into compact bodies such as mirror stars leads naturally to an explanation 
for the mysterious Massive Astrophysical Compact Halo Objects (or MACHOs). 

These objects are revealed only by their gravitational effects. If we look at a star in 
the sky and at a certain time a compact invisible object passes near our line of sight, it 
acts as a gravitational niicrolens and doubles the image we are seeing. If the resolution 
of the telescope is not large enough to resolve the two star images, we see only an 
enhancement of its brilliance, producing a typical symmetric, achromatic and unique 
light curve, which depends on the mass of the invisible object and on its distance from 
the line of sight (impact parameter). In figure 7.1 we show an explanatory scheme of 
this phenomenon and typical light curves of microlensing events. 



Figure 7.1: Left panel. Scheme of a microlensing in presence of a MACHO. Right panel. 
Typical light curves of a microlensing event as a function of the impact parameter. 

In the galactic halo (provided that it is the elliptical mirror galaxy) the mirror stars 
should be observed as MACHOs in gravitational microlensing [16, 26, 67, 123, 124]. The 
MACHO collaboration [2] has been studying the nature of halo dark matter by using 
the gravitational microlensing technique. This experiment has collected 5.7 years of 
data and provided statistically strong evidence for dark matter in the form of invisible 
star sized objects, which is what you would expect if there was a significant amount 
of mirror matter in our galaxy. The MACHO collaboration has done a maximum 
likelihood analysis which implies a MACHO halo fraction of 20% for a typical halo 
model with a 95% confidence interval of 8% to 50%. Their most likely MACHO mass 
is between 0.15Mq and O.QMq (depending on the halo model), with an average around 
M ~ 0.5 Mq, which is difficult to explain in terms of the brown dwarves with masses 
below the hydrogen ignition limit M < O.IMq or other baryonic objects [78]. These 
observations are consistent with a mirror matter halo because the entire halo would 
not be expected to be in the form of mirror stars. Mirror gas and dust would also 
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be expected because they are a necessary consequence of stellar evolution and should 
therefore significantly populate the halo. Thus, leaving aside the difficult question of 
the initial stellar mass function, one can remark that once the mirror stars could be 
very old and evolve faster than the ordinary ones, it is suggestive to think that most of 
massive ones, with mass above the Chandrasekhar limit Mch — 1.5 Mq have already 
ended up as supernovae, so that only the lighter ones remain as the microlensing 
objects. Perhaps, this is the first observational evidence of the mirror matter? 

Other mirror astrophysical objects 

Mirror globuleir clusters. It is also plausible that in the galactic halo some fraction 
of mirror stars exists in the form of compact substructures like globular or open clusters, 
in the same way as it happens for ordinary stars. In this case, for a significant statistics, 
one could observe interesting time and angular correlations between the microlensing 
events. 

Supernovae and gamma ray bursts. Given that M baryons form the stars, some 
of them must also explode as M supernovae, with a rate in principle predictable after 
the study of the mirror star formation and evolution. Their explosion in our galaxy 
cannot be directly seen by ordinary observer, however it could be observed in terms 
of gravitational waves, and possibly revealed by their next generation detectors. In 
addition, if the M and O neutrinos are mixed [1, 23, 71, 76, 178], it can lead the 
observable neutrino signal, which could be also accompanied by the weak gamma ray 
burst [27, 179]. 

Supermassive black holes. Another tempting issue is whether the M matter itself 
could help in producing big central black holes, with masses ~ 10^ Mq, which are 
thought to be main engines of the active galactic nuclei. 

Substellar scales. If mirror matter exists in our galaxy, then binary systems consist- 
ing of ordinary and mirror matter should also exist. While systems containing approx- 
imately equal amounts of ordinary and mirror matter are unlikely due to e.g. differing 
rates of collapse for ordinary and mirror matter (due to different initial conditions: 
chemical composition, temperature distribution, etc.), systems containing predomi- 
nately ordinary matter with a small amount of mirror matter (and viccvcrsa) should 
exist. Remarkably, there is interesting evidence for the existence of such systems com- 
ing from extra-solar planet astronomy [68, 75]. 
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7.2 The mirror star models 

Given all the above considerations on the importance of the study of mirror stars in 
order to explain the MACHOs, and the galaxy evolution and formation, we now turn 
to the study of the mirror star evolution. 

As we know (see § 4.1), the microphysics of the mirror sector is exactly the same 
as the visible one, the only changes are due to the boundary conditions. This is a very 
favourable condition for the study of M stars, because the necessary knowledge is the 
same than for the O ones, that we know very well. This means that M stars follow 
the same evolutionary stages than visible ones. A very brief review of stellar evolution 
will be given at the beginning of the next section. 

The same physics for both O and M sectors means that the equations governing 
the mirror stellar evolution and the physical ingredients to put inside them (namely 
the equation of state, the opacity tables, and the nuclear reactions) are the same than 
for visible stars. The only change regards the composition of the M star. In fact, 
while the typical helium abundance for O stars is y ~ 0.24, for the M stars we have 
Y' — 0.40-0.80. This interval is obtained considering that its lower hmit is given by 
the primordial helium abundance coming from the mirror Big Bang nucleosynthesis 
studied in § 4.5. In next section we will evaluate its impact on the evolution of M 
stars. 

If we consider a single isolated star^, its evolutionary and structural properties 
depend only on the mass and the chemical composition. In particular, the latter is 
expressed by the abundances of hydrogen (X), helium (Y), and the so-called heavy 
elements or metals (Z), i.e. all the elements heavier than H and He. ^ 

We computed mirror star models using the evolutionary code FRANEC {Frascati 
RAphson Newton Evolutionary Code), a numerical tool to solve the equations of stellar 
structures. As inputs for this code we chose the opacity tables of Alexander & Ferguson 
(1994) [3] for temperatures lower than 10000 K and those obtained in the Livermore 
laboratories by Rogers & Iglesias (1996) [150] for higher temperatures, the equation of 
state of Saumon, Chabrier and Van Horn (1995) [154], and the grey approximation for 
the integration of stellar atmospheres^. These inputs are valid over the entire ranges 

^We are practically neglecting the interactions existing in systems of two or three stars. 

^In this section we use the prime (') to indicate mirror quantities only if they appear together with 
ordinary ones; otherwise wc don't use it, taking in mind that high y-values refer to mirror stars and 
low F- values to the ordinary ones. 

*The "grey atmosphere" approximation assumes local thermodynamical equilibrium with opacity 
independent of frequency. In this case the temperature in the stellar atmosphere is given by 
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of temperatures and densities reached by our models. 

We had no need to modify the code because, as explained before, mirror stars are 
evolutionary equivalent to ordinary stars with a high helium abundance. Thus, we 
computed stellar models for large ranges of masses and helium contents, and for a low 
metallicity Z. A low Z-value means that we are treating a so-called stellar population 
II, i.e. an old stellar population, coming soon after the first one (the population III, 
without metals). 

We summarize the values used for stellar parameters in table 7.1, noting that the 
model with M — O.SMq has been computed also for the not listed intermediate values 
Y — 0.30, 0.40, 0.60, 0.80, in order to better investigate the dependence of evolutionary 
properties on the helium content for a mass interesting as MACHO candidate. 



parameter 


values 


Z 


10-^ 


Y 


0.24 - 0.50 - 0.70 


M/Mq 


0.5 - 0.6 - 0.8 - 1.0 - 1.5 - 2.0 - 3.0 - 4.0 - 5.0 - 7.0 - 10 



Table 7.1: Parameters and their values for mirror star models. The models are for all the 
combinations of the parameter values. In addition the model with M = O.SMq has been 
computed also for the not listed intermediate values Y = 0.30, 0.40, 0.60, 0.80. 

We note that those which we are now calling mirror stars can perfectly be ordinary 
stars in an advanced evolutionary stage of the Universe, when the galactic matter is 
enriched by many elements processed by the previous stellar populations, thus forming 
stars with more helium in their structures. 

7.3 Evolution of the mirror stars 

First of all we remember that in stellar astrophysics the evolution of a star is studied 
in the so-called H-R (Hertzsprung-Russell) diagram, where we plot the luminosity L 

where T(t) is the temperature of an atmospheric layer located at the optical depth r, and Tf. is the 
effective temperature of the star. 
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and the effective temperature Tg of the star^. In order to understand the evolutionary 
differences between ordinary and mirror stars, it is also necessary to give a very brief 
review of basic stellar evolutionary theory. 

During the first fast phase of gravitational contraction at nearly constant effective 
temperature and decreasing luminosity, the star goes down along his Hayashi track^, 
negligibly slowing down its contraction only while the structure is fastly burning the few 
light elements (D, Li, Be, B) present. Contraction increases the central temperature 
Tc, until stars with masses M ^ 0.1 Mq ignite the hydrogen burning^ in their cores at 
a temperature Tc ~ 6 x 10^ K, while the ones with lower masses do not ignite hydrogen 
and die as brown dwarfs. After depletion of hydrogen in the core, the burning passes 
to a shell at the boundary of the core, which is now made of He. At this stage the star 
starts decreasing its effective temperature (turn-oSt). Meanwhile the He core contracts 
until masses greater than ~ 0.5 Mq reach a central temperature Tc ~ 10^ K and 
start He-burning in C and O; stars with lower masses die as white dwarfs and start 
their coohng sequences. The nucleosynthesis process continues burning elements into 
heavier and heavier nuclei, but only the heaviest stars (M > 6 — 8 Mq) complete 
the advanced evolutionary stages and die exploding as type 11 supernovae and leaving 
in their place a neutron star or a black hole. An evolutionary feature common to all 
burning phases is that shell burnings induce convective phenomena which push the star 
toward his Hayashi track, namely toward lower effective temperatures, while central 
burnings push it toward higher Tg. 

A key point is the evaluation of evolutionary times. The most lasting phase of hfe of 
a star is surely that corresponding to the central hydrogen burning (the so-called main 
sequence), with time-scales of 10^° yr for masses near the solar mass. Thus, we can 
approximate the lifetime of a star with its main sequence time. Since both luminosity 

^The effective temperature of a star is defined by 

L = AnR^aT^ , 

where a is the Stephan-Boltzmann constant, L is the luminosity, and R is the radius at the height of 
the photosphere. Thus, Tg is the characteristic temperature of the stellar surface if it emits as a black 
body. 

^Hayashi track is the evolutionary track of a totally convective stellar model. It is the coldest 
possible track for a star of a given mass, and it is located at the extreme right of the H-R diagram. 

^Thc exact vahic of the hydrogen burning minimum mass Adhbmm is dependent on the metallicity. 
For our models Z = 10"^ and Mhbmm — 0.1 Mq, while for solar metallicity Z = 0.02 and Mh^mm — 
0.08 Mq. 

^Thcrc arc two possible ways of burning hydrogen: the first one, called PP or proton-proton chain, 
becomes efficient at Tc ~ 6 x 10® K, while the second one, called CNO chain, at Tg ~ 15 x 10® K. 
Since their efficiencies are dominant at different temperatures, the PP chain provides energy for 
smaller stellar masses, and the CNO does it for bigger ones. 



§7.3 Evolution of the mirror stars 



155 



and effective temperature depend on the mass and chemical composition, clearly the 
lifetime too depends on them. We use now the proportionality relations [46] 

L OC i^^-^M^-^ (7.1) 

and 

« 1^'-' , (7.2) 

where fj, is the mean molecular weight. Prom eq. (7.1) we obtain that bigger masses 
need higher luminosities, so that they use all the available hydrogen earlier than the 
hghter ones. Prom both eqs. (7.1) and (7.2) we know that an increase of helium 
abundance corresponds to an increase of the mean molecular weight and consequently 
in both luminosity and effective temperature. The increase in luminosity means that 
the star needs more fuel to produce it, but at the same time its amount is lower, because 
higher y-values necessarily imply lower X-values. Both these events act to shorten 
the lifetime of a mirror star, which has a high He content. This can be formalized in 
the following relation [46] 

tMS OC ^ , (7.3) 

where Ims is the main sequence lifetime. 

After these predictions on evolutionary properties of He-rich stars, we analyse the 
quantitative results of our models. They can be divided into two groups. The first one 
is made of models with mass M = 0.8Mq and many different F-values. The second 
one is made instead of models with only three different He contents and a large range 
of masses. 

We start from figure 7.2, where we plot the models of O.SMq in the H-R dia- 
gram. The models arc followed until the He-burning ignition, i.e. along their main 
sequence, turn-off and red giant^ phases, which practically occupy all their lifetimes. 
Our qualitative predictions are indeed confirmed. Models with more helium are more 
luminous and hot; for example the main sequence luminosity ratio of the model with 
Y = 0.80 to the one with Y = 0.24 is ~ 10^. Other consequences of an He increase 
are a longer (in the diagram, not in time) phase of decreasing temperature at nearly 
constant luminosity, and a shorter red giant branch. 

From these models we computed the evolutionary times until the He-ignition, i.e. 
for the entire plotted tracks, and we summarize them in table 7.2. As expected, the 
ages decrease for growing Y, but we see now how much high is this correlation. Por 

red giant is a cold giant star in the phase of H-burning in shell before the He-ignition. 
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Y — 0.40 the lifetime is already about one third compared to a visible {Y — 0.24) 
star, while for the highest value, Y — 0.80, it is roughly 10^ times lower. We can 
approximately say that an increase of 10% in helium abundance roughly divides by 
two the stellar lifetime. In figure 7.3 we plot the evolutionary times listed in the table. 
We see that, using a logarithmic scale for the stellar age, we obtain a quasi-hnear 
relation between it and the helium content for this range of parameters. 



o 




4.2 4.0 3.8 3.6 

log Te 



Figure 7.2: Evolutionary tracks in the H-R diagram of stars with M = O.SMq, Z = 10 ^ 
and different helium contents Y = 0.24, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80. 

Let us now extend the analysis to models covering a large range of masses, from 0.5 
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Y 


0.24 


0.30 


0.40 


0.50 


0.60 


0.70 


0.80 


age (lO^yr) 


13.2 


8.53 


4.50 


2.17 


1.01 


0.417 


0.169 



Table 7.2: Ages computed for stars of mass M = O.SMq and the indicated helium contents. 
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Figure 7.3: Evolutionary times (listed in tabic 7.2) for models with M = 0.8 Mq, Z = 10 
and different helium contents Y = 0.24, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80. 

Mq to 10 Mq. Since the dependence on the helium content has been already studied 
for the 0.8 Mq case, we concentrate on only three F-values. Figure 7.4 shows the 
evolutionary tracks for all the masses and only two helium contents, and again the 
models are followed up to the He-ignition. For every mass the F-dependence is the 
same as for the above discussed 0.8 Mq model. For models with masses M > 2 Mq the 
growth in Y causes a considerable increase of the He-ignition effective temperature, 
together with the disappearance of the red giant branch. 

In table 7.3 we list the lifetimes for all masses and the three indicated helium 
contents. The ratio of an ordinary star (Y = 0.24) evolutionary time to the high 
He-content mirror one {Y = 0.70) is between ~ 30 for 0.5 Mq and ~ 10 for 10 Mq. 
These data are plotted in figure 7.5, where we see that the same dependence on the 
star masses holds for every helium content, with a shift toward lower ages for higher 
y-values. 

This is an evidence that, under large mass ranges and different boundary conditions 
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(in terms of temperatures of the mirror sector, and thus stellar helium content), the 
lifetimes of mirror stars are roughly an order of magnitude greater than the ones of 
visible stars. 

This means that, compared to O stars, M stars evolve faster and enrich earher 
the galaxy of processed mirror gas, with implications for galaxy evolution. Prom the 
detailed study of this evolution together with the necessary information of the initial 
mirror stellar mass function, we could predict the expected population of mirror stars, 
in order to compare it with current MACHO observations. In addition, we could 
evaluate the amount of gravitational waves expected from mirror supernovae. These 
are just some interesting future applications of the present study. 



mass 


age (yr) 


age (yr) 


age (yr) 


{M/M. ) 


{Y = 0.21) 


{Y = 0.50) 


{Y = 0.70) 


0.5 


7.06 X 10^° 


1.12 X 10^0 


1.92 X 10^ 


0.6 


3.71 X 10^° 


5.80 X 10^ 


1.04 X 10^ 


0.8 


1.32 X 10^° 


2.17 X 10*^ 


4.17 X 10^ 


1.0 


6.07 X 10^ 


1.05 X 10^ 


2.19 X 10^ 


1.5 


1.60 X 10^ 


3.01 X 10^ 


9.02 X 10^ 


2.0 


6.43 X 10^ 


1.40 X 10^ 


4.45 X 10^ 


3.0 


2.14 X 10« 


5.64 X 10'^ 


1.78 X 10^ 


4.0 


1.15 X 10^ 


3.13 X 10^ 


9.41 X 10^ 


5.0 


7.26 X 10^ 


2.05 X 10^ 


6.56 X 10^ 


7.0 


3.81 X 10^ 


1.17 X 10^ 


4.14 X 10^ 


10 


2.09 X 10^ 


7.18 X 10^ 


2.78 X 10^ 



Table 7.3: Ages computed for stars of the indicated mass and helium content, with a metal- 
licity Z = 10-^. 
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Figure 7.4: Evolutionary tracks in the H-R diagram of stars with different masses M = 
0.5,0.6,0.8,1.0,1.5,2.0,3.0,4.0,5.0,7.0, 10 Mq, Z = 10"^ and two different hehum contents 
y = 0.24 (red hne) and 0.50 (blue). 
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Figure 7.5: Evolutionary times (listed in table 7.3) for models with M = 0.5, 0.6, 0.8, 1.0, 
1.5, 2.0, 3.0, 4.0, 5.0, 7.0, 10 M©, Z = 10"^ and three different helium contents Y = 0.24 
(red dots), 0.50 (blue) and 0.70 (green). Are also indicated the models of figure 7.3 (black). 



Chapter 8 
Conclusions 



The aim of this thesis is contained in a question: "is mirror matter a reUable dark 
matter candidate?" Its emergence arises from the problems encountered by the stan- 
dard candidate, which is now the cold dark matter, in some aspects, as for example the 
central galactic density profiles or the number of small satellites. In order to answer 
this question, we studied the cosmological implications of the parallel mirror world 
with the same microphysics as the ordinary one, but which couples the latter only 
gravitationally, and its consistence with present observational data, in particular the 
ones coming from the Big Bang nucleosynthesis, the cosmic microwave background 
radiation and the large scale structure. 

The nucleosynthesis bounds on the effective number of extra light neutrinos demand 
that the mirror sector should have a smaller temperature than the ordinary one, T' < T, 
with the limit x ~ T'/T < 0.64 set by the constraint AN^, < 1. By this reason its 
evolution should be substantially deviated from the standard cosmology as far as the 
crucial epochs like baryogenesis, nucleosynthesis, baryon-photon decoupling, etc. are 
concerned. 

Starting from an asymmetric inflationary scenario which could explain the different 
initial temperatures of the two sectors, in the context of both the GUT or electroweak 
baryogenesis scenarios the condition T' < T yields that the mirror sector should pro- 
duce a larger baryon asymmetry than the observable one, rj'^ > r)B- 

Therefore, the temperature bound implies that the mirror sector contains less rel- 
ativistic matter (photons and neutrinos) than the ordinary one, Q'^. <^ Qr, so that in 
the relativistic expansion epoch the cosmological energy density is dominated by the 
ordinary component, while the mirror one gives a negligible contribution. However, 
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for the non-relativistic epoch the complementary situation can occur when the mirror 
baryon matter density is bigger than the ordinary one, ft'^ > fts- Hence, the mirror 
baryonic dark matter can contribute the dark matter of the Universe along with the 
cold dark matter or even constitute a dominant dark matter component. We know 
also that mirror world must be a helium-dominated world, since in its sector the Big 
Bang nucleosynthesis epoch proceeds differently from the ordinary one, and it predicts 
the mirror helium abundance in the range Y' — 0.5-0.8, considerably larger than the 
one for observable helium, Y ~ 0.24. 

Since the existence of a mirror hidden sector changes the time of key epochs, there 
are important consequences in the structure formation scenario for a Mirror Universe. 
We studied this scenario in presence of adiabatic scalar density perturbations, which 
are now the most probable kind of primordial fluctuations, in the context of the Jeans 
gravitational instability theory. 

Given that the physics is the same in both sectors, key differences are the shifts of 
fundamental epochs, namely the matter-radiation equality occurs in a Mirror Universe 
before than in a standard one, (aeq)mir < (aeq)ord, and the baryon-photon equipartition 
and the matter-radiation decoupling do the same: a(,^ < ab-y; ajec < ^dec- The first step 
is given by the study of the mirror sound speed and its comparison with the ordinary 
one and with the velocity dispersion of a typical cold dark matter candidate. Prom 
this study we obtain the mirror Jeans length and mass, again to be compared with 
the same quantities obtained for the ordinary sector and for the cold dark matter. 
There are two different possibilities, according to the value of x, which can be higher 
or lower than Xeq ~ 0.046(Qin^^)~^, the value for which mirror decoupling occurs at 
matter-radiation equality time. 

The values of the length and mass scales clearly depend on the mirror sector tem- 
perature and baryonic density, but we found that M'j is always smaller than Mj, with 
a typical ratio ~ 10 for x > Xeq, while for cold dark matter it is ~ 10^^. 

Another important quantity to describe the structure evolution is the dissipative 
scale, represented in the mirror sector by the mirror Silk scale. We found that it is 
much lower than the ordinary one, obtaining Mg ~ 10^° M© for x ~ a^eq, a value similar 
to the free streaming scale for a typical warm dark matter candidate, and much higher 
than the one for cold dark matter. 

We put together all these informations to build two different mirror scenarios, for 
X > Xeq and x < Xeq. In the latter case we obtain a maximum mirror Jeans scale 
similar to the Silk mass, so that practically all perturbations with masses greater than 
the Silk mass grow uninterruptedly, just as in a cold dark matter scenario. 
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After this, we modified a numerical code existing for the standard Universe in order 
to take into account a hidden mirror sector, and computed the evolution of perturba- 
tions in the linear regime for all the components present in a Mirror Universe, namely 
the ordinary and mirror baryons and photons, and the cold dark matter. We did 
this for various mirror temperatures and baryon densities, and for different perturba- 
tive scales, finding all the features predicted by our structure formation study (as for 
example the mirror decoupling and the CDM-like behaviour for low x- values). 

Using the same numerical code we were able to predict the expected power spectra 
of cosmic microwave background and large scale structure for a fiat Mirror Universe. 
We analysed the dependence of power spectra on mirror parameters x and and 
also on other cosmological parameters, as Qm, h, n, N^. 

In CMB spectra we found various differences from a so-called standard concordance 
model for x > 0.3 and a dependence not linear in x, specially evident in the first and 
third peaks. The dependence on the mirror baryon density is instead very low. We 
computed also the power spectrum of mirror CMB photons, even if this study is now 
only an academic exercise, because by definition we can't reveal them. 

Turning to the LSS power spectra, wc showed an influence of the mirror sector 
bigger than for the CMB, with a great dependence on both mirror temperature and 
baryonic density. Both of them cause oscillations in power spectrum, but the flrst one 
influences the scale at which they start, while the second one their depth. In this case 
we see the mirror sector effects also for low a;- values, if we don't take a too small value 
for Vt'^. 

Wc extended the models also to smaller (non linear) scales, in order to show the 
cutoff present in the mirror scenario. We demonstrated the existence of this cutoff, 
specially dependent on value, but modulated also by Q'^. This is an important feature 
of the mirror structure formation scenario, because it could explain the observed small 
number of satellites which is a problem for cold dark matter, and at the same time 
this is valid also for low x-values, that give mirror baryons equivalent to CDM for the 
CMB and LSS at linear scales. 

The next step was the comparison with observations. In this phase it is important 
the joint analysis of both CMB and LSS data, which gave us an important limit on the 
mirror parameter space. In fact, wc obtained the conclusion that mirror models with 
high X and high fl'^ are excluded by LSS observations, because they generate too deep 
oscillations in power spectra. This is an important bound, which limits to the three 
following possibilities to have a mirror sector: 
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• high X and low Q'y^ (differences from the CDM in the CMB, and oscillations in 
the LSS with a depth modulated by the baryon density) ; 

• low X and high Q'^ (completely equivalent to the CDM for the CMB, and few 
differences for the LSS in the linear region); 

• low X and low Q{j (completely equivalent to the CDM for the CMB, and nearly 
equivalent for the LSS in the linear region and beyond, according to the mirror 
baryon density). 

Thus, with the current experimental accuracy, we can exclude only models with high 
X and high fl'^, but with the soon available high precision data on CMB (Map, Planck) 
and LSS (2dF and SDSS) we will be able to choose between the CDM and mirror 
cosmological scenarios. 

In last chapter we qualitatively discussed the implications of the mirror baryons 
representing a kind of self-interacting dark matter for the galactic halo structure and 
MACHOs. Both of them can be explained as mirror consequences, the first one as 
mirror galaxy, the second ones as mirror stars. Finally, we computed the evolutionary 
properties of mirror stars using a numerical code. Stars in mirror sector are different 
only for their helium content {Y' = 0.4-0.8) from the visible ones (Y ~ 0.24). This 
difference implies much faster evolutionary times, dependent on the exact He content 
and thus on the temperature of the mirror sector. 

Therefore, at the end of this work of thesis, we reached a partial answer: on the hght 
of current observations of BBN, CMB and LSS in linear scales the mirror baryonic dark 
matter is not only fully in agreement with observations, but in some case it could be 
even preferable to the CDM scenario. In addition, we obtained some useful constraints 
on the parameter space, which can address our future efforts to understand other 
aspects of the Mirror Universe. 

Suggestions for future works 

Let us conclude by briefly describing some of the planned future developments of this 
work. They are all linked each other and help to complete the picture of a Mirror 
Universe. 

Fit of all cosmological pcirameters. Since the program used in this thesis for 
computing CMB and LSS power spectra privileges the precision to the performance, 
it is too slow to make us able to build a parameter grid thin enough to fit the ordinary 
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and mirror cosmological parameters. Thus, we will adapt a numerical code faster than 
the one now used (as for example the popular CMBFAST), in order to drastically 
decrease the computational time. 

Mirror stair formation. An important question that still needs to be resolved is the 
mirror star formation. The different primordial chemical content of the mirror sector 
causes differences in opacities and then in fragmentation of protostellar clouds. The 
initial mirror stellar mass function, together with the mirror evolutionary data (that 
we obtained in this thesis), are the necessary inputs to insert as stellar feedback in 
N-body simulations of galaxy formation and evolution. 

Estimation of present mirror star population. With the knowledge of mirror star 
formation and evolution we could in principle predict the present stellar population of 
mirror stars, in order to compare it with current MACHO observations. 
Gravitational waves from the mirror sector. After the study of the present 
mirror star population, we can evaluate the spectrum of gravitational waves expected 
from the mirror sector as produced by supernovae explosions or binary systems of 
compact objects. 

Extension to smaller (non linectr) scales and galctxy formation. Using the 
transfer functions obtained here in the linear approximation, we can use N-body sim- 
ulations to include non-hnear effects, as merging and stellar feedback. In this way we 
can also simulate the galaxy formation in presence of a mirror sector, compute the den- 
sity profiles, and verify if, as we expect, the mirror scenario solves the open problems 
placed by the cold dark matter scenario. 

Supermassive black hole formation. Another important issue worthy of explo- 
ration is the influence of the mirror sector on supermassive black holes formation in 
galaxy centers. Maybe mirror matter could help in producing them. 
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Appendix A 

Useful formula for cosmology and 
thermodynamics 



A.l General Relativity 

The essence of Einstein's theory is to transform gravitation from being a force to being 
a property of space-time, which may be curved. The interval between two events 



(A.l) 



is fixed by the metric tensor g^p which describes the space-time geometry^ {q'^'^QhP 
S^). For the Riemannian spaces, the tensor of curvature is 



or 



up 



-pfM -pa 
J- ^flJ- 



dx'^ dxl^ '"^ '^^ ' 

where the F's are ChristoSel symbols 

1 



(A.2) 
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dxf^ dx"^ dx'^ 



(A.3) 



The equation of motion of a free particle is determined by the space-time metric 

d'^x"' dx^ dx'^ 

so that free particle moves on a geodesic. 



(A.4) 



^Repeated indices imply summation and a, (3 run from to 3; = f is the time coordinate and 
{i = 1, 2, 3) are space coordinates. 
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On the other hand, the metric Qap is itself determined by the distribution of matter, 
described by the energy-momentum tensor Tap, according to the Einstein equations 



Ga/3 = Rap — -^Rgap 



aP 1 



where 



RajB — R. 



a-y/3 1 



R — g^'^Ro.p 



(A.5) 



(A.6) 



are respectively the Ricci tensor and the Ricci scalar, and A is the cosmological con- 
stant. For the FRW metric, the non zero components of the Ricci tensor and the value 
of the Ricci scalar are 

^2 2k' 
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For a perfect fluid the energy-momentum tensor has the form 

Tap = (P + p)UaUp - pgap , 



(A.7) 



(A.8) 



where p and p are respectively the energy density and pressure of the fluid, Ua — 
Qap dx^/ ds is the fluid four- velocity. Considering the symmetries of the FRW metric 
(uniformity and isotropy), which demand that w° = 1 and — in the comoving 
coordinate system, we obtain T^p — diag{p, —p, —p, —p). This is valid also in presence 
of the cosmological constant, if we substitute p and p as indicated in the following 

A 



p -^p-pA 



p ^ P + Pa; 



PA 



-PA = 



SttG 



(A.9) 



Therefore, the Einstein equations for a Universe described by the FRW metric are 
reduced to equations (1-13) and (1.15), where the relations between the energy and 
pressure densities are in general related as p = wp. In particular, for the dominance of 
relativistic and non relativistic matter we have respectively w — 1/3 and w — 0, while 
for the vacuum energy dominance one has w — —1. 



A. 2 Flat models 

Here we report two special cases of the relationships (1.28)-(1.33): dust or matter 
dominated Universe {w — 0) 

/ t \ 

a{t) = ao(^-j (A.IO) 
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t = to(l + z)-^/'' (A.ll) 

2 
M 
1 
2 



= ^^=Hoil + zf' (A.12) 
90 = ^ (A.13) 



= (A.14) 



and radiation dominated Universe {w — 1/3) 

t \ ^2 



a(t) = aoi^-J (A.16) 
i = to{l + z)-^ (A.17) 
= l = i/o(l + ^)' (A.18) 



90 = 1 (A.19) 
1 

3 

32^t2 



^0 = (A.20) 
Pr = ^7t4t^ (A.21) 



A. 3 The synchronous gauge 

Since our interests lie in the physics in an expanding Universe, we use comoving co- 
ordinates = (t, x), with the expansion factor a(T) of the Universe factored out. 
The comoving coordinates are related to the proper time and positions t and r by 
dx^ — dr — dt/a{T), dx — dr/a{T). Dots will denote derivatives with respect to r: 
d = da/ dr. 

The components goo and qqi of the metric tensor in the synchronous gauge are by 
definition unperturbed. The line element is given by 

ds^ = a^{r){-dT'^ + {5ij + hij)dx'dx^} . (A.22) 

The metric perturbation hij can be decomposed into a trace part h = ha and a tracclcss 
part consisting of three pieces, hjj , /i^ , and hfj , where hij — hSij/3 + hjj + hf, + /i^- ; h}lj 
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can be written in terms of some scalar field /i and hjj in terms of some divergenceless 
vector A as 

hij = diAj + djAi, diA^O. (A.23) 

The two scalar fields h and n (or hjj) characterize the scalar mode of the metric 
perturbations, while A^ (or hjj) and hfj represent the vector and the tensor modes, 
respectively. 

We will be working in the Fourier space k. We introduce two fields /i(k, r) and 
77 (k, r) in /c-space and write the scalar mode of hij as a Fourier integral 

/ii^(x, r)^ J d^/te*-^ |^i%/i(k, r) + (4% - ^Sij) 6r){k, r)| , k ^ kk . (A.24) 

Note that h is used to denote the trace of hij in both the real space and the Fourier 
space. 

In spite of its wide-spread use, there are disadvantages associated with the syn- 
chronous gauge: spurious gauge modes contained in the solutions to the equations 
for the density perturbations, due to the arbitrary choice of the initial hypersurface 
and its coordinate assignments; coordinate singularities arising when two observers' 
trajectories intersect each other (a point in spacetime will have two coordinate labels, 
since the coordinates are defined by freely faUing observers) , so that a different initial 
hypersurface of constant time has to be chosen to remove these singularities. 

Einstein equations and energy-momentum conservation 

We find it most convenient to solve the linearized Einstein equations in the Fourier 
space k, and a cosmological constant is allowed through its inclusion in p and p using the 
(A. 9). In the synchronous gauge, the scalar perturbations are characterized by /i(k, r) 
and 77(k, r) in equation (A.24). In terms of h and 77, the time-time, longitudinal time- 
space, trace space-space, and longitudinal traceless space-space parts of the Einstein 
equations give the following four equations to linear order in /c-space: 

inGa^ST^ , (A.25) 

A7rGa^(p + p)e , (A.26) 

-87rGa^ST\ , (A.27) 

-24:7iGa^{p + p)a . (A.28) 



k^T] -h = 

h + 2-h-2k'^r] = 
a 

h + 67i + 2-(h + 6f)) - 2k'^rj = 
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The variables 9 and a are defined as 

{p + p)e = ik'6T°j , (p + p)a = -{hkj - ^Sij)E'j , (A.29) 

and Sj- = T j — denotes the traceless component of T j. When the different 

components of matter and radiation (i.e., CDM, HDM, baryons, photons, and massless 
neutrinos) are treated separately, {p+p)0 = Y,i{pi + Pi)9i and (p+p)cr = Y,i{pi+Pi)(^i , 
where the index i runs over the particle species. 

Now we derive the transformation relating ST"p, making use of the definition of 
energy-momentum tensor (A. 8). For a fluid moving with a small coordinate velocity 
V* = dx^dr, can be treated as a perturbation of the same order as Sp — p — p, 
Sp = p — p, and the metric perturbations. Then to linear order in the perturbations 
the energy-momentum tensor is given by 

n = -ip+sp), 

T'j = (p + 5p)5'j + E^,- , S\ = , (A.30) 

where we have allowed an anisotropic shear perturbation S in T y 

The conservation of energy-momentum is a consequence of the Einstein equations. 
Let w = p/p describe the equation of state (see §1.1). Then the perturbed part of 
energy-momentum conservation equations 

r"^;« = dj:"^'^ + rVT^'^ + r^r^'' = o (A.31) 

in A;-space implies 

e = -^-^l-2,w)9-^e+^-^k-5-k^a. (A.32) 

' 1 + w l + w ^ ' 

These equations are valid for a single uncoupled fluid, or for the net (mass-averaged) 5 
and 9 for all fluids. They need to be modified for individual components if they interact 
with each other. An example is the baryonic fluid in our model, which couples to the 
photons before recombination via Thomson scattering. An extra term representing 
momentum transfer between the two components needs to be added to the 5 equation 
for the baryons. 
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For the isentropic primordial perturbations considered in this thesis, the equations 
above simphfy since Sp — c'^Sp, where = dp/ dp — w + pdw / dp is the adiabatic sound 
speed squared. For the photons and baryons (the only coUisional fluid components 
with pressure), w is a constant = 1/3 for photons and w ~ for baryons, since 
they are non relativistic at the times of interest). Thus, Sp/5p — w — 0. 

A. 4 Thermodynamics of the Universe 

In a gas of a given species with g internal degrees of freedom and energy E — 



ylpj^ + m^, kinetic equilibrium is established by sufficiently rapid elastic scattering 
processes; in this case, for an ideal gas, the equilibrium phase-space density is 



where +/— refers to Fermi-Dirac/Bose- Einstein statistics and p, is the chemical po- 
tential. In general each species has its own equilibrium temperature T, and the entire 
Universe can be represented as a plasma with different temperatures. However, if sev- 
eral species strongly interact among them, they will reach a mutual equilibrium and 
a common temperature; this is indeed the situation at early times. As the Universe 
expands and cools down, some species may start interacting more and more weakly 
and eventually decouple. We can consider the photon temperature as the plasma 
reference temperature T of the Universe. 

In chemical equilibrium, established by processes which can create and destroy 
particles (differently from kinetic equilibrium), the chemical potential is additively 
conserved. So it is zero for particles such as photons and Z'^, which can be emitted and 
absorbed in any number, and consequently opposite for a particle and its antiparticle, 
which can annihilate into such bosons. 

The quantities of interest are the number density n, energy density p and pressure 
p of a given species, defined for a generic species of mass m with chemical potential p 
at temperature T as 





-1 



(A.33) 



n 




(A.34) 
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(A.35) 
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67r2 Jm exp [{E - /x) /T] ± 1 



(i£ . 



(A.36) 



For non relativistic species {T <^ m) we have (for both Fermi-Dirac and Bose-Einstein 
statistics) 



n ^ 9[—) exp ^ 



(A.37) 

(A.38) 
(A.39) 



p n + 2 ^ nm ^ 
p ~ nT -C p . 

The average energy per particle (E) = p/n is instead given by 

(E) + -T . (A.40) 

For relativistic species (T ^ m) we obtain in the non degenerate case (T ^ //) 



p ~ < 



f (3/4) my)9T' 

[ (C(3)M 

f (7/8) (7rV30) 

[ (7rV30)^r^ 

P - (1/3) P 

'(7/6) (7r730) C(3)T ~ 3.15T FD 

(ttVSO) C(3)r ~ 2.70 T 



SF 



(A.41) 

(A.42) 
(A.43) 
(A.44) 



In the relations above ({x) is the Ricmann zcta-function, and ({3) ~ 1.2021. 

Another important quantity for the evolution of a thermodynamic system is its 
entropy S. Assuming zero chemical potentials^ the entropy per comoving volume is 
defined, up to an additive constant, by ^ 



V 



P + P 



.3 P + P 

' T 



(A.45) 



Note that throughout most of the history of the Universe (in particular the early 
Universe) the reaction rates of particles in the thermal bath, F, were much greater 



^This is a very good approximation, as all evidence indicates that <C T. 
^With a chemical potential we find S = a^{p + p — fin)/T. 
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than the expansion rate, H, and thermal equihbrium^ should have been maintained. 
In this case the entropy in a comoving volume is conserved 

dS = 0. (A.46) 

It is also useful to define the entropy density in the comoving volume, 

s = - = oc a ^ (A.47) 

Note that if 5 = 4pr/3Tr is the radiation entropy, then S = S/ksri}) is the entropy 
per baryon, where is the radiation temperature, /cb is the Boltzmann constant and 

is the baryon number density. Given that Pr oc and pb oc riy,, the entropy per 
baryon satisfies the relation 

S oc . (A.48) 



A. 5 Particle numbers at key epochs 

From decoupling till present days, neutrinos remain relativistic and therefore continue 
to retain their equilibrium distribution; hence the degrees of freedom characterizing 

the present energy density and entropy (7 + 3 z/) are (using eqs. (1-45) and (1-47)) 

4 

5° = 2 + -(3x2) ^ 3.36, (A.49) 

^ 2+^- (3x2) (^^ 3.91 . (A.50) 
At BBN epoch T ~ 1 MeV and the present particles (7 + e+e~ + 3 u) give 

= = 2 + ^ • 4 + ^ • (3x2) = 10.75 . (A.51) 

o o 

At T > 100 GeV all SM particles - 8 gluons, W^, Z and 7, 3 families of quarks 
(3 colors per each) and charged leptons, 3 families of v and the Higgs doublet - are 
relativistic, and we have 

= gf^ = 106.75 . (A.52) 
In a minimal supersymmetric SM (MSSM) this number should be nearly doubled 

SSM ^ ^MSSM ^ 221.5 . (A.53) 

At T ^ 100 GeV there could be some other particle in the thermal bath, so that 
g* ^ 5**^ or g^ > g^^^^ in the supersymmetric case. 



This is not true during phase transitions, when entropy is not conserved, in general. 



Appendix B 
Structure formation 



B.l Linear Newtonian theory 

The general fluid equations. We model the Universe as a fluid, so that all the 
relevant quantities are described by smoothly varying functions of position. After the 
first bound structures, such as galaxies, form, they are treated like particles along with 
the genuine particles that remain unbound. 

For an expanding fluid it is convenient to adopt a comoving reference frame, speci- 
fied by space coordinates and a universal time t] if we consider a fluid with density 
p and pressure p, moving with "peculiar" velocity Ua in a gravitational potential its 
evolution is governed by the standard Eulerian equations for a self gravitating medium 

^ + ?>Hp+h^{fm'') = 0, (B.l) 
d^a du 1 11 

-r-:,Xo, + ^ + Hua + -u^dpu^+ —dc,p+ -da^ = 0, (B.2) 
di^ at a ap a 

d'^^-A-KGa^p = 0. (B.3) 

Expressions (B.1)-(B.3) are respectively known as the continuity, the Euler and the 
Poisson equations. They describe mass conservation, momentum conservation and the 
Newtonian gravitational potential. 

The unperturbed background. The simplest non-static solution to the system 
(B.1)-(B.3) describes a smoothly expanding = 0), homogeneous and isotropic fluid 
(i.e. po — po{t) , po — Po{t)). In particular, the unperturbed background Universe is 
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characterized by the system 



with solutions 



^ + ^Hpo = 0, (B.4) 

^Xa+-da^O = 0, (B.5) 

a^$o - 471(70% = 0, (B.6) 



2 

Po oc , Vq = Hra and $o = --kGo^PoX^ . (B.7) 



The lineeir regime. Consider perturbations about the aforementioned background 
solution 

p^Po + 5p, p^po + 5p, = < + $ = $o + ■ (B.8) 

where Sp, Sp, 5v'^ and 5$ are the perturbed first order variables with spatial as well as 
temporal dependence, Sp — Sp{t, Xa). In the linear regime the perturbed quantities are 
much smaller than their zero order counterparts, Sp <ti po- During this period higher 
order terms, for example the product SpSva, are negligible. This means that different 
perturbative modes evolve independently and therefore can be treated separately. Note 
that Sva = Ua is simply the peculiar velocity describing deviations from the smooth 
Hubble expansion. Also, the fluid pressure is related to the density via the equation 
of state of the medium, which is p = p{p) for "barotropic" fluids (see §1.1). 

In perturbation analysis it is advantageous to employ dimensionless variables for 
flrst order quantities. Here, we will be using the dimensionless "density contrast" 
S = Sp/ Pq. Throughout the linear regime Sp <^ po meaning that S <^ 1. Substituting 
eq. (B.8) into (B.l) and keeping up to flrst order terms only, we obtain the relation for 
the linear evolution of density fluctuations 

% + -dJv''^0. (B.9) 
at a 

Also, substituting eq. (B.8) into (B.2) and linearizing we obtain the propagation equa- 
tion for the velocity perturbation 

3Si) iP' 1 

— ^ + HSv^ + -^dj + -dj^ = , (B.IO) 

at apo a 

where Vs is the adiabatic sound speed, defined, according to eq. (1.19), as 

'^V'^w^/\ (B.ll) 



s 
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where S denotes the entropy. Thus, using the considerations exposed in §1.1, we find 
that in a dust fiuid Vg — 0, while in a radiative fiuid Vg — l/-\/3. Note that eq. (B.ll) 
indicates that ty > 1 is impossible, because it would imply Vg > 1- Otherwise, iiw < 0, 
then it is no longer related to the sound speed, which would be imaginary. 

Finally, inserting eq. (B.8) into the Poisson equation and linearizing we arrive at 
the relation 

8^5^ - AnGa'^poS = (B.12) 

for the evolution of the perturbed gravitational potential. Results (B.9)-(B.12) deter- 
mine the behaviour of the density perturbation completely. When combined they lead 
to a second order differential equation that describes the linear evolution of the density 
contrast. In particular, taking the time derivative of (B.9) and employing eqs. (B.IO) 
and (B.12), we find to linear order 

The Jeans Length. Equation (B.13) is a wave-like equation with two extra terms 
in the right-hand side; one due to the expansion of the Universe and the other due to 
gravity. It is, therefore, natural to seek plane wave solutions of the form 

k 

where 5(k) = ^(k)(^) and k^ is the comoving wavevector. We obtain 

which determines the evolution of the k-th perturbative mode. The first term in the 
right-hand side of the above is due to the expansion and always suppresses the growth 
of (5(k). The second reflects the conflict between pressure support-^ and gravity. When 
AnGpo 3> v'lk'^/a^ gravity dominates. On the other hand, pressure support wins if 
vlk"^ /a^ > AttGpq. The threshold AirGpo = v^k"^ /a^ deflnes the length scale 




(B.16) 



The physical scale Aj, known as the Jeans length, constitutes a characteristic feature 

of the perturbation. It separates the gravitationally stable modes from the unstable 

^For a baryonic gas the pressure is the result of particle collisions, while for "dark matter" collisions 
are negligible and pressure comes from the readjustment of the orbits of the coUisionless species. In 
both cases it is the velocity dispersion of the perturbed component which determines the Jeans length. 
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ones. Fluctuations on scales well beyond Aj grow via gravitational instability, while 
modes with A <^ Aj are stabilized by pressure. The Jeans length corresponds to the 
Jeans mass, defined as the mass contained within a sphere of radius Aj/2 



Mt 




(B.17) 



where p is the density of the perturbed component. 

Multi-component fluids. When dealing with a multi-component medium 
(e.g. baryons, mirror baryons, photons, neutrinos or other exotic particles), the k-th 
perturbative mode in the non-relativistic component evolves according to 



d^ 



a"' 



:b.i8) 



where the index i refers to the component in question. The sum is over all species and 
= Pi/ Po provides a measure of each component's contribution to the total background 
density po — J2iPi- To first approximation, H is determined by the component that 
dominates gravitationally, while Vg is the velocity dispersion of the perturbed species 
which provide the pressure support. 

Solutions. Wc will now look for solutions to eqs. (B.15) and (B.18) in the following 
four different situations: (i) perturbations in the dominant non-relativistic component 
(baryonic or not) for t > tcq, (ii) fiuctuations in the non-baryonic matter for t < teq] 
(iii) baryonic perturbations in the presence of a dominant collisionless species. 
Perturbed Einstein-de Sitter Universe. Consider a dust dominated (i.e. p — 
= Wg) FRW cosmology with fiat spatial section. This model, also known as the 
Einstein-de Sitter Universe, is though to provide a good description of our Universe 
after recombination. To zero order a oc t^^^, H = 2/3t and po = l/67rG't^ (see §A.2). 
Perturbing this background, we look at scales well below the Hubble radius, where the 
Newtonian treatment is applicable. Using definition (B.16) and the relation A = 27ra/k, 
eq. (B.15) becomes ^ 



2d2(5 4 d(5 2 
I — f 



1- ^ 




6 = 0. 



:b.19) 



For modes well within the horizon but still much larger than the Jeans length (Aj <^ 
A ^ Xh), we find 

S = Cit^/^ + C2t-^ (B.20) 
^Note that from now on we will drop the tilde and the wavenumber (k) for convenience. 
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for the evolution of the density contrast. As expected, there are two solutions: one 
growing and one decaying. Any given perturbation is expressed as a linear combina- 
tion of the two modes. At late times, however, only the growing mode is important. 
Therefore, after matter-radiation equality perturbations in the non-relativistic com- 
ponent grow proportionally to the scale factor (recall that a oc i^/^ for dust). Note 
that baryonic perturbations cannot grow until matter has decoupled from radiation 
at recombination (we always assume that teq < tree)- Dark matter particles, on the 
other hand, are already coUisionless and fluctuations in their density can grow immedi- 
ately after equipartition. After recombination, perturbations in the baryons also grow 
proportionally to the scale factor. 

On scales well below the Jeans length (i.e. for A <^ Aj), eq. (B.19) admits the 
solution 

which describes a damped oscillation. Thus, small-scale perturbations in the non- 
relativistic matter are suppressed by pressure. Also, before recombination baryons 
and photons are tightly coupled and {v'^)h oc T5 — T't cx: a~^, after decoupling {v'^)h oc 
Tb oc a-^. 

Mixture of radiation and dcirk matter. Consider the radiation dominated regime 
when a oc i^/^ and H — l/2t. On scales much smaller than the Hubble radius we can 
employ the Newtonian theory to study perturbations in the non-relativistic matter. 
Applying eq. (B.18) to a mixture of radiation and coUisionless particles {vs — 0) prior 
to equipartition (pdm <S p-y — Po), given that the small-scale photon distribution is 
smooth (i.e. {5j) ~ 0), we have 

.d^^DM^d^^Q^ (B.22) 



dt2 dt 



that admits the solution 



Sdm = Ci + Cslnt . (B.23) 

Thus, in the radiation epoch small scale perturbations in the coUisionless component 
experience a very slow logarithmic growth, even when A > Aj. The stagnation or 
freezing-in of matter perturbations prior to equilibrium is generic to models with a 
period of rapid expansion dominated by relativistic particles and is sometimes referred 
to as the "Meszaros effect" [122]. 

Mixture of dcO^k matter and baryons. During the period from equilibrium to 
recombination perturbations in the dark component grow by a factor of Qrec/cieq — 
Teq/Trec — 21 fl^h^ . At the same time baryonic fluctuations do not experience any 
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growth because of the tight couphng between photons and baryons. After decouphng, 
perturbations in the ordinary matter will also start growing driven by the gravitational 
potential of the coUisionless species. Consider the post-recombination Universe domi- 
nated by non-baryonic dark matter. For baryonic fluctuations on scales larger that Aj, 
eq. (B.18) gives 



where C is a constant. The initial conditions at recombination arc 6\y = 0, because 
of the tight coupling between the baryons and the smoothly distributed photons, and 
Sbm 7^ given that the dark matter particles are already coUisionless. The solution 



shows that 5^ — > 5dm as a ^ Orec- In other words, baryonic fluctuations quickly catch 
up with perturbations in the dark matter component after decoupling. Alternatively, 
one might say that the baryons fall into the "potential wells" created by the coUisionless 
species. 

B.2 Linear relativist ic theory 

The relativistic equations. The basis of the relativistic analysis is the Einstein field 
equations (A. 5), describing the interaction between matter and spacctimc geometry 
(see appendix A.l). The Einstein tensor Gais has the extremely important property 
of having an identically vanishing divergence, that is V^Gas = 0. When applied to 
eq. (A. 5) neglecting the cosmological constant, the latter yields the conservation law 



For a perfect fluid the stress-energy tensor takes the simple form (A. 8). Here, similarly 
to the Newtonian analysis, we assume a barotropic fluid with p — p{p). The relativistic 
analogues of the continuity and Euler equations are obtained from the timelike and 
spacelike parts of the conservation law (B.26). In particular, we obtain the energy 
density conservation law 




(B.24) 




(B.25) 



V^T^fS = . 



(B.26) 



xp + 3H{p + p) = , 



(B.27) 



and the momentum density conservation equation 



{p + p)ua + 'DaP^O 



(B.28) 
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Equations (B.27) and (B.28) are supplemented by 

Rabu'^u'' = ^k{p + 3p) , (B.29) 

which is the relativistic analogue of the Poisson equation and relates the spacetime 
geometry to the matter sources. Note that the cosmological constant has been set to 
zero. 

The linear regime. Perturbing eq. (B.27) and defining 6 = Sp/po, the perturbed 
continuity equation to first order gives 

S' -3wHoS + 3{l + w)5H ^0 , (B.30) 

where the dash indicates differentiation with respect to t and w = po/po determines 
the equation of state of the medium. Also, starting from eq. (B.28) we have 

4 

{SH)' + 2HodH + -TiGpod + ' D^d = (B.31) 

to first order, where 6H describes scalar deviations from the smooth background ex- 
pansion. Using eqs. (B.30) and (B.31), and Fourier decomposing, in terms of the scale 
factor a we arrive at 



[l + 8w- - 6v^) - - 



3 



5(k) = (B.32) 



for the evolution of the k-th perturbative mode. 

Solutions. We will now seek solutions to the relativistic perturbation equations to 
supplement the Newtonian results of the previous section. The cases to be considered 
are: (i) super-horizon sized perturbations in the dominant non-relativistic component; 
(ii) fluctuations in the relativistic matter before matter-radiation equality both inside 
and outside the Hubble radius. 

Perturbed Einstein - de Sitter Universe. On scales beyond Ah one needs to engage 
general relativity even when dealing with non-relativistic matter. For pressureless dust 
w — — and for modes lying beyond the Hubble radius (A » Ah and \?vl/a'^Hl <^ 
1), eq. (B.32) becomes ^ 

^d?5 3 d5 3, 
a^— + -a— --5 = 0, (B.33) 
da2 2 da 2 ' ^ ' 

with 

(5 = Cia + C2a-^/2 . (B.34) 



^Again we have dropped the tilde and the wavenumber for simpHcity. 
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Appendix B. Structure formation 



Thus, large-scale perturbations in the non-relativistic component grow as 5b oc Q oc i^/^. 
The radiation dominated era. Before equahty radiation dominates the energy 
density of the Universe and w — 1/3 = v^. During this period eq. (B.32) gives 

On large scales, when A ^ Ah and k'^/a^H^ -C 1, the above reduces to 

-25^ = 0, (B.36) 

with a power law solution of the form 

= CiO^ + Caa"^ . (B.37) 

Hence, before matter-radiation equality large-scale perturbations in the radiative fluid 
grow as 6y oc a^. Note that cq. (B.36) also governs the evolution of the non-rclativistic 
component (baryonic or not), since it does not incorporate any pressure effect. There- 
fore, solution (B.37) also applies to baryons and collisionless matter. 

On sub-horizon scales, with A <^ Ah and k^/a^i^Q ^ 1, eq. (B.35) becomes 

and admits the oscillatory solution 

= Ce^^H/^ . (B.39) 

Thus, in the radiation era small-scale fluctuations in the relativistic component oscillate 
like sound waves. Given that Ah /A ^ 1, the oscillation frequency is very high. As a 
result, (6^) ^ on scales well below the Hubble radius. In other words, the radiative 
fluid is expected to have a smooth distribution on small scales. 

Note that in the radiation era the transition from growing to oscillatory modes 
occurs at A ~ Ah, which implies that before equality the role of the Jeans length is 
played by the Hubble radius. 
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